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Beniitzung von Singularititen bei der Berechnung magnetischer 
Felder mit Hilfe des Differenzenverfahrens 


Von Miroslav Kotal, Prag 
Mit 8 Textabbildungen 


Zusammenfassung. Magnetische Felder, die partiellen Differentialgleichungen mit verinder- 
lichen” Koeffizienten geniigen, kénnen leicht mit Hilfe des Differenzenverfahrens (Relaxations- 
methode) behandelt werden. Zur Berechnung der Luftgebiete mit ferromagnetischen Elementen 
benutzt man das Rechenverfahren der singuliren Punkte, wie an einem Beispiel vorgefithrt wird. 


I. Einleitung 


Die Untersuchung von magnetischen Feldern in Luftgebieten und in ferromagneti- 
schen Gebieten mit veranderlicher Permeabilitaét ist beim Entwurf einer elektrischen 
Maschine eine der wichtigsten Aufgaben des Berechnungsingenieurs. Die Ermittlung 
dieser Felder stot im Falle veranderlicher Permeabilitét auf groBe Schwierigkeiten und 
die bisher bekannte Lehmannsche graphische Methode! ergibt nur sehr angenaherte 
Ergebnisse. Die praktische Bedeutung dieser Methode wird auch durch die Tatsache 
vermindert, dafi numerische Werte durch Interpolation gewonnen werden miissen. 
Die Beibehaltung einer konstanten Permeabilitat wiirde zwar eine gewisse Vereinfachung 
der Berechnung bringen, aber die gewonnenen Werte wiirden nicht den Eigenschaften 
des ferromagnetischen Materials entsprechen. 

Bei der Berechnung von Feldern mit veranderlicher Permeabilitat mu8 eine direkte 
Relation zwischen FluBfunktion und Permeabilitat hergestellt werden, und dies erschwert 
die Rechenarbeit. Das gilt sowohl fiir die Lehmannsche Methode als auch fiir die 
Naherungsmethode mit Hilfe des Differenzenverfahrens, welche in der Literatur als 
Relaxationsmethode bekannt ist und in den nachsten Abschnitten der vorliegenden 
Arbeit behandelt wird. Die Grundziige dieses Relaxationsverfahrens wurden in einer 
friiheren Abhandlung besprochen?. 

Weitere Komplikationen entstehen auch bei der Einfiithrung der Randbedingungen 
an der Grenze Eisen—Luft, wo Refraktion der magnetischen Linien die Kinfiihrung 
singularer Punkte in die Berechnung erfordert, wie spater gezeigt wird. Das Verfahren 
mit Hilfe der singularen Punkte, das in anderen Gebieten der technischen Wissen- 
schaften, wie z. B. Aerodynamik und Festigkeitslehre, mehr bekannt ist, findet auch 
in der Elektrotechnik seine Anwendung. Auf diese Weise kénnen die Probleme der 
Potentialfelder, die bisher als unlésbar betrachtet werden, mit Hilfe der Relaxations- 
methode gelost werden. 


II. Luftgebiete 


Bei der Berechnung von Luftfeldern wird die Laplacesche bzw. Poissonsche Diffe- 
rentialgleichung 


1 Th. Lehmann: Spectres ferromagnétiques. Revue Générale de V’Electricité XIX, 1926, 


S. 43, 85. 4 
2 M. Kotal: Relaxationsmethode. Osterr. Ing.-Archiv X1, 1957, 58. 93. 
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ae te Y) (1) 


naherungsweise durch eine Differenzengleichung 


ap 
Ou? 5 


2 2 2 : 2 Xp 
a, (a, + As) oa Gg (Ay + Ay) Mie a3 (as + ay) oe Ay (44 + Ap) Pa 
1 1 
—2(5 ds ae = hf (x, y)o (2) 


ri 


Z 
| ersetzt. Diese Differenzengleichung gilt fiir ein 
azht hy unregelmaBiges Relaxationsbild (Relaxations- 

kreuz, Abb. 1). Wenn ein regelmaBiges Relaxa- 


a3lv aly tionsbild Ay = Ag = Az = A, = M angenommen 
Jd . .- 
| if ; | wird, erhalt man 


ayly Git P2t+ Gxt Pa 4p, =—h' f (x, Y)o- (3) 


Die oben angefithrten Beziehungen (2) und 
(3) benutzt man auch bei der Berechnung der 
ferromagnetischen Gebiete, allerdings unter der 
Abb. 1. UnregelmaBiges Relaxationskreuz Voraussetzung einer konstanten Permeabilitat. 


Ill. Ferromagnetische Gebiete 


Bei der Ableitung allgemeiner mathematischer Beziehungen geht man von der 
Poissonschen Differentialgleichung mit verainderlichen Koeffizienten 


og (A5e) tay (Bog l =P @ yw) (4) 
aus. Wenn A = B=1 gesetzt wird, bekommt man die Poissonsche Gleichung in ihrer 
gewohnlichen Gestalt, und fiir f (x, y) = 0 entsteht die Laplacesche Gleichung. Fiir 
A, B und f (x, y) werden meist bekannte Funktionen der x, y genommen, was leider 
bei der Berechnung der elektrotechnischen Felder nicht méglich ist*. 

Die mathematische Formulierung der magnetischen Felder in der Elektrotechnik 
kann in Form einer Laplaceschen Differentialgleichung 


cd] _, OP ] _4 &Y 
me =e (5) 


gegeben werden, worin u die Permeabilitat, g der magnetische Induktionsflu8 ist. 


Mit Hilfe eines Kinheitsdurchschnittes wird fiir unsere weiteren Betrachtungen der 
InduktionsfluB gm aus der magnetischen Induktion B abgeleitet werden. Weiterhin 
wird in Ubereinstimmung mit den Lehmannschen Empfehlungen das dreidimensionale 
Problem eines magnetischen Feldes einer elektrischen Maschine auf ein zweidimensiona- 
les reduziert. 

Nach Ausfiihrung der Differentiation in Gl. 5 und mit 


dinuw dz 1 
ee Th Sih (6) 


erhalt man folgende Differenzengleichung fiir ein unregelmafiges Relaxationskreuz 


* H.Emmons: The Numerical Solution of Partial Differential Equations. Quart. Appl. 
Math. 1, 1944, 8S. 173. 
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2 . 2 de 2 2 ig, ; 2 
a, (a, + a3) id Ay (Ay + Ay) Pe sz (ds + 4) Ps + My (Aq + Gy) ate Ge Gg | a, 4 - 
i Ag Oe ek 1 ay 1 : 
Z (@, +43) fy ' dg ay " Ho ag (43 + a) » Us 


. ( Cee fe Cae Oy a 
q, (a, + a5) 717 a,a, 49 As (Az + a4) Ps 


1 


( Ay 1 | a&—Y% 1 Ay ) 
— n 4 
Ag (Gz + Ay) [Me Ay Ay Mo My (Ag + Ay) fg 


: Xs Ay — A> 
[ (az + 44) Par Ug Uy ae Aq (M4 + Ag) vr) ee (7) 


und fiir ein regelmaBiges Relaxationskreuz 


1) u Ly bs Ly fy Ly pe’ 
(1 int), : (1 reWake | (1 4 ni) Ps : (1 | gn) Ps — 49 0= o- (8) 


In diesen Gleichungen stellt Ry ein sogenanntes Residuum dar. 

Ks ist vorteilhaft die bei y,...g, stehenden Koeffizienten mit Hilfe eines Nomo- 
gramms mit zwei aneinandergehefteten Doppelleitern fiir magnetische Induktion und 
Permeabilitaét zu losen. Auf diese Weise kann man unmittelbar aus den g-Werten die 
Korrektionswerte ableiten. Die Nomogramme ebenso wie auch die spater erwahnten 
Diagramme, die im praktischen Beispiel in Abschnitt VII benutzt worden sind, wurden 
wegen Raummangels weggelassen. 


IV. Randgebiete mit Beriicksichtigung der Grenzbedingungen 


Die Berechnung eines kombinierten Eisen-Luft-Gebietes und die gleichzeitige 
Beriicksichtigung der veranderlichen Permeabilitat erfordert zunachst eine Losung der 
zweiten Randwertaufgabe (Neumannsches Problem). Dann mu auch die Refraktion 
der magnetischen Induktionslinien betrachtet werden. Dieses Verfahren kann als viertes 
Randwertproblem der Potentialtheorie bezeichnet werden. In den bisher veroffentlichten 
Arbeiten werden diese Bedingungen nur unzu- 
reichend beriicksichtigt. Inden Southwellschen 
Arbeiten* wird nur konstante Permeabilitét in 
Betracht gezogen und keine Beziehung zwischen 
veranderlicher Permeabilitat und Refraktion der 
magnetischen Linien eingefiihrt. Auch Govor- 
kov® legt solche Relationen nur teilweise zu- 
grunde. 

Eine genaue analytische Fassung der oben 
angefiihrten magnetischen Erscheinungen geben 
die folgenden Gleichungen 

Bet ok, CO" ap! 4 ap” 
[ij Ce Wis CI: tii, et? 
worin yu, yp’ Permeabilitat bzw. Funktionswert 

in den Luftgebieten, 
u,, p Permeabilitat bzw. Funktionswert 
in den ferromagnetischenGebieten, 


vy Normale zur Randkante Abb. 2. Unregelmafiges Relaxationskreuz 
t Tangente zur Randkante. an einer Luft-Eisen-Grenze 


(9, 10) 


u=7 es 


4 R.V. Southwell: Relaxations Methods in Theoretical Physics, Oxford 1946. 
5 V. A. Govorkov: Theorie des magnetischen Feldes, Moskau 1957 (in russischer Sprache). 
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Wir betrachten weiter nach Abb. 2 das Relaxationskreuz an der Luftseite und legen 
die in Differenzenform ausgedriickten Randbedingungen Gl. 10 und 9 zugrunde. Dann 
folgen fiir die Werte p,’ und ¢,’ eines fiktiven Luftfeldes im Eisen die Ausdriicke® 


P i re _,\ #, 8in B cos w + p,, cos B sin « My Pa Po (2 sin a cos B , 
ae & ae ae ty, Sin (« + B) Tos See, My sin(a +B) | %° 
(11 
Lay LI, at 
, Mn ” rr\ Sin % COS B ae My sin B COs % 3 te ‘S ( \ sin « cos B , 
< = (fe %. ae i, Sin (« + B) aks a by hecgea = Me 
12) 


Dabei muB betont werden, daB die Werte y,’ an der Luftseite und y,’’ an der ferro- 
magnetischen Seite in unserem Falle nicht von gleicher GréBe sind. Diese Tatsache 
wird durch Refraktion der magnetischen Linien verursacht. 

Die auf diese Weise gewonnenen Werte ¢,’ und gy,’ benutzt man fiir die Berechnun- 
gen, die mit Hilfe der unregelmaBigen Relaxationsbilder durchgefiihrt werden. Kine 


kraftige Vereinfachung erreicht man mit «+ p= a Weiters wird «= 0 vorausgesetzt, 


d. h. die Luft-Eisen-Grenzlinie fallt mit dem Relaxationsstrahl 2—0—4 zusammen. 
Die Gleichungen fiir die Unbekannten ¢y,’ und gy,’ nehmen dann folgende Gestalt an: 


, , Hy, / / 
Ps = lL, (vs AeA ) r Po (13) 
Pr ” a? Sok a” uv a,” aa a” i a,” a a,” if 
Paces ae ee ae 0 * gave Ge are iat as (14) 


Diese Gleichungen kénnen als Sonderfalle der Gl. 11 und 12 betrachtet werden. 
Nach Einsetzen in die Gl. 2 bekommt man 


2, ng 
ay (Ay + As) eden 


2 | 2 (= y i 2 i 
My (Ay + Ay) Az (Az + M4) \My eee Ay (Ap + Ay) Pe a2 


| 2 P rr By aE 2 (<)' ” 2 fan | 2 (Ay a y) mr 
' Az (44+ a3) 73 bg Ay (Ay + My) \y, Pa dg (dg +) By, ~~ Ay Ay? Po 


2 2 2 (dy + G4) ’ 
E (a; + a) As Ug | Ay Ay” | coe Ko CS, 
und fiir ein regelmaiBiges Relaxationskreuz 
G1 £Ye FP. + 294 + ms = me P — 399 = Ry. (16) 
Die Zeichen + oder —, die schon in Gl. 9 und 10 auftraten, ergeben sich aus der 


Richtung des Feldes. Regelmafig ist y,'’ — p,"’ positiv fiir negatives gy,’ — my,’ und 
umgekehrt. Diese Tatsache liefert einen gewissen Zusammenhang zwischen den Werten 
des Luftgebietes p,’ und ~,’ und den Werten des ferromagnetischen Gebietes q,’’und 
pa. Fir diese Betrachtungen ist auch der Drehungssinn der Numerierung der Re- 
laxationsstrahlen maBgebend. 

Auf ahnliche Weise werden auch die unbekannten Werte ¢,"’ und y,”’ fiir das Relaxa- 
tionskreuz an der ferromagnetischen Seite bestimmt 


ee ee ay bn sinacosf + u,cosasinB ay ; on (Bo sin B cos B YY 

Po = (P2 Po) mt, Sin (« + B) ay (Y1' — Po )(@ 1) sin (« ++ B) +o" (17) 
ae ; 1, 4, Sin «cos B+ w,,cosasinB | a ‘ atin sin « COs % ms 

Pi = (% Po) m, Sin (a + B) Ay (Y2' — Po ) (Be \am («+ B) +o" (18) 


oder, wenn die oben eingefiihrten Vereinfachungen zugrunde gelegt werden, 


° M. Kotal: Methoden der Einheitsrelaxation. Strojirenstvi 9, 1959, S. 296 (in tschechischer 
Sprache). 
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rr My , 


ge ate 9 (G1 — Go) + Po" (19) 


as? — a? Cn? — 2 


{i , | io ool 2 , = ey 2 a \e vr 
Pata a az 0 = (#) A alae a a +( *) V4 is (20) 


Uy 


Wenn diese Ausdriicke in die Gl. 7 eingesetzt werden, bekommt man 


ae LS | 2 | 2 (3) | ee eli 2 Ba} HL, 2 / 
3 (4, +43)73 " la, (a,+a,) ' a, (a,+a,) \a,) |?4 @, (aay) P21 Ib, @, (% +a, % = 


ae 2 (@2)" Bees 2 Me 2 (dy 3 a) / 
Oy (Ag + Ag) \Ay Pa a, (a, + as) u,, My A? DO 
1 1 

) 5) 9 In Fe le Sia yt aa 
CS (dy as “| eas R Ls My P33 — Po 
As (a, + ag) Ay Ay Gig Gy Gore eG As As 

Oy 1 Ay, — a 1 a 1 

In — + +—ln- 2 

& (dg+ 44) “bg | My ay Mo My (My + Gy) Fa M4 

Ag wr Az — Uy vr A v7 

‘ el E..| | hse Ee 9 

: (a, + a) Cals A, A Po My (M+ Gg) 4 (21) 


und fiir ein regelmaBiges Relaxationskreuz 


" fa rr (< 1) oe Ts a ee ' / 
P3 (1+ in“) +9, (2 z inet) — 29, Inv. nu, 1° F Ye t%_ 4 ei 


; Mo) on < 
=o (3 + In2) po" = Bo (22) 

Die Zeichen + oder — sind wieder durch die Feldrichtung gegeben: p,' — 4’ wird 
positiv genommen, wenn @;'’ — ¢,’’ negativ ist und umgekehrt. 


V. Kreisférmige Gebiete 


Kreisformige Gebiete kommen bei der Untersuchung magnetischer Felder in elek- 
trischen Maschinen, besonders bei den sogenannten Querschnitten sehr haufig vor. 
In diesem Falle bietet die Benutzung der Zylinderkoordinaten manche Vorteile. Die 
Laplacesche bzw. Poissonsche Differentialgleichung, welche die Verteilung einer Stré- 
mungsfunktion darstellt, lautet unter Voraussetzung einer veranderlichen Permeabili- 
tat in diesen Koordinaten 


i ie ole @1 1 @& 9 
r le ub a t 39 wr ad = f(*, y)- (23) 


\ 


Wenn die Kriimmung der Relaxationsstrahlen im Relaxationskreuz vernachlassig- 
bar ist und diese daher durch Linienabschnitte ersetzt werden konnen, erhalt man nach 
einfacher Umformung aus Gl. 23 folgende Differenzengleichung (Abb. 3) 


2 (2 fa) i Pa fy Ay *) Po (2 Uy *) 
Gat Oy\ 0," a, (iy (Ay + Ay) ( T) | | Qy (Ao + a) (Pay ee 
wee [eee h io | 2 oie 
lay re: as) : : 1am ‘: : As ee Ay) ;. 2 
aes, Ciara, Fa)” 
: ie me M4) a : y : = Mo Ag ee a) a 
ashe es fo a laayt|t Be 


wobei 7,0; = ah gesetzt wurde. 
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Wenn die Luft-Eisen-Grenze mit dem Relaxa- 
tionsstrahl 2—0—4 zusammenfallt, bekommt man 
folgende Formel fiir das Relaxationskreuz an der 
Luftseite 


oe My 
; 9 Ms | 2 = h | aa h (2) re 
a (4, + as) 1 Gy (dy + Ay) | Gy (My + Gy) \Oy Pa 
Os 
a "o i ” 2 rr 
ly (Ay + Gy) Pe 3 (a, +3) 7? py 
Ay 


in) (@: ‘ Un 
Gy (Ag + Ay) \M Pa 
2 UH, A, — Ag (2 “*h)]— 


Po (dz + 4) M3 M,, Ag” Ay To 
| Se eee, 2 | 2 (a, — “| = 
rie Ay My | Ay (dy + az)” Ay A, a3 
Abb. 3. Relaxationskreuz fir kreis- 
formige Gebiete cme O (25) 


VI. Rechenverfahren mit Hilfe der singuliren Punkte in den Grenzgebieten 
Kisen—Luft 


In den vorangehenden Abschnitten wurden die Grenzgebiete der magnetischen 
Luftfelder und der ferromagnetischen Felder mit veranderlicher Permeabilitat, die 
voneinander durch gerade Grenzlinien getrennt werden, behandelt. Das Verfahren 
kann auch im Falle maBig abgerundeter Kanten beibehalten werden. Gewisse Schwie- 
rigkeiten entstehen an den Stellen, wo die Kanten in spitzen Winkeln gebrochen wer- 
den. Die Funktionswerte dieser Winkelscheitel diirfen nicht in regelmaRigen Relaxa- 
tionsverfahren verwendet werden. 
In diesen Punkten erreichen nam- 


Toe m2 “¥ 
< / 3 lich die Funktionswerte unend- 
\ P ue liche oder unbestimmte Groen 
“N 6 70 1 7 und miissen daher als singulire 
‘f i a g 


Punkte angesehen werden. Das 
magnetische Feld mu in der 
Nahe dieser Punkte durch beson- 
dere Funktionen dargestellt wer- 
den. Wegen der erwahnten physi- 
kalischen Erscheinungen an der 
Luft-Eisen-Grenze benutzt man 
das Rechenverfahren mit Hilfe 
der singularen Punkte sowohl an 
der Luftseite als auch an der 


Tv =7 


2 


Raa 3 ts +5 Send. 3 ferromagnetischen Seite des be- 
wie : a trachteten Gebietes. 
i met me f. Die Lésung des Problems 


k6nnte auch so durchgefiihrt wer- 
Abb. 4. Quadratisches Netz mit Singularitét im Punkt 0 den, daB in der nachsten Um- 
gebung des singuliren Punktes 
ein verdichtetes Netz benutzt wird, allerdings ohne Hinzunahme dieses Punktes in 
das Rechenverfahren. Diese Methode wird durch die Vermehrung der Anzahl der 
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Netzpunkte sehr miihsam und deshalb nur selten verwendet ( 
lehre). Das Singularitaétenverfahren bietet dagegen manche Vorteile. Es wurde von 


z. B. in der Festigkeits- 


a 
Zu 

Ze 25 eA 
PSN ey) i : 
Sc0s() pees} 
Zsint } Zi 

ae GOS: 
cos f } i “4 2h O} 
2S/1 (mane 
22 } ost Uf 42 
ne BY sinim E(Zn+1)p 
gm # 

Vie aes 

le 10 +— G/ 

a5 S f{—7| | 

t lie 
CRO >) 
4 G65 y 
7G9F OI GEE WS OTE G7 GB G6 W? SS V2 71 
160° Zor 210° IG? 360° 


Tod 
cosfm@(én+1/p 


IS 


J 
gos me (Zn+iyp 


Abb. 5. Hilfskurven fiir Gl. 28 


L. C. Woods? benutzt und wird fiir unsere Zwecke hier abgedndert. In der Nahe des 
singularen Punktes wahlen wir als Ersatzfunktion die Funktion re’, worin r und 6 


7 L.C. Woods: The Relaxation Treatment of Singular Points in Poisson’s Equation. Quart. 
Journ. Mech. and Appl. Math. VI, 1953, S. 163. 
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Polarkoordinaten sind. Mit Hilfe dieser Funktion, deren Real- und Imaginarteil har- 
monisch sind, bilden wir Beziehungen, die erwiinschte Eigenschaften haben. Wir 
kénnen also nach Abb. 4 fiir Luftgebiete 


(rem — (c sin m 6 + d cos m 8) (5-)" (26) 
und fiir Gebiete mit ferromagnetischem Eisen 
m In (=) (27) 


(die sogenannte logarithmische Permeabilitat) schreiben. Die Konstanten c und d 
werden mit Hilfe folgender Gleichungen 


+ Get os+G1—4M = f(x, Y)o ta ein (mn 5) 


i (2n +1)) 
s+ Pe + Gr t+ % — 4% = (2h)? (2, Yo +> sin {m= ae Ny 
sk 
1 3 ) _ 


bestimmt. Die Beiwerte der goniometrischen Funktionen dieser Gleichungen wurden 
fiir einzelne Werte von m in Abb. 5 eingetragen. Die Konstanten c und d kénnen aller- 
dings nur fiir eine bestimmte Konfiguration der g-Werte (einen Relaxationsschritt) 
benutzt werden. Andere Konfigurationen der g-Werte erfordern auch eine neue Bestim- 
mung der Konstanten c und d. m ist em Exponent, dessen GréBe vom Winkel 6 abhan- 
gig ist. 

Die Art der Stromung nach Gl. 26 hangt von der GroBe des Exponenten m ab. Diese 
Singularitat ist identisch mit einem bereits friiher gelésten Singularitatsproblem des 
Stirnfeldes emes Synchrongenerators?. 


Geht man von der Gl. 26 aus, so bekommt man z. B. fiir das Residuum des Netz- 
punktes 5 folgenden Ausdruck 


Rs = Gis + Pia + V2 + Pi 4 ys; h? (2. Ys — 


r / 
= c {5% |sin m(5 —) +sinm | + sin ms —4- es 


—a{5%"| cosm (5 — 6) + cosm o| + cosm + 1— 4:2%"cogsm- 


= Vis + Pia + 2 + G1 — 495 — h? f (x, y)s — c sin; — d cos,, (29) 
worin 6 = arctg (0,5). 


Ahnliche Gleichungen erhalt man auch fiir die iibrigen Netzpunkte. Die Beiwerte 
der zustaéndigen trigonometrischen Funktionen kénnen zur Erleichterung der Berech- 
nungen in Diagrammen eingetragen werden. 


Fir die Behandlung der Gebiete mit veranderlicher Permeabilitat werden diese 
Gleichungen nach Gl. 8 abgeaindert. Man bekommt z. B. fiir das Residuum des Netz- 
punktes 5 folgenden Ausdruck 


2'g..09. 2, Anm.2 
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1 i! 
esa zine) pra(I oy ae: al ps(14 pint) 
: 2 


ta( + Int) — 4 9,— 2? f (x, 9) — 


. 4 7 - : 5 J 
—€ [5% sin m(F —0) + 52" sinmoé + sinm>— 4:2%™ sin mt | — 


+ 


bo! 


— gd \p%n ae Vy m i Y, z 
cosm zie cos m 6 + cos m=z + 1 — 4-24"cosm 7 — 


- 


c Hag |_- Tt 1 : c be s : Tt 
saya pe So Ba aay ake Ye m i 
q im - ee ms OAL SNOWY 6| 4 In ie 5 sin (3 6)| 
Toe 0 I d l 
In“? Ym Mia Vm ey 
qin |cosm 3 — 5*2™ cos m 6| Z In re D302 C68 \*5 6}|. (30) 


Wie schon erwaihnt wurde, ist es vorteilhaft, die trigonometrischen Funktionen im vor- 
aus zu errechnen. 


Bei der Ausrechnung der Residuen in den Randpunkten miissen Gl. 9 und 10, die 
fiir diese Punkte maBgebend sind, erfiillt werden. Nach Differentiation von Gl. 26 
bekommt man, wenn gleichzeitig m0 = 0, ; = 0} 5 7 1, h=1 gesetzt und 
Gl. 16 beachtet wird, folgende Formel fiir das Residuum eines Relaxationskreuzes an 
der Luftseite 


by, / 4) 3 EF 4 / i / , tr t 
k,=m (c+ d)+ 9, AO, = Og Mea 2 0, ee Aa Sag ea 


My, Li 


as Dm at ao DV m Mer j a AN — 
c|2 sin m7 + 2 Hsin m( Z 


—d 2%" cos ma — 9m 1 9. 9m1 Q%ml cog m(— 4) — 3}. (31) 
M3 4 
In diesen Gleichungen nehmen die g-Werte des singularen Netzpunktes 2 nicht- 
singulare Werte an. Wenn schon dieser Punkt singulare Werte besitzt, benutzt man in 
der Gl. 31 anstatt y,’’ — m,"" den Ausdruck 2 (gy) — 4’’). 


Der oben beschriebene Rechenvorgang wird fiir den Fall des Luftgebietes mit ferro- 
magnetischem Hisen und fiir die Bedingung 0,5 <m <1 verwendet. Er gilt nicht fiir 
den Sonderfall im ferromagnetischen Eisen, wo m > 1 ist. Dieser erfordert das Rechen- 
verfahren auf Grund der logarithmischen Singularitat nach Gl. 27. Wir ersetzen zuerst 
den Wert in dem singularen Netzpunkt 0 durch einen nichtsingularen Wert gor. Das 
Residuum in diesem Punkte ist dann 


Ryo = G1 + G2 + Ys + Pa — Por — WF (H, Y)o (32) 
und die Residuen in den tibrigen Punkten sind 
Rs = is + Mia + 2 + %1 — 495 — A? f (a, y), — bm (In 5 — In 4). (33) 


In diesen Gleichungen ist & ein Faktor, dessen Bedeutung spater erklart wird. Die 
Funktionswerte in den Netzpunkten auBerhalb des betrachteten Gebietes werden 
durch Extrapolation bestimmt. 

Die auf solche Weise abgeleiteten Gleichungen werden jetzt nach den fiir Hisen- 
gebiete geltenden Bedingungen (Gl. 8) abgedndert. Man erhalt also 


ky= wl —4in@)+ pa(l ~~ )+ Ps (1 + ¢mn)+ ps(1 +4) — 4H, (34) 


und in ahnlicher Weise fiir die iibriggebliebenen Netzpunkte 
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By = oa (1— ging) + u(t ging) +oa(t+ ginge) +m(1+ gles) 
7 oe La Haat) (35) 
— 49, — m|Ing pe . 


1) uy lsy 1) 4s ( Sti \ = 
Ry =95(1— yn) + 10/1 ne a) vs(1 + gine) + gor aa ol hay) 
— 49, —_m[2In2— Fin 2 In] - &. (36) 


Die logarithmischen Ausdriicke kénnen ahnlich wie friiher vorteilhaft mit Hilfe 
der Nomogramme errechnet werden. 

Die Werte in den Randnetzpunkten werden durch einen ahnlichen Rechenvorgang 
wie im vorangehenden Falle bestimmt. Schreibt man zuerst eine Ableitung in der y- 
Richtung i 


3 1 
h(3) =S00—2¢.+eo0— MF (a Wo (37) 
bekommt man nach Gl. 27 
1 ay | 

por =e [402 — Pro + 2H f (aw, yo + 2h (52) + bm (— In 2h}. (38) 

Ahbnlich auch in der x-Richtung 

1 op 

por = [4P. — M1 + 2K E(w, yo + 2h (52) — km In 2h. (39) 
3 . Numerische Werte, die nach Gl. 38 


und 39 ermittelt werden, diirfen nur 
kleine numerische Differenzen auf- 
weisen. Bei gréBeren Abweichungen 
| wird ein Mittelwert genommen. Fiir 


: 0 é 
die Werte der Ableitungen a und we 


benutzt man die entsprechenden Aus- 
driicke, die schon fiir die Bedingungen 
an der ferromagnetischen Randseite 
eingefiihrt wurden. 
Da die Werte von go, schon bekannt 
| sind, kénnen auch die tibriggebliebenen 
z tee 76 Werte in den Randnetzpunkten an der 
ferromagnetischen Eisenseite bestimmt 
werden. Geht man von der Gl. 22 aus, 
so bekommt man 


Tel 2 wey ; 
SEMENVELNEITINS 08 LEWES. 


Abb. 6. Abhangigkeit des k, vom Seitenverhiltnis 
eines elektrischen Leiters 


ut Ho rreg 1 Ms 1 (a ee) LU, re , , 1 by 
P3 (1 +n) +9 (2 zn) — to Inu, 7° F Pa + Po + Py nos 
Ll ” ¢ 1) 4 1 ll 
—(3 +n), —km|(2+7mn@)n2+5In2in@ =F. (40) 


Der schon frither eingefiihrte Faktor k = k, k, ist in erster Linie von dem sogenann- 
ten Ahnlichkeitskoeffizienten , und in zweiter Linie von der geometrischen Form des 
betreffenden Gebietes k, (Abb. 6) abhangig. 

In dem anschlieBend betrachteten Beispiel wird der Wirkungskreis des singuliren 
Punktes in der GréBe von 2 bis 34 angenommen. An dieser Stelle betragt der Mittel- 
wert des magnetischen Feldes rund 15000 Gau8*. Wenn dort folgende Randbedingungen 


* In diesem Aufsatz werden die Berechnungen in der Einheit G fir die Induktion aufgestellt, 
da dieses Rechenverfahren manche Vorteile bietet. 
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p — 15000 log + = 0 (41) 


erfiillt werden, bekommt man nach Gl. 27 den. Wert des Faktors k,: 
— 15000 log,, e = — 6514,3 ~ — 6514. (42) 


Der Faktor k, kommt in diesem Beispiel nicht in Betracht. 
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Abb. 7. Beispiel der Losung eines Gebietes mit ferromagnetischem Eisen und Luft 


VII. Beispiel 


Als Beispiel wird ein aus ferromagnetischem Eisen und Luft zusammengesetztes 
Gebiet nach Abb. 7 betrachtet. Die Eigenschaften des ferromagnetischen Materials 
werden durch die Magnetisierungskurve Abb. 8 gegeben. Abb. 7 stellt einen Nutenteil 
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des Rotors einer Gleichstrommaschine dar. Die Bestimmung dieses Feldes wurde bisher 

nur mit der Lehmannschen graphischen Methode durchgefiihrt. Im Punkte 28 befindet 
sich eine rechtwinkelige Ecke der Grenze zwischen Luft und Eisen. 

Die urspriinglich gewaihlten Funktionswerte des magnetischen Feldes werden in 

Gaui angegeben. Das Herunterdriicken der Residuen wird im Iterationsverfahren in 

Luft nach GI. 3 und in Hisen 

nach Gl. 8 durchgefiihrt. Die 

GU a ee ¢ qe _0 7 Randwerte werden von der 


2B ¥ Luftseite nach Gl. 16 und 
ee it von der ferromagnetischen 
i or @ Seite nach Gl. 22 errechnet. 
Z | Lal Den Wert im Netzpunkt 28 
B eo eS halt man in dieser ersten 
4 |. zz Phase der Berechnung nicht- 

fie singular. 
Je Man bestimmt weiter 
nach Gl. 39 und 38 den 
13 |p + fiktiven Wert im Netz- 
ta aoa : ze punkt 28: 18000 G, und 
17 4 : 4 + nach einer Interpolation des 
8 a linia > cape Wertes im Netzpunkt 38 
49 rie al werden die Werte von c und 
a 71 | ee d der ersten Wertekonfi- 
te ] } | j guration nach Gl. 28 festge- 

Ub r 

ai | | stellt: c= —7190, d=0. 
a | | Nach der sechsten Werte- 
is St konfiguration konnen schon 
s | Bee | die Werte c und d vernach- 
bi) | Pope | lassigt werden. Dann wird 
| eal eine Korrektur der Induk- 


CN ee ee ee ae tionswerte in der Nahe des 

singularen Punktes 28 durch- 

Abb. 8. Magnetisierungskurve gefiihrt und nach _ dieser 

werden die Induktionswerte 

in diesem Punkte an der Luftseite (16300 G) und an der ferromagnetischen Seite 
(18550 G) ermittelt. 


VIII. Zusammenfassung 


Die in den vorangehenden Abschnitten angegebene Berechnung eines magnetischen 
Feldes mit ferromagnetischem Eisen steht in guter Ubereinstimmung mit der Lésung 
die Lehmann mit Hilfe seiner graphischen Methode erhalten hat. Die Rechenarbeit 
bei der Relaxationsmethode kann leicht mit Hilfe von Nomogrammen und Diagrammen 
oder mit Rechenmaschinen bewaltigt werden. 

(Bingegangen am 30. Mai 1960) 


ig. S. 1, Amm,. 4, 
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Low-Frequency Parametric Amplification 


Using the Variable Reactance of an Induction-Motor Winding 


By Richard Eleo*, Pittsburgh, Pennsylvania (U.S. A.) 
With 15 Figures 


Abstract. Although parametric amplifiers are commonly considered as tools for use at microwave 
frequencies there is no inherent frequency limitation. Parametric amplifiers operating at power 
frequencies utilizing an a.c. induction machine winding as the variable reactance are described. 
The gain and bandwidth relationships for a negative-resistance parametric amplifier using a 
time-varying inductance L (t) are derived and it is shown that this type of amplifier can be used 
at a signal frequency of 60 eps. When an induction machine winding is employed as L (t), rotational 
voltages induced in the rotor of the machine alter the characteristics of the variable reactance. 
This problem is discussed and related to the experimental results which indicate that the amplitude 
of the variable reactance is a function of the rotor speed. This additional effect limits the power 
gain to small values when wp = 2 ws but for ws+ 2 ws the maximum theoretical gains can be realized 
in practice. The advantages of this type of rotating amplifier are the absence of brushes and the 
simple economic construction that is possible. 


Zusammenfassung. Parametrische Verstirker werden vielfach als Bausteine der Héchst- 
frequenztechnik betrachtet; grundsiatzlich gesehen sind sie jedoch nicht an ein bestimmtes Frequenz- 
band gebunden. Der in der vorliegenden Veroffentlichung beschriebene parametrische Verstiarker 
arbeitet im Frequenzbereich der Starkstromtechnik und benutzt die Wicklung eines Wechselstrom- 
Induktionsmotors als veranderlichen Blindwiderstand. Die Beziehungen zwischen Verstarkungs- 
grad und Bandweite werden fiir einen parametrischen Verstiirker mit negativer Widerstands- 
kennlinie abgeleitet, der eine zeitverainderliche Induktivitat LD (¢) hat. Es wird nachgewiesen, 
da diese Art Verstarker mit einer Signalfrequenz von 60 Hz betrieben werden kann. Benutzen 
wir die Wicklung einer Induktionsmaschine als L (t), so verindern im Rotor der Maschine induzierte 
Rotationsspannungen die Kennmerkmale des verdnderlichen Blindwiderstandes. Diese Frage 
wird erdrtert und mit den experimentellen Ergebnissen verglichen. Die letzteren deuten an, dab 
der Scheitelwert des verinderlichen Blindwiderstandes von der Drehzahl des Rotors abhangt. 
Dieser zusaitzliche Effekt begrenzt den Leistungsverstérkungsfaktor auf kleine Werte, wenn 
@p = 2 ws. Der maximale theoretisch mogliche Wert kann jedoch praktisch erreicht werden, wenn 
@p + 2s. Diese Art des rotierenden Verstirkers benétigt keine Biirsten und ermoglicht eine 
einfache und preisgiinstige Bauform. 


List of Symbols 


a real part of impedance, ohms 
b imaginary part of impedance 

c capacitance, micro-farads 

€ transformer voltage, volts, instantaneous 


7 frequency, cycles per second (cps) 

ho pump frequency 

Gp power gain 

I amplitude of ‘difference’ frequency current, amperes, peak 
i amplitude of current 

as amplitude of “signal” frequency current 

a current, amperes, instantaneous 

Pay ts transformer currents in windings 1 and 2 respectively 
L inductance, henries 

ae amplitude of time-varying inductance, henries, peak 
25 average value of time-varying inductance 

L (t) inductance, henries, instantaneous 


* Carnegie Institute of Technology. 
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Ly, Lo. self inductance of winding 1 and 2 respectively 


Ms.» mutual inductance between windings 1 and 2, henries, instantaneous 
M., amplitude of mutual inductance, henries, peak 
Hiss power, watts, average 
Py input power 
Ya output power 
ee power at “‘difference” frequency 
Q quality of tuned circuit 
R resistance, ohms 
Re critical load resistance 
Rr load resistance 
Iie resistive losses, ohms 
yi negative resistance, ohms 
fy period of time function, seconds 
t time, seconds 
Va amplitude of “‘difference’”’ frequency voltage at terminals of non-linear 
element, volts, peak 
Ve amplitude of “‘signal’’ voltage 
Va ‘difference’ frequency voltage at terminals of non-linear element, 
volts, instantaneous 
Un terminal voltage of non-linear element 
Up pump voltage 
OP signal voltage 
Z impedance, ohms 
a phase angle of 7 with respect to v,, radians 
he phase angle of 7 with respect to vz, radians 
B bandwidth factor 
0 mechanical angle of rotation, radians or degrees 
@ phase angle of L(t) with respect to v,, radians 
Wa “difference” angular frequency, radians per sec. 
Om rotor shaft angular velocity, radians per sec. 
Wo resonant angular frequency of tuned circuit 
Wp “pump” angular frequency 
Ws “sional” angular frequency 
Introduction 


Parametric amplifiers are normally considered as tools for use at microwave fre- 
quencies, however there is no inherent limitation on the frequencies at which they can 
be used and operation at power and audio frequencies is possible. The time varying 
inductance needed for such an amplifier can, at low frequencies, be obtained by using 
the windings of an induction machine in such a manner that the terminal inductance 
is a function of time when the machine is driven at some angular speed. One particular 
arrangement which is treated is the single phase short-circuited rotor connection which 
offers the advantage of being an inexpensive, brushless, rotating a. c. amplifier that is 
capable of operating at relatively high power levels (1 kW or more) with substantial 
power gains (20 db). One possible application is as a power amplifier for a 60 cps 
carrier servo system. The ease with which this type of amplifier can be constructed 
at 60 cps makes its use as an educational aid extremely attractive since focus can be 
placed upon the basic principles of parametric amplification without becoming involved 
in the special techniques of making measurements at very high frequencies. 
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I. Principles of Parametric Amplifiers 


The basic principles of parametric amplification were proposed by Lord Rayleigh, 
in connection with non-linear acoustic phenomena!, and by Hartley? in 1930 to de- 
scribe the operation of certain types of modulators. However, in the past ten years this 
type of amplifier has attracted much attention because of favorable low-noise proper- 
ties? and useful gain at microwave frequencies. 

The general character of the phenomenon indicates that the mechanism which provi- 
des gain at microwave frequencies can provide gain at low frequencies. For instance, 
a simple circuit configuration incorporating the parametric mechanism comprises a 
lossless resonant circuit containing a capacitor with plates that can be moved together 
and separated periodically®. Assuming an initial charge Q on the capacitor, movement 
of the plates apart causes the capacitance to decrease and therefore, from Q=VC, the 
voltage rises and hence the energy in the circuit (}@QV) increases as a result of the 
mechanical work on the plates. The resonant circuit allows this energy tc be applied 
to a load under the control of a signal voltage. Thus a parametric amplifier is basi- 
eally a modulator employing a non-linear or time-varying energy-storage element. 
This structure is shown in Fig. 1. The parameter variation is usually produced electri- 
cally rather than mechanically. 
The non-linear element com- 
prises either a capacitor, or an 
inductor, that exhibits a capa- 
citance, or an inductance, which 
is proportional to the terminal 
voltage or current in the device. 
Typical non-linear elements at 
high frequencies are the reverse- 
biased semiconductor diode 
and the saturable inductor. 

The three branches of the 
circuit shown in Fig. 1 are Fig. 1. Schematic diagram of a parametric amplifier using 
termed the ‘signal’, the a general nonlinear element, where wp + ws 
“pump” and the ‘“‘output’’. The 
“signal” circuit comprises an alternating voltage source of a frequency o,/2” and 
a filter, tuned to the signal frequency. Similarly, the ‘“‘pump” employs a source and a 
filter operating at a frequency ,/2 a which is different from the signal frequency. The 
“output” branch contains a filter tuned to either the sum or difference of the “‘signal”’ 
and “pump” frequencies. Terminating the “output” branch is the load R; to which 
power at the “output” frequency is to be delivered. 

When the “signal” and “pump” voltages are applied to the circuit, voltage compo- 
nents at the sum and difference frequencies appear across the non-linear element 
because of the modulation effect of the non-linearity. If the filter in the “output” 


1 Lord Rayleigh: The Theory of Sound. 2nd Ed., Vol. 1, p. 78. New York, N. Y.: Dover 
Publications Inc. 1945. 

2 R.V.L. Hartley: Oscillations in Systems with Non-linear Reactance. Bell Sys. Techn. 
Jour. 15, 424 (1936). 

3 H. Hans and R. Adler: Circuit Theory of Linear Noisy Networks, p. 68. New York, N. Y.: 
The Technology Press of the Massachusetts Institute of Technology and John Wiley and Sons 
a eal Theory of the Ferromagnetic Micro-wave Amplifier. J. of App. Phys. 28, 1225 (1957). 

5 L. Brillouin: Wave Propagation in Periodic Structures. 2nd Ed., p. 186. New York, N. Y.: 
Dover Publications Inc. 1953. 


16 R. Elco: 


branch is tuned to the difference frequency, a current at the frequency (w,— @,)/2 7 
flows in the load. Hence, there is power flow into the non-linear element from the signal 
and the pump sources; and power flow out of the non-linear element to the load. A 
general analysis of this circuit, given by Manley and Rowe® ’, shows that the output 
power can be greater than the signal power and that the circuit can be used as an amplifier. 

From the difference-frequency concept, in the case where w,=2 @,;, the output 
power is at the same frequency as the signal power and the circuit acts as a negative- 
resistance amplifier. This is the class of parametric amplifier considered in this paper. 
An example of this type of amplifier is shown in Fig. 2, in which the non-linear element 
is a saturable reactor and the output branch is combined with the signal branch (since 
the difference frequency is equal to the signal frequency). If the currents in the reactor 
at the signal frequency are small compared with the pump-frequency currents, the pump 
branch and the reactor can be replaced for purposes of analysis by an equivalent 
negative resistance f°. 

This negative resistance is a function of the amplitude (Z,,) of the inductance 


change of the reactor, and the signal frequency, or, 
R, =} Ca Ws). (1) 


To determine the parametric 
variation, Lp, a small signal 
analysis is normally used since 
Coturable dp > ts, and the inductance of the 
Reactor reactor is then a direct function 
of the “pump” voltage. Since the 
“pump” voltage is sinusoidal, the 
inductance can also be expressed 
approximately as a sinusoidal 
function of time: 


L() = L, + 1, sin(o,0+ Oe) 


Signal 


; Th 1 
Fig. 2. Example of a negative resistance parametric 6 torte Ly espe bee sne: ye 


amplifier using a saturable reactor as the nonlinear value of the time-varying induc- 
element where wp = as tance, the amplitude L,,, the de- 


gree of non-linearity and @, the 
phase angle between the pump voltage and the signal voltage. The time-variation of 
energy storage can be produced electrically (as for example in Fig. 2), or mechanically, 
by varying the parameter of the element as a function of time. At power frequencies, 
such as 60 cycles, it is possible to use the windings of a rotating a. c. machine as the 
time-varying, energy-storage element, because of the variation of reactance that can 
be obtained as a function of angular position of the machine rotor. 


II. Analysis of Negative-Resistance Parametric Amplifier with 
Time-Varying Inductance 


The negative-resistance parametric-amplifier equivalent circuit, with a time- 
varying inductance and the output and signal branches combined, is shown in Fig. 3. 


* J.M.Manley and H. E. Rowe: Some General Properties of Non-linear Elements. Proc. 
IRE 44, 904 (1956). 

’ H. E. Rowe: Some General Properties of Non-linear Elements, II, Small Signal Th 
Proc. IRE 46, 850 (1958). i a 

8 $.Bloom and K.K.N.Chang: Theory of Parametric Amplification Using Non-li 
Reactances. RCA Review 18, 578 (1957). ‘ aie 
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The signal and the time-varying inductance are sinusoidal functions of time, and the 
filter comprising the series combination of L and C is resonant at the signal frequency. 
Included in the total inductance of L is the average value of the time-varying induc- 
tance Ly. Any losses that are present in 
the filter and the time varying inductance Load 
are represented by a resistance, R,,, and 
the load presented by the resistance, R; 
(the signal is assumed as having zero 
source impedance). To facilitate the de- 
scription of the circuit, the time-varying Resonant 
inductance will be termed the ‘“pump’’. op 
Because of the presence of the filter, 
the only current that can flow is a current 
at the signal frequency w,/2 a, or one at a 
frequency very close to the signal fre- 
quency. Hence, the current 7(t) can be 
expressed as 


a (t) = TJ, sin (@, t + a). (3) 


: ws Fig. 3. Negative resistance parametric amplifier 
Two important conditions are to be using a time-varying inductance 


considered: 1) when the pump frequency 
is exactly twice the signal frequency; and 2) when the pump frequency is not twice the 
signal frequency. 


Condition (1) w, = 2 a; 
From Kirchhoff’s voltage law, the signal voltage v,, in Fig. 3, is: 


tol iat Li). (4) 


Since the power relationships in the circuit are of primary interest consider the average 
power of each of the terms in (4). 


vy =1(Rr+ #,) + L 


Because the filter losses have been included in the resistance R,,, the terms L = 


and = [i dt vanish when averaged over the period of the signal or an integra] multiple 


of this ‘period. The resulting average power terms are: 
T if % 
t ym [ « (Gad : 
plied =a dt +o li Goi (5) 
0 0 0 
in which 7’ =1/f,. If the power dissipated in the load R, is to exceed the power supplied 
by the signal (a condition which must apply if a power gain is to be obtained) the second 


term of the right side of equation (5) must be less than zero, or: 


5 

1 (reed : 

z)ig lL (4 dt <0. (6) 
0 

Under this condition the pump does not absorb power at the signal frequency but 

rather supplies power. 


The signal voltage and inductance variation in Fig. 3 are given by: 
Op == Og BLO, ES E(t) = L, sin (w,t + ®). (7) 
For the condition a, = 2 w,, the difference frequency (wp — «,)/27 is equal to the signal 
frequency and the expression for the current given by equation (3) is valid. Sub- 
stituting (3) and (7) in equation (5) and integrating gives: 


Ingenieur-Archiv XVI/1 2 
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T 


as cos a =“ (Ry Rn) 5p | Ln Tn? o, sin (2 w,¢ + ®) sin (2@,¢ + 2a) dt + 


0 
x 
cos (2wst+ 2 


ae z| En i? oe Ws : ha 9 [cos (2 Ds t —- ®)| dt. (8) 
0 


The expression : a "cos « is the power supplied by the signal, and J,,? Rz/2 is the 


power dissipated in the load R,. The integrals in (8) will be written respectively as 
P, and P,. Trigonometric substitution in P, yields: 
I 
Pe a (cos ® cos 2% + sin ® sin 2 «) dt. (9) 
0 


Integrated over the limits, (9) is 


2 
Pi = OA” gO = 2.0], (10) 
Similarly, for P,: 
2 
p, = —2“**"' [eos (6 — 2). (11) 
Equation (8) then becomes 
m m me Ss Lim ¢ 
Jn cos a = Bp + Rn — 5 cos (® — 2a) (12) 
where the quantity a cos (® — 2 «) has the units ohms and is an apparent negative 


resistance R, in the circuit. The expression for the negative resistance has the same 
form as the negative conductance, cos ® given for a time-varying capacitor 
parametric amplifier by Bridges?. 

The power gain of the circuit, when w, = 2 @, is: 


Pout Ln? Rr 


Gp ct Pin Deli Vir cosnce (13) 
and substitution of (12) into (13) gives 
1 
Gp “ \ Em @s Lin . © 2 : (14) 
ee 2 Rr COS ( a) 


Thus the gain is dependent on the signal frequency «,/2 x, the amplitude of the time- 


varying inductance L,,, and the load R;. If the term yor 
1+ &,,/Rz the gain is infinite and the circuit becomes unstable and may oscillate if 
there is a suitable disturbance (such as noise) tc start signal currents flowing. 

The phase angle, «, in (14) is determined by the phase angle of the total circuit im- 
pedance. Since the pump is a negative resistance and the filter is resonant at «,/2 2, 
a equals zero for the condition that w, = 2 w;. Then to maximize the gain with a con- 
stant load f;, the phase angle ® of the inductance variation with respect to the signal 
frequency should be zero. 


cos (® — 2 «) is equal to 


Condition (2) wy, + 2 a; 


In this case, the difference frequency, wz, is no longer equal to the signal frequency 
and if the difference frequency is much greater than, or much less than the signal 
frequency, no current at the difference frequency can flow because of the filter. The 
analysis presented in Appendix I shows that the integrals P, and P, are zero and the 


* T. J. Bridges: A Parametric Electron Beam Amplifier. Proc. IRE 46, 494 (1958). 
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pump no longer represents a negative resistence. The gain therefore becomes less than 
unity: 
i 


fe fat ares (15) 


However, if the difference frequency is very close to the signal frequency, current 
at the difference frequency can flow because of the finite bandwidth of the filter. This 
current produces a secondary modulation effect which causes the amplitude of the 
current given by (3) to vary slowly with time. If the frequency at which the amplitude 
varies is small compared to the signal frequency, the method of equivalent linearization 
makes (3) a valid approximation. However, if equation (14) is to be used to describe 
the power gain in this case, the effective negative resistance must be corrected by some 
factor which includes the impedance of the circuit as seen by the pump at the difference 
frequency. 

The impedance seen by the difference-frequency voltage that appears across the 
pump is the impedance of the load and the filter. Because of the filter, the difference- 
frequency current, J,, is reduced in amplitude; and a, its phase with respect to the 
difference-frequency voltage across the pump, is increased. The approximate power 
flow out of the pump at the difference frequency is: 


Py = Iq Va COS da, (16) 
which decreases as the absolute value of (ws — w;) increases. Since the magnitude of 
the negative resistance is proportional to P,_, the negative resistance given in (14) is 
multiplied by the factor £, 

8 iy ji : ie — ws) L — 1/(@p — @,) aie |e — =" (17) 


Ws 


Ri + Rm 
which is derived from (16) in Appen- 
dix IT. 

The average powers in (12) were 
obtained by integrating the instan- 
taneous power over one cycle of the 
signal, treating ®, the phase angle 
of the pump, as a constant. How- 
ever, when the pump frequency is 
not exactly equal to twice the signal 
frequency, @ changes from the be- 
ginning of one cycle of the signal to 
the beginning of the next cycle. If 
‘the pump frequency is greater than 
the signal frequency, ® leads the 
signal and increases with time at a 
rate which is proportional to the 
incremental frequency, (wy — ;)/2 7. 
Hence the phase angle of the pump 09 0.95 io 08 nF 
is a function of time given by 


Normalized Gain (IGmoxl/IGI) 


Normalized Pump Frequency wpe/2w, 


D = (wp — 2 ws) t. (18) Fig. 4. Calculated response of filter (A) and gain 

bs cad of negative resistance parametric amplifier (B) as a 

Substitution of (17), (18) and function of the pump frequency where: 1/ws Rt ¢= 
(III 4) (Appendix ITT) into (14) gives =10.0 and ws Im/2 Rr = 0.8 


10 N. Kryloff and N. Bogoliuboff: Introduction to Non-linear Mechanics. Princeton, N. J.: 
Princeton University Press. 1943. 
2° 
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the general expression for the gain of a negative-resistance parametric amplifier 
with a timevarying inductance as the pump: 


GS : ; (19) 


a 
tans | eee 
Rr + Rm 
The normalized peak gain (since G, is now a function of time), and the normalized 
impedance of the circuit, |Z,,|/|Z| excluding the pump, for a filter of Q=10 and 
wo; Lm|2 Ry of 0.8, are shown in Fig. 4 as functions of frequency. The bandwidth 
between the points of half power gain is seen to be approximately 0.02 w,/2 @,. Hence 
for a signal frequency of 60 c/s the expected bandwidth is 0.6 cps. 

The decrease in the effective bandwidth from that of the filter alone is caused by 
the negative resistence R, that appears in the circuit, since this cancels some of the cir- 
cuit loss and raises the effective Q factor. 

The asymmetry of the gain curve is a result of the frequency dependence of the nega- 
tive resistance as given by the bandwidth factor # in (17). 


Rin pe Ws Lim 
Rr 2 Rx 


1+ B cos | (mp — 2@s) t — 2 


III. Use of A.C. Machines as the Time-varying Inductance 
Induction Machines 


The feasibility of constructing a negative-resistance parametric amplifier to operate 
at signal frequencies of 60 cps is largely determined by realizable values of the ratio 
of w; L,, to R, which control the gain in (14). If 6=2.« the gain becomes infinitive 
when 

ws Lm vm 


= > 20 
Stator Axis 2 Rr as Rt (20) 


Equation (20) indicates that with a time-varying 
~~ @ Rotor Axis inductance of amplitude 5-10-2 henrys and 
negligible losses (R,, ~ 0) used in the circuit of 
ua Fig. 3, reasonable gains should be obtained with 
/ a load resistance, R;, of 10 ohms at a signal fre- 
if Se quency of 60 cps. Another important point to be 
seen from (20) is that the ratio of w, LZ, to the 
Zt losses #,, should be large. For instance if the in- 
Lag ductance given above actually has an internal 
3 4 Ro resistance cf 5 ohms, the load resistance for maxi- 

a ee / | mum gain is reduced to 5 ohms. 
as An induction machine can be described as a 
transformer in which the primary winding (stator) 
is stationary and the secondary winding or wind- 
piianer Ghantcelnchited irolent auaaeyon ings (rotor) rotate about an axis perpendicular 
machine where it Sail wenn (Thoen the axis of both windings’. From Fig. 5 there is 
magnetizing currents are not included @ mutual inductance M,, between the stator and 
in this representation) rotor windings, which is a function of the angular 
position of the rotor. If the rotor is short circuited 
and rotated slowly, the inductance seen at the stator terminals as a function of O can 

be found from the solution of the equations 


1 A. E. Fitzgerald and C. Kingsley, Jr.: Electric Machinery, Chap. 2. New York, N. Y.: 
McGraw-Hill Book Co. 1952. 


e| 


of Rotation / | 


Fig. 5. Transformer equivalent circuit 
for a single-phase stator and single- 
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di, , ad ‘ 
Cae Ly Hy + ay (Mie te) 
: dixand ; 
= Us Jaa +- ee + ay (M2 11), (21) 


where L,, is the self-inductance of the stator and L,, is the self-inductance of the rotor. 
Solution of (21) for the ratio of e,/i, when R, ~ 0, gives for the effective inductance 
of the stator 


Silene 
L(@).= 51. (1 + cos 20). (22) 


For confirmation of this expression, measurements were made of stator inductance 
for various angular positions of the rotor (at rest) with a single-phase short-circuited 
rotor in the M. I. T.-Westinghouse Generalized Machine! recently available in the 
laboratory. The parameters 
of the machine with this 
connection were: M,, = 0.15 
pose hy; L,, = 0.4hy;.L5.= 
0.05 hy and the static mea- 
sured inductance (Fig. 6) 
agreed with that predicted 
by (22). 

The amplitude of the static 
inductance variation obtained 
in this fashion is of the right 


order of magnitude for practi- 0 60 120 i830 240-300 
cal use in a parametric am- Angular Position of Rotor @ (degrees) 


plifier. However, ei diffi- Fig. 6. Measured inductance as a function of angular position, 
culties arise because of addi- O, for the ‘Generalized Machine”’ with a single-phase short- 
tional induced voltages and circuited rotor 


currents in the rotor when 

the rotor is driven at some angular speed. These rotational currents will affect, to some 
degree, the nature and order of magnitude of the variation of inductance seen at the 
stator terminals. For an angular speed of , the mutual inductance becomes 
M,.cos w»,t and if the reactance (w, L.5) is assumed to be small compared to the 
rotor resistance, from (21) the voltage across the stator is 


2 Oe OC: © 
NM. “poo 


Inductance (henries) 


0.! 


di, 
dt 


eg fea 


dM,\2] , di [3dM,, 
es = — fh ae Maa t+ (Sat) | + ae ae 


a2) | + Gt | Got Ma — Ba Das] + 


[Mre}*- (23) 
The explicit form of the time variation of the inductance is not clearly evident 


from (23) which indicates that the machine can not be exactly replaced by an equi- 


: é : eee d . a 
valent time-varying inductance since it is not of the general form e = Lo+t a 


which gives the voltage across a time-varying inductance. However, physical reasoning 
(and the experimental results of IV) suggest that at any speed other than w,, equal to 
@;, which will be termed the synchronous condition, the machine can be approximated 
by a time-varying inductance of the form given in (22) in which the angle @ is a function 
of time (OQ = w, t). When the rotor is driven at a speed less than synchronous speed 
the rotating fields of the rotor cause changes in the stator fields in much the same 
manner as in the static case, when the rotor is rotated slowly. 


12 D.C. White and H. H. Woodson: A New Electromechanical Energy Conversion Labora- 
tory. AIEE Paper 57-603. 
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At synchronous speed the machine is similar to a synchronous machine which is 
operating with a short-circuited field winding and with no direct-current excitation 
of the field. Hence there can be no transfer of mechanical or electrical power in the 
machine. Since a time-varying inductance implies a transfer of energy stored in the field 
to an external circuit, the machine when driven at synchronous speed is not a time- 
varying inductance and appears only as a positive resistance (losses in rotor and 
stator) which dissipates power. Further justification of this reasoning follows from the 
torque-speed characteristic of a single-phase induction machine, which shows that at 
speeds less than or greater than synchronous speed power transfer does take place. 
The phenomenon of self-excitation of induction machines (induction generator)’, 
also implies such power transfer considerations. 

From this reasoning and the results of Appendix IV which compares the difference 
frequency voltage, Vz, for a time-varying inductance to that for an induction machine, 
the L (t) seen at the stator terminals of a 2 pole single-phase induction machine with a 
single-phase short-circuited rotor will be taken as: 

L (t) ~| #2 M,,*| sin (at + ®) (24) 
where w, = 2 @m. The amplitude of (24) was determined for the condition L,,. ~ 0 
and for any practical machine this condition is not very realistic. However, (24) is 
useful qualitatively because it shows that the variation of the amplitude of L (t) with 
@» is linear and for wy, close to 2 w, it can be used as an approximation to L (?). 


Consideration of Other Possible Induction-Machine Connections 


a) Single-phase stator and polyphase short-circuited rotor. This connection provides 
an L(t) similar to that of the single-phase rotor connection except at speeds near zero 
at which L(t) becomes a constant. 

b) Single-phase stator and single-phase rotor; series connected. If the rotor and the 
stator of the machine are connected in series, the rotor and stator currents of equation (21) 
are equal (when e, + 0) and the L (t) seen at the terminals of this series connection is 


L(t) = (Log + Ly) + 2 My, COS Wm b, (25) 


which is definitely a time-varying inductance, as long as R, ~ 0. However, to obtain 
the proper frequency for L (¢), when w,/2 2 is equal to 60 cps, the machine must be 
driven at twice normal speed, which for most machines is inadvisable. 


Salient-Pole Synchronous Machines 


In a salient-pole a.c. machine* the reluctance of the stator magnetic circuit is a 
function of the angular position of the rotor and if the rotor is driven at an angular 
speed of w,, the inductance of the stator is, to a first approximation 


EL (t)== Ly +L,8in.%, 03,, 4. (26) 


where m equals the number of poles. When used as a time-varying inductance the 
machine is operated without field excitation. 

The practical use of such a machine for our purpose is, however, limited by the 
relatively small amplitude of the time-varying inductance of the stator (5- 10-4 hy) 
for a selected salient-pole machine. The resistance of the machine winding was 0.25 ohms 


** H.M.McConnell: A Self-Excited Induction Generator with Regulated Voltage. ATEE 
Paper 54-314. 


* The machine is assumed to have no amortisseur windings. 
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which, taken into account with other power losses, made it impractical to obtain gain 


23 


in such a machine system. However, it may be possible to design the machine such that 
the inductance variation is large enough to provide practical power gains. 


IV. Experimental Investigation of Parametric Machine Amplifiers 


Two different machines were used as the time-varying inductance in the experimental 
mvestigation. The first was the M. I. T.-Westinghouse “generalized machine’’* 
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Fig. 7. Experimental circuit connections for Fig.8. Experimental circuit connections for ma- 
machine “A” where: [y= 0.22 hy; Im=0.2hy; chine “B” where: Ly, = 0.056 hy; Im= 0.045 hy; 
L’=0.04hy; C= 28.0uf; Rm =3.00; Re-=34.00 L’= 0.056hy; C= 61.6uf; Rn = 5.0Q; Re = 3.50 


(Machine “A’’) with the single-phase rotor connection, and the second was a commercial 
3 phase wound rotor induction motor (Machine “‘B’’). The circuit connections and para- 


Mechanical Power P., (watts) 


0.9 0.95 10 1.05 Ll 
Normalized Pump Frequency “P/2w, 


Fig. 9. Mechanical power supplied by the drive 
motor for machine “‘A”’ where Rr = 28 Q 


Fe Seen 12) pera. 


meter values for each machine are given 
in Figs. 7 and 8 in which the average 
value of the measured JL (t) is included 
with the inductance of the filter (LZ 
and C). 

From the theoretical expression (19) 
the critical load Rc needed for the power 
gain to become infinite, was calculated 
for the circuit values shown in Figs. 7 
and 8. Then to check these predictions 
and other theoretical conclusions the 


i(t) 


r 1/60 sec. 


Fig. 10. Load current for machine “B” at a 
pump frequency of maximum gain, 130 eps, 
when Rr> Re 
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0 4 8 12 16 ; : 
Load Resistance R, (ohms) Signal Power in P. (watts) 
Fig. 11. Theoretical and experimental power Fig. 12. Theoretical and experimental output 
gains for machine “‘B’’ as a function of the powers as a function of the input power for 
load resistance, Rz, where fp = 130 cps machine ‘“‘B”’ where fp = 130 cps 


Power Gain (R/P,) 


0.95 1.0 105 1.10 y 115 Normalized Pump Frequency wp/2w, 
Normalized Pump Frequency “P 
pe, y 2u, Fig. 14. Power gains of experimental para- 


Fig. 13. Output power as a function of the metric amplifier as a function of the normalized 
normalized pump frequency when Ri< Rk. pump frequency where ws is constant and cp 
and v;=0.0: machine “A”, Rz = 210; equals 2 wm [] machine “‘A’’, Rr =28 0 


machine “B”, Rr=3.6Q A machine “B’’, Rr = 4.50, 
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following test cycle was conducted for each machine. Initially the input signal power 
was held constant, with a particular R,, and the power-gain determined as a function 
of the (angular) mechanical speed w,, of the driven induction machine. Then the 
mechanical speed was held constant at the value which gave maximum gain with the 
particular R;, and the gain determined as a function of the signal input-power. Finally, 
Rr was reduced below Re so that the circuit became oscillatory and, in this condition, 
the output power was observed as a function of the angular speed of the induction 
machine. The signal frequency for these tests was 60 cps. 

At a pump frequency (w= 2 ,) of 120 eps, the point at which maximum gain 
should normally (when L,, is constant) occur, the gain for both machines was a mini- 
mum, which is expected since L,, is zero at 
®p = 2 @;. The mechanical power supplied to 
the shaft of L(t) by the drive motor was also 
a minimum at this speed (Fig. 9), which indi- 
cates that the dip in gain was not a consequence 
of extra losses arising at this frequency. 

The pump frequency at which maximum gain 
was obtained was 135 cps for machine “‘A”’ and 
130 cps for machine “B” in the single-phase 
short-circuited rotor connection. Theload current 
at this pump frequency for machine “‘B”’ is shown 
in Fig. 10, in which the frequency of the modu- 
lation envelope is 10 cps. For given values of R, 
the maximum gain obtained corresponded to 
the theoretical values expected from (19). The 
load resistance for oscillatory conditions checked 
quite closely; for example, Rc calculated for (b) 9 10 m 
machine ““B” was 3.54 ohms and the experi- 
mental Re was 3.6 ohms. The theoretical and 


IL ml (henries) 


Relative Gain 


experimental gains of the machine ““B” amplifier 

as a function of the load resistance FR; are given 4 oy yan 

in Fig. 11. Fig. 12 shows the theoretical and ex- 2 AG 

perimental curves of the output power as afunc- 2 i st 

tion of the input power for machine “B” and ©& 

indicates that the amplifier is reasonably linear (c) "9 LO rm 

over the operating range shown. Normalized Pump Frequency wp/2w, 


For values of R; less than Rc the circuit 
always behaved as an oscillator. This mode ot ee 
operation was initiated by adjusting the pump ,, 2 constant are shown in (a) andcby. 
to the frequency at which maximum gain was Combining (a) and (b) gives the gain curve 
obtained and then applying a signal voltage of in (ec) where Im is a function of wp 
one or two volts. The load current then in- 
creased as oscillation began, and once this occurred the signal voltage could be reduced 
to zero and the signal source short-circuited without affecting the current in the load. 
The output power as a function of pump frequency when the circuit is an oscillator is 
given in Fig. 13. 

The modulation of the load current shown in Fig. 10 is caused by a secondary modu- 
lation effect between the signal frequency and the “‘difference” frequency. This is 
explained by the time-varying phase angle ®, given by (18), when the pump frequency 
is not equal to twice the signal frequency. 

The experimental pump frequency at which maximum gain is obtained is displaced 
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from «,/2 a, = 1 (Fig. 14) because of several factors. The first of these is the change 
in the amplitude (L,,) of L (¢) in the region where ,/2 , = 1. The experimental values 
of L» as a function of pump frequency are shown in Fig. 15a. [These values were ob- 
tained by solving (19) for ZL, and substituting the experimental values for Gysiabtgs 
R,.] The variation of L,, with w, is approximately linear as was expected from (24). 
Since R, is proportional to L,,, the gain curve of Fig. 15b will be modified, as shown in 
Fig. 15c, in which the notch is caused by the decrease in Ly at w,/2 wo, = 1. The diffe- 
rence in the height of the peaks AG is due to the asymmetry of the gain and might 
also be caused, in part, by a change in the average value of L(t) which would cause a 
change in the resonant frequency of the filter. 


V. Conclusions 


1. The experimental results obtained with an induction machine as the time- 
varying inductance indicate that the theoretical treatment of Parts II and III of the 
paper does in general describe the physical circuit and that predicted gains can be ob- 
tained at signal frequencies of 60 cps. The derivation and reasoning on the equivalent 
circuit for the induction machine hinges on the assumption that the machine can be 
treated as a linear system. This description is sufficient when relatively small powers 
are dissipated in the load. However, when the machine is operating at high power 
levels it is obvious that the non-linearity due to saturation of the iron will cause the be- 
havior of the amplifier to deviate from that which is predicted by the linear analysis. 
This is evident when the machine is operated as an oscillator. The angular speed at 
which maximum power output is obtained is shifted slightly from the speed at which 
maximum gain is obtained when the machine is operated as an amplifier (Figs. 13 and 
14). The small differences in speeds and the similarity of the curves obtained in each 
condition indicate that the linear analysis is accurate over the range of power levels 
investigated experimentally. 

The amplitude (L,,) of the time-varying inductance that the machine presents 
at its terminals is important only over the range of angular speeds for which gain 
would be obtained if L,, were constant, since beyond this range the resonant circuit 
prevents the necessary currents from flowing and the gain is less than unity regardless 
of the value of L,,. Hence the description of L,, as a function of the pump and signal 
frequencies as presented is sufficient. 

2. This general demonstration of 60 cps parametric amplification presents a new 
and unusual application of an induction machine. The advantages of such an amplifier 
are that it is extremely simple and it does not require brushes or slip rings. The only 
disadvantage is the dependence of the power gain on the kind of load that is presented 
to the machine. However for fixed resistive loads optimum gain can be obtained by 
using a transformer to adjust the effective load presented to the machine. 


Appendix I 
Evaluation of the integrals P, and P, when w, + 2 ws 


When @, is not equal to 2 w,, the integrals P, and P, from equation (8) become, 


IG 
a @s Lm I?m 


Py = 577" | sin (n @,t + ®) sin 2 w, t dt (I 1) 


0 


and 


Ti 
@s Lim 1% m 


Bs ae sin? w, t cos (n w, t + ®) dt ({ 2) 


0 
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in which 7 is the period of the lowest frequency component contained in the integral 


and the pump frequency is given by » w,, where n is any integer that is not equal to 
two. 


Then making trigonometric substitutions in on 


4 . 
@s Lm Ln? 
P= 79 | {leos (n w,t — 2 @,t) — cos (n ast + 2 w, t)] cosB + 


0 
+ [sin (2 ot + nw; t) + sin (2 wt — n a; t) sin ®]} dt. (i3) 


Since 7’ is equal to the period of the lowest frequency component inside the integral 
the value of the integral is zero. Similarly for P,, 
P —_ 2 Os Lm In? S 1 
- jel {cos n w,tcos® — t [cos (2 wt — nw; t) + cos (2 a,¢t + nw, t)] cos® — 
0 


—sinn w,tsin® + $ [sin (nw, ¢ + 2 @,t) + sin (n wo, t — 2, ¢t)] sin®} dt, (I 4) 


which also becomes zero. Hence for the conditions given P, and P, are zero. 


Appendix II 


Determination of the negative resistance R, for wy + 2 as 


When @, is not equal to 2 w, the voltage across the time-varying inductance V, 
will have components at the signal frequency w,/22, the difference frequency 
(@, — o;)/2 2 and the sum frequency (w, + @;)/2 2. Since currents at (w, + @;)/2z 
cannot flow because of the filter only the wm, component and (w, — w,) component 
need to be considered. If the difference frequency is close to the signal frequency there 
will be two components of the current flowing in the circuit, one at w,/2 a and the other 
at (w, — o;)/2 x. Then the current 7’ (t) is 


i’ (@) = I, sin w, t + Iysin (@» — o,)t. (II 1) 


The inductance variation is 


TAG) oe iene tage (II 2) 
and so substitution into the term, se! {ZL (é) a’ (t)} 


of equation (4) gives 
Vin = Lm {I, [os sin wp t cos wt + wy SiN @; t COS wy, t] + 
+ Iq [(w@p — ws) Sin ws t COs (@p — ws) t + wp COS w, t sin (wy, — os) t]} (II 3) 


and the difference frequency component of the voltage is obtained by trigonometric 


substitutions as 


Ui xe a an (wp — ws)t. (II 4) 


Thus the rms value of V, is proportional to the product of the signal frequency 
component of current and the difference frequency reactance (wp, — s) Ln. 

The negative resistance R, appearing in the circuit of Fig. 3 is proportional to 
the power P,— supplied by the time-varying inductance at the difference frequency 
(@», — @s)/2 7% or 

Py = La Va COS qd ~ R, (II 5) 
where J, and V, are rms values and the phase angle «, due to the impedance Z which 
the circuit presents to Va, is 


wd L — lod 4 (IT 6) 


— -1 
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When the difference frequency w,/2 a is equal to w,/2 7 the negative resistance R, 
is given by equation (12) but when «,/2 x is not equal to the signal frequency the nega- 
tive resistance R,’ is given by 


Heyy iced Wigs B. (II 7) 
The ratio of the component powers at the condition wz equal to w, and at , not 
equal to a, is: 
TeV ae ee in 
Ta’ Va’ cos aa Rn B 
From equation (II 4) it is seen that for w, near w, the current J, can be treated as 
a constant. Then J,’ is in terms of J, 


(II 8) 


La = A (@p x Ws) G (II 9) 
and V,' in terms of Vz is a 
Vii Lees Lone, (II 10) 
Substitution of (II 9) and (II 10) into equation (II) 8 gives 
we (Z| 1 
(ap ae eowas (IT 11) 
or 
B = [Fi cos ag | P— ; (II 12) 
where 
|Z| = (#, + Rr) and 
|Z’ — VRn +- Rr)? + (wma L — 1/aac)?. 
From equation (II 6) 
cos G7 = ia IL 
disee Pa ( 13) 
and substituting in (II 12) gives 
(Ri + Rm)? (wp — ws)? 
Simplify; B= (Rpt Rn? + (oa L —VYooue)*] of” a 
Simplifying 
wo wd L — 1/mac)? |-1 [wu]? 
pj ee ee (TT 15) 
Mas ®q = Wp) — w, and @, is constant. When w, is constant and «, is variable f 
ecomes 
wa L — 1/wae}2 |-1 [22 
pol eae ey ee (Ul 16) 


Ap pendix ot 
Determination of the phase angle « given in equation (5) 


Since the current 7(¢t) and the signal V, are assumed to be sinusoidal the circuit 
in Fig. 3 can be represented by an impedance Z which has the form: 


Z=-a+ 7b (IIT 1) 
Then, from Fig. 3, a and 6 are: 


a= (Bees ao Ry rr R,) 
b = (a, L — 1/a,c) for: Ry = (Be 2R,) (IIT 2) 


in which £, is the negative resistance given in (12). 
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The phase angle, «, of the current ¢ (t) is given by 
a= tans. (IIT 3) 


When the filter is resonant at w, and w, = w», b = 0 and « is also equal to zero. From 
the physical reasoning of Appendix II the negative resistance is a function of the diffe- 
rence frequency, and if w, is different from the resonant frequency of the filter w, and 
@p = 2 w,; the negative resistance will be decreased by the factor f. 

As a first order approximation for « the negative resistance can be assumed to be 
very small when @, is not equal to w», hence « becomes 


os DT — l/wse 


a = tan=! eae (III 4) 


Then for small changes in the pump or signal frequencies about the frequencies 
@p = 2 w, and w, = w, equation (III 4) may be applied. 


Appendix IV 


Determination of L (t) as a function of w, and wy 


For the condition when the difference frequency is very close to the signal frequency 
the current (7,) in (23) can be taken as 
(HU, Sin; (IV 1) 
and the mutual coupling M,, is 
Ma MN. COS Wy t (IV 2) 
where: wp = 2 @m. 
Substitution of (IV 1 and 2) into (23) gives 


Im M 2 Wm" os 3 . 2 Or . 
ns z {| - “gait g Os On| sin (2 @m + os) ¢ + —5- sin ot — 
m be 3 ‘ 
[2s +S — 2 @ oq] sin (2 om — 0.) t} (LV 3) 


whence the difference frequency voltage is: 


Qs. M,,2 Im 
— (Om a8 Ws) On — 3) az) R, 


sin (2 Wn — ws) ?. (IV 4) 


For a pure time-varying inductance the difference frequency voltage from 


e=-Lo oi (where 1 = J,, sin w,¢ and L = L,, sin 2 w,, t) is 


i (o — 3) Lim Im S10 (2 Om — @;) t. (IV 5) 


Then equating (IV 4) and (IV 5) % 
(@p — 2 ws) M,,? 
2 R, 
which indicates that Z,, varies linearly with w, and at w, = 2 as, Lm is zero. Since the 
gain relationships (19) holds for positive or negative values of L,, (because of the 
time variation of ®), the amplitude of the inductance is: 

(wp —2 Ws) Mi, 
Lg = oR, 
This expression (IV 7) is only approximate and the error increases as «» is changed 
from w, = 2 ,, but in the region where w» is close to 2 w, (IV 7) does provide a suffi- 
ciently accurate description of Ln. (Received September 12, 1960) 


Tip a 


(IV 6) 


(IV 7) 
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Wirmespannungen in Dreiecksfliigeln* 
Von L. Béswirth, Solothurn 
Mit 18 Textabbildungen 


Zusammenfassung. Die Warmespannungen in einem aus vier flachen Kegelschalen be- 
stehenden Deltafliigel werden berechnet. Das Temperaturfeld wird als symmetrisch beziiglich 
der vertikalen Symmetrieebene des Fligels vorausgesetzt, kann aber sonst beliebig sein. Ein 
Zahlenbeispiel wird gerechnet. 


I. Problemstellung 


Die Steigerung der Fluggeschwindigkeit bei modernen Flugzeugen bringt eine 
immer gréBerere kinetische Aufheizung der Fliigel mit sich. 

Gegenstand der vorliegenden Arbeit ist es, die in Zusammenhang damit auftreten- 
den Warmespannungen unter den Voraussetzungen der linearisierten Elastizitats- 
theorie zu ermitteln. Die Arbeit beschrankt sich auf deltaformige Fliigel. Diese werden 
durch ein Modell, das aus vier geeignet verbundenen flachen Kegelschalen konstanter 
Wandstarke besteht und im Innern keine Versteifungsrippen besitzt, angenahert. Das 
Temperaturfeld in der Schalenmittelflache wird als symmetrisch beziiglich der verti- 
kalen Symmetrieebene des Fliigels angenommen, kann sonst aber beliebig sein. Die 
Temperaturverteilung iiber die Schalendicke wird durch eine Gerade angenahert. 


II. Grundgleichungen fiir flache Kegelschalen 


Die Grundgleichungen nach Marguerre in der Schreibweise von Wlassow ? 
mit Beriicksichtigung eines tiber die Schalendicke linearen Temperaturfeldes lauten fiir 
flache Kegelschalen in einem Polarkoordinatensystem 7, y mit dem Ursprung in der 
Kegelspitze 

1 aw 


(Eh)A yk =e = =aVe P 
(1) 
yy2w+N-06+ SF =i ap) Ver. 
Dabei bedeuten 
EL... Elastizitatsmodul 
h Schalendicke 
0 fF ... Airysche Spannungsfunktion 
ty . Parameter fiir die Flachheit der Schale, 
Abb. 1 
w ... Verschiebung senkrecht zur Schalen- 
mittelflache (in Richtung z, Abb. 2) 
c-rk /e=% qa ... Warmedehnungskoeffizient 
7 krmmung T ... Temperatur in der Schalenmittelflache 
N ... Biegesteifigkeit 
uw ... Querdehnzahl 
I’... Gradient der Temperatur in Richtung der 


Abb. 1 Schalennormalen (in Richtung z) 


z The research reported in this paper has been sponsored by the United States Government 
(Principal investigator: Prof. Dr. H. Parkus). 


* K.Marguerre: Zur Theorie der gekriimmten Platte groRer Forminderung. P i] 
. 5 f - Proc. 5th Int. 
Congr. of Appl. Mech., 8. 93—101. Cambridge, Mass. 1938. : ey 


W.S.Wlassow: Allgemeine Schalentheorie und ihre Anwendung in der Technik. § 
‘ Z * » 9. 239— 255. 
Berlin: Akademie-Verlag. 1958. (Ubersetzt von A. Kromm.) ; : eee 
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Aus den Lésungen der Gl. (1) F (r, 7), w (r, gy) kann man die interessierenden Schnitt- 


groBen ermitteln?: 


Schnittkrafte 
“3 ler 1 oF 
m= | 0 de ~- * eg? r or 
agree h ghaylt 3 
my = | Op Oe = “ap (2) 
on 1 oF 1 #F 
tre = | Freq I = 5 aps r:~<Ordy 
Schnittmomente 
<a aw 1 dw 1 ow 
m, = —{o,2dz = N| 5S tu (= ag? - Fr) a(l+y)T| 
here 1 ew 1 dw aw 
Mg = — | ogzdz = NE, ag? Phareiitinais le ee) r| (3) 
eal Eh 1 @w 1 dw 
Mo = | Org? de = 12 (1+ wp) ( r orep 7 a 


NV 


? = Const 


Abb. 2 
Querkrafte 
= |oede=NS[-Pwta (lty)l] 
ea (4) 
de =| ey d2= N= [—Vw+e (+) 7] 


Das Integrationsintervall geht jeweils von — bis + a . Diese SchnittgroBen sind in 


Abb. 2 an einem Element der Schalenmittelflache in positivem Sinn eingetragen. 
z ist senkrecht zur Schalenmittelflaiche und so gerichtet, daB 7, gy, z ein Rechts- 
koordinatensystem bilden. 


III. Die Randbedingungen 
Vereinfachtes Modell eines Deltafligels 


Im folgenden beschranken wir uns auf ein Gebilde, das sich aus vier flachen Kegel- 
schalen in der in Abb. 3 skizzierten Anordnung zusammensetzt. Es soll naherungsweise 
die Verhaltnisse bei einem Deltafliigel wiedergeben. Die Verbindung der einzelnen 


2 Beztiglich einer zusammenfassenden Darstellung der Theorie der Warmespannungen sei auf die 


zwei Biicher 
E. Melan und H. Parkus: Warmespannungen infolge stationarer Temperaturfelder. Wien, 
Springer-Verlag 1953, 
H. Parkus: Instationdre Warmespannungen. Wien, Springer-Verlag 1959, 

verwiesen. 
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flachen Schalen untereinander wird als biegungssteif angenommen. Ferner soll das Ge- 
bilde symmetrisch in bezug auf die x-y-, und z-y-Ebene sein. Infolge der Flachheit der 
Schalen unterscheidet sich ein Koordinatennetz in der Schalenmittelflache nicht von 
seiner Projektion auf die Grundflache. 

Fiir das Temperaturfeld machen wir die Annahme, da8 es beziiglich der y-z-Ebene 
symmetrisch ist. Die GréBen 7’ und I’ des Temperaturfeldes werden im allgemeinen auf 
der oberen und unteren Schale verschieden sein. Wir spalten sie in einen beziiglich der 
v-y-Ebene symmetrischen und einen schiefsymmetrischen Anteil auf und betrachten 
die Randbedingungen getrennt fiir jeden Fall. 


Symmetrisches Temperaturfeld 


Die Spannungen und Verformungen sind in jeder der vier Kegelschalen in entspre- 
chenden Punkten gleich. Es geniigt daher, wenn wir eine Schale allein untersuchen. 
Wir verwenden jetzt ein Koordinatensysten x, y, z (Abb. 3). Die Schale wollen wir 
nicht durch die Kriimmungen charakterisieren, sondern durch eine Gleichung 


= We (x, y); 
Y NPALA 
- Tinn 
DB 
INO ae ee 
a7 23 ZW 
Y 
hand I 
Son UN , Lins - 
Fal Os 
or Abb. 5 
—> 
ZL 
e@ kandi 
Zi 
pee ate 
Abb. 3 


wobei W den z-Abstand der Punkte der Schalenmittelflache bedeutet. Fiir ein symme- 
trisches Temperaturfeld lauten die Randbedingungen: 


Rand I, II: 


aw 
Ww = 0, a= 0, Nns = 0; Nan = 0, (5) 


ap = 
wobei » und s Normale und Bogenlinge lings des Randes bedeuten (Abb. 3). 


Rand III: 


ow Aw 
PS Ala OUae ten uaa. 


=0, tm=0, u=0. (6) 


u bedeutet die Verschiebung in w-Richtung. Abb. 4 verdeutlicht den Zusammenhang 
zwischen Querkraft und Normalkraft auf dem Rand. 
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Schiefsymmetrisches Temperaturfeld 


i Die Spannungen und Verformungen sind hier ebenso wie das Temperaturfeld be- 
_ziglich der a-y-Ebene schiefsymmetrisch. Die Randbedingungen lauten: 


Rand I, I: 
M=90, dat ees eo fen, 6a, B= 0, (7) 


v bedeutet die Verschiebung in Richtung y. Den Zusammenhang der Ersatzquerkraft 
mit der Schnittkraft lings der Rander I und II zeigt Abb. 5. 


Rand IIT: 


fe SF gyal >, = 9, Mns=90, w=0. (8) 


W ist fiir flache Schalen klein im Verhaltnis zu den tibrigen Abmessungen des Fliigels; 
_ebenso ist 


ow 
too 


Ks ergibt sich sowohl beim symmetrischen wie beim schiefsymmetrischen Tempera- 
_turfeld ein gemischtes Randwertproblem. 

In vielen Fallen wird die Annahme einer starren Binderscheibe (Abb. 6) in der y-z- 
Ebene der Wirklichkeit niher kommen als die Annahme, da8B die Schalen dort biege- 
_ steif verbunden sind. 

Wir wollen daher noch kurz auf die Randbedingung fiir diesen Fall eingehen. 


Starre Binderscheibe in der y-z-Ebene 


Die Randbedingungen langs der Rander I und II bleiben unverandert; lings des 
Randes III haben wir jetzt zu erfiillen: 
aw 


mae OD w= 0; oh =), p=; (9) 


Wir konnen dies folgendermafen erreichen: Zunachst behandeln wir das Problem so, 
als ob die Schalen in der y-z-Ebene biegesteif verbunden waren (wie frither). Die beiden 
Bedingungen wu = 0 und ~~ = 0 werden infolge der Symmetrie des Temperaturfeldes 
 beziiglich der y-z-Ebene automatisch erfiillt. Die Erfiillung der beiden anderen Bedin- 
-gungen kénnen wir erzwingen, indem wir eine zunachst noch unbekannte zusiatzliche 
_ Linienquer- und Schublast aufbringen und deren Verteilung so bestimmen, daB v= 0 
und w=0O wird (Abb. 6). 
Die hier genannten Randbedingungen gelten auch mit guter Naherung fiir drei- 
eckige Stabilisierungsflossen, wie sie haufig an Flugkérpern Verwendung finden. 


IV: Lésung der Differentialgleichungen 


Iterationsmethode fiir die Grundgleichungen 
Eine Lésung durch Iteration ist méglich, wenn wir die Gleichung der Schalenmittel- 
fliche in der Form 


W =cf (x,y) 


schreiben kénnen’. J bedeutet dabei irgendeine charakteristische Langenabmessung 


3 A.A.Nazarov: On the Theory of Thin Shallow-Shells. NACA TM 1426, 1956. (Ubersetzt aus 
Prik. Mat. i. Mekh., Vol. 13, pp. 547— 550, 1949.) 
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der Schale; c ist der bereits eingefiihrte dimensionslose Parameter, der die Flachheit 
der Schale kennzeichnet, Abb. 1. Die Funktion f (x, y) ist dimensionslos; ihr Betrag 
hat ein Maximum in der GroBenordnung von eins. 

Wir entwickeln nun die Unbekannten w und F nach Potenzen des Parameters ¢ 


w=w,tcw,+c? ws+...) 
P2F Act eer a 
Fiihren wir diese Entwicklung in die Grundgleichungen (1) ein und vergleichen Terme 


mit gleichen Potenzen von c, so erhalten wir folgendes System nicht mehr gekoppelter 
Gleichungen 


(Eh VV oe Va \ (10) 
V2 7? w,=—a(1+n) Vf 
1 2 
(E h)- 7? y?P,=— so 
1 @F, ie 
eS a tee 


Da wir bei der Ableitung der Grundgleichungen fiir die flachen Schalen bereits alle 
GréBen vernachlassigt haben, welche in der GroéBenordnung von c? waren, ist es nicht 
sinnvoll, die Iteration iiber den zweiten Schritt hinaus fortzusetzen. 

Die Randbedingungen kénnen ebenfalls nach Potenzen von c geordnet geschrieben 
werden. Zu den Gln. (10) gehéren dann die folgenden: 


Symmetrisches Temperaturfeld: 


Rand I, II Rand III 
YN = Grin 0 
OW, Ow, 
Nnrsi = 0 Nasi = 0 
Nani = Uy => 
und zu den Gln. (11) 
Rand I, II Rand III 
We == () Qn2 =—tas 
OW, =i OW, 
on é 7 
e (13) 
Nns2 = Nns2 = 0 
Nnn2 = Us = 0 
Fir ein schiefsymmetrisches Temperaturfeld haben wir: 
fiir die Gln. (10) 
Rand I, II Rand III 
Mn1 =0 qui = 0 
Co] ns 
ga +2 =o i =o 
(14) 


u% =0 Nnsy =9 | 


v, =0 uy; =0 


_ chungen 


~ Man kann solche eventuell auf ex- 
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und fiir die Gln. (11) 


Rand I, JI Rand ITT 
0 
My 2 aaa) Qn2 ciliee Boas 
OMnsz oe of ow 
Yn2 ae os La Mnna a =—.() (15) 
Rie) nse =0 


Die Gln. (10) sind identisch mit den Gleichungen, welche das Verhalten einer Scheibe 
bzw. Platte unter der Einwirkung eines Temperaturfeldes 7 bzw. I" beschreiben ‘. 
Die zugehérigen Randbedingungen entsprechen im Falle eines symmetrischen Tem- 
peraturfeldes (Rand I, IT) denen eines freien Scheibenrandes bzw. denen eines eingespann- 
ten Plattenrandes; im Falle eines schiefsymmetrischen Temperaturfeldes denen eines 
unverschieblichen Randes einer Scheibe bzw. freien Randes einer Platte. Da wir ein 
beziiglich der vertikalen Symmetrieebene des Fliigels symmetrisches Temperaturfeld 
vorausgesetzt haben, kénnen wir die Bedingungen auf dem Rand III automatisch 
erfiillen, wenn wir statt des bisherigen Bereiches den gesamten Bereich des Fliigels be- 
trachten. Wir erhalten dadurch statt eines gemischten ein einfaches Randwertpro- 
blem. Das gemischte Randwertproblem la8t sich auch bei den Gl. (11) umgehen. Wenn 
wir ebenso verfahren wie bei den Gln. (10), werden wir langs des Randes III aus Sym- 
metriegriinden ¢,2—0 erhalten. Wir brauchen also nur noch lings dieses Randes eine 
zusatzliche Linienquerlast von der GréBe 


p= — 21 ans (16) 

aufbringen. 

Auffinden von Partikularintegralen 

Fiir die Gln. (10) suchen wir zunachst Partikularintegrale und lésen dann die zu- 
gehorigen homogenen Gleichun- 
gen mit entsprechend modifizier- 
ten Randbedingungen. Einfache 
Partikularintegrale erhalten wir 
durch Lésen der folgenden Glei- 


Vy? F,=—EhaT | 


V2 w,=—a(l+u)T. f a 


perimentellem oder numerischem 
Wege gewinnen. Eine andere Még- 
lichkeit besteht darin, das Tem- 
peraturfeld in eine Reihe zu ent- 
wickeln. Wir wollen dies hier 
durchfiihren. Fiir die Reihe verwenden wir ein Koordinatensystem 7, g, Abb. 7. 


Abb. 7 


T 7, ¢) =D Pan?™ 008 ~" | 
ae Ga, n= 0, 1; 2,3, ...).4 (18) 
oe Wome Ly 7 = NM — 
P(r, ¢) = 2/2) Gam ™ CO8 —— | 


4 #.Melan und H.Parkus: Warmespannungen infolge stationdarer Temperaturfelder, 
S.56—70. Wien: Springer-Verlag. 1953. 
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Die Koeffizienten p,m Und drm kénnen noch zeitabhangig sein. Das Temperaturfeld 
haben wir beziiglich der y-z-Ebene symmetrisch vorausgesetzt, daher enthalten die 
Reihen keine Sinus-Glieder. Wegen der Linearitét der Grundgleichungen gentgt es, 
von den Reihen nur je ein Glied zu betrachten. 


T= pr" €0s “2 Y 
F=qr" cos “s Q. 
Fir das Partikularintegral machen wir den Ansatz 
Fye=pr"+* cos 9. (19) 


Tragt man diesen Ansatz in die erste der Gln. (17) ein und fiihrt einen Koeffizienten- 
vergleich durch, so erhilt man 


k=32 

NY akhp 

Bt ne 72 
(ms 2) A 

und fiir einen analogen Ansatz fiir w,’ 

Kea 2 

ey a(l+ )q 

U Me ne 72 
Ce he earaar 


2 72 


Fiir den Sonderfall (m+ 2)? = “* machen wir fiir Ff,’ den Ansatz 


a 
= Tr NT — 
Feapret* In| cos @, 


wobei ¢ eine charakteristische Lange ist. Wir wahlten dafiir die Fliigeltiefe, Abb. 7. 
Einsetzen in Gl. (17) und Koeffizientenvergleich ergibt 


bare 

— akhp 

PF (na 24 
und analog fiir w,’ 

k=2 

Ai _a(l+p)q 

1 2 (m+ 2) t 


Lésung der homogenen Differentialgleichungen 
Wir miissen nun die homogenen Differentialgleichungen lésen 


“vet =o 
ae (20) 
Vy? V7 =O 
ripe (21) 


Es sind dies die bekannten Bipotentialgleichungen. Wir wollen zu ihrer Lésung die 
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Methode von Muschelischwili heranziehen®. An den Randern haben wir jetzt die 


um die entsprechenden Werte aus dem Partikularintegral verminderten Werte zu 
nehmen; z. B. 


Dye —— =e 
ow,! dw? 
dig tha evden 


Die Lésungen der Gln. (20) lauten in komplexer Schreibweise 


2P i =z D, (z) +2 D, (2) +X, (z) + XY (2) 


2wy"=2 9, (2) +29, (2) + %1 (2) + x1 (2) 


Dabei sind 9, (z), X, (z), ~, (2), % (z) analytische Funktionen der komplexen Variablen 
z7=x-+17y. Ferner ist z= x—iy. 
Spater schreiben wir ¥ (z) fiir X’(z) und p (z) fiir y’(z). Wir miissen nun die Rand- 


bedingungen in den Funktionen © und Y bzw. @ und y formulieren. Fiir die Gln. (20) 
haben sie die folgende Form® 


(22) 


a) im Falle eines symmetrischen Temperaturfeldes 


® 2) +20 @+¥@)=—(F +iF) (23) 
e@+2¢ @+ v@=—(f tise): (24) 


b) im Falle eines schiefsymmetrischen Temperaturfeldes 


— ® (z)+2@' (z) +P (2) =2G4 (ut+iv). (25) 
—Ag (2) +z (2) +y (= 
1 x(s),9(s) s : . 
=—xaam J[— me tif (a+ Ge) ds | de +a ay). (26) 
one ny 


: : : 3 
Dabei sind w und v die Verschiebungen in Richtung x und vy; ferner ist x = me , und 
4. = it . Die Groen m’ und (¢ + * i‘ lassen sich durch w? ausdriicken. 


dMn 
Os 


Fiir die Gln. (21) lauten die Randbedingungen in den komplexen Funktionen ebenso 
wie fiir die Gln. (20), nur ist bei einem schiefsymmetrischen Temperaturfeld 


OMns2 ’ 


(a 2 — as 


zu ersetzen durch 


OMns2 \P 0 
(du2+ = 4 +] pone 


Da wegen des Dreiecksrandes wenig Aussicht besteht, Losungen in den Koordinaten 
x, y zu finden, bilden wir den Dreiecksbereich (z-Ebene) konform auf das Innere des 
Finheitskreises (¢-Ebene) mittels der Funktion 


z= (C) 


5 N. I. Muskhelischvili: Some Basic Problems of the Mathematical Theory of Elasticity. 
Groningen (Holland): P. Noordhoff Ltd. 1953. (Englische Ubersetzung durch J. R. M. Radok.) 

6 §. G. Lechnizki: Von einigen mit der Biegung diinner Platten verbundenen Fragen. Prik. 
Mat. i. Mekh., II, 2, S. 181—209 (1938). 
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ab. Die Variable ¢=oe!® auf dem Einheitskreis bezeichnen wir mit o =e'*. Die Rand- 
bedingungen lauten dann, z. B. fiir w bei symmetrischem Temperaturfeld 


 (0)-+ 22 g(a) + ¥(2)=— (+i) = (0). (27) 


Dabei ware statt y (c) praziser g (oc) = [w (c)] zu schreiben. Der Index 0 kann aber 
weggelassen werden, da eine Verwechslung kaum méglich ist. Die anderen Randbedin- 


gungen findet man in analoger Weise. Die Ausdriicke le und G (co) entwickeln wir 
: @ (o 


ebenso wie die unbekannten analytischen Funktionen in Potenzreihen von o 


<r \? 7 
CAUSE Dar G4(ao=D) Ci. vp(o)= Lar yl(o)= ib se% (28) 
w’ (c) — co — co 0 0 


Setzt man obige Reihen in die Randbedingung (27) ein und fiihrt einen Koeffizienten- 
vergleich durch, so erhalt man ein System von unendlich vielen gekoppelten Gleichun- 
gen fiir die unbekannten Koeffizienten a, und b,. Wir wollen jedoch, nach Musche- 
lischwili, die Gleichungen fiir g und y entkoppeln. Zunachst multiplizieren wir beide 


Seiten der Gl. (27) mit s-7 ; —. und integrieren dann iiber do; wir erhalten so 


Se OR ed Fae, do mary pt do ; 

Ot sa OS Ee (0) sae + ¥ = a7 OE (ge: (29) 

Wenn wir in diese Gleichung die Potenzreihen fiir » (c), ae und G (oc) einfiihren und 
@ (oO 


die Integrale ausrechnen, erhalten wir ein System von linearen Gleichungen, in denen 
als Unbekannte nur mehr die Koeffizienten a, vorkommen 


Vorausgesetzt, da es sich um gleichmafig konvergente Reihen handelt, konnen wir 
schreiben 
ay Oy CP S (Sin coes-1} C+y(0)=D Oe. (30) 
0 n=0\m=1 0 
Durch Koeffizientenvergleich erhalten wir zwei Systeme reeller linearer Gleichungen 
fiir Real- und Imaginarteil der a,. Wir wollen mit einer Naherungslésung zufrieden sein, 
die wir dadurch erhalten, daB wir nur eine endliche Anzahl 2n + 1 von Koeffizienten c; 
und C; beriicksichtigen. Wir erhalten dann n Koeffizienten a;. Die Gleichungen fiir 
diese lauten 


Oy Ey Cy 2G, yee 3h Sone ee ee Hide C=C, 
Os = Gy Cy 2a Ogee. ene ae (71 —1) NGn-10,=C, 
Hida hivie Usteks Sed Ls dh SRW a (31) 
Oy + Ola eau ad OS dee ee eee e=(6)- 


Fir y(0) haben wir die Gleichung 


Oy 0g 2dnty heasecdow veg nein NAn Cn—1 + y (0) =O. (32) 
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Zur Ermittlung der Koeffizienten der Funktion y (¢) gehen wir von der konjugiert kom- 


_plexen Form von Gl. (27) aus 


_ Wir multiplizieren wieder mit 


a+ 20 


9 (3) + Grey % (0) + (0) =G (0). 


1 1 ? ; x 
Fe und integrieren iiber do 


co 


PY + 2"( ahaa ee mas] abe ge!) GR STOLE 
b= m=1 (=) fen) 


_ Daraus erhalten wir die unbekannten Koeffizienten b;, 


ay 


As 


b, =O_, aa Cp ek M Am. (33) 


1 


_ Um die Lésung eindeutig zu machen, wurde ¢ (0) = 0 festgelegt. Nun kénnen wir w* 
~ anschreiben 


2wt =o (LC) (C) +0 (6) y (+74 ( +4 (8), 
wobei 7 (2)= | y (0) 0" ()dt +0. 


Im allgemeinen werden jedoch nicht so sehr die Verschiebungen w als die SchnittgréBen 
interessieren. Fiir diese erhalt man aus w* in einem den Radien und Kreisen der ¢- 
Ebene entsprechenden krummlinigen orthogonalen Koordinatennetz der z-Ebene® 


(34) 


(m, + mg) = 2N (1+) [p* (6) + 9* (2)] 
2N 
0 


Cat grins atacenirta hi Teall 


= 94 h 
(my —m, + 20m,5)'= 2aTD 


wobei 
yi (omg tind, yh (Sea a 


Die Indizes o und # bezeichnen die Richtungen zunehmender Werte von @ und #. 
Fir die Querkrafte erhalt man 


: —4N¢ ; 
(Go— 190 )* = ae (¢). (35) 


Ganz analog kann man beim ersten Schritt der Iteration fiir das Scheibenproblem vor- 
gehen. Mit einer analogen Bezeichnungsweise ergibt sich fiir die Schnittkrafte 


(m+ ng)! = 2 [B* (6) + O* (2)] ne 
(5 — Np + 2iNgs)"= Seo [wm (£) B* () +o’ (6) B* (2). 


Fiir die Verschiebungen uw, v erhalt man 


2G (u+iv)t=x D (2) —w (¢) O* (¢)— ¥ (6). (37) 


Da die Partikularintegrale fiir w und F bei symmetrischem Temperaturfeld ebenso wie 
die Randbedingungen bis auf einen Faktor gleich sind, kénnen wir in diesem Fall die 


Side Sao ise AMM | °. 
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Liésung F* aus w* durch Multiplikation mit einem Faktor gewinnen. Bei schiefsymme- 
trischem Temperaturfeld sind die Randbedingungen fiir w und f verschieden. 

Die Gleichungssysteme, die dann zur Bestimmung der Koeffizienten der Funk- 
tionen y (C), y (¢), ® (6), W (¢) aufzulésen sind, unterscheiden sich von den Gln. (31) 
nur dadurch, da8 jeweils der erste Summand auf der linken Seite den Faktor — x bzw. 
— A aufweist. 

Haben wir auf diese Weise alle interessierenden GroBen ermittelt, konnen wir zum 
zweiten Schritt der Iteration tibergehen. 


Der zweite Schritt der Iteration 


Wir brauchen zuerst Partikularintegrale der Gln. (11) 


‘ 2 
(Eh) 72 7? Fy = = a 
1 e238; 
IN 2 V2 Wn ae oe So 
wobei jetzt w,, F; und p bereits bekannt sind. Aus w, und F’, bilden wir die Ausdriicke 
1 2 1 #&F ' he 
= =a atid ae Diese entwickeln wir in Reihen von der Form 


1 CPi de & VW 


r or <= 
m 


al =. NM TG — 
ee Pant” COS — op. 


Fir das Partikularintegral F,’ machen wir wieder den Ansatz 
FP = Qn T™+* COS “= P- 


Einsetzen in die Gln. (11) und Koeffizientenvergleich liefert 


k=4, Qna= a E we 
(m +2)?" 2) |(m +42 == 
oe a 
Fy Vollkkommen analog erhalten wir ein Partikular- 
integral fiir w,. Fiir den von p herriihrenden Anteil 
von w,? kénnen wir z. B. die Lésung fiir den frei 
drehbar gelagerten Plattenstreifen mit Linien- 
belastung verwenden’, Abb. 8. Wir miissen zuerst 
p(y)in eine Fourierreihe von der Periode 2 ¢ ent- 
= wickeln 
Abb. 8 7 
p(y)= Sa, sin “72 
Der von p herriihrende Anteil von w,? lautet dann : 
. ant nee) ==. 
fiir y = 0 wit = Saag (1+ ; Je t sin “= Y, 
we eer Qn 8 nm y pee : 
fiir y <0 w= tN iw ( ee ‘sin y 


7K. Girkmann: Flachentragwerke. 8. 170—173. Wien: Springer-Verlag. 1948. 2. Auflage 
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Nun kénnen wir lings der Rander I und IT wieder die zur Lésung der homogenen 
Gleichungen erforderlichen Ausdriicke bilden und die Gin. (21) mit den entsprechenden 
Randbedingungen lésen. Dazu kénnen wir wieder die Methode von Muschelischwili 


verwenden. Der in den Randbedingungen vorkommende Ausdruck 17! kann aus der 
Dreiecksform und dem Parameter c ermittelt werden. 

Wir miissen nun fiir alle Koeffizienten Pn, drm der Reihe des Temperaturfeldes 
Finm, Fenm Winm, Woanm bestimmen. Dann kénnen wir die Losung fiir Fliigel aller Dicken 
anschreiben. Infolge der Linearitiat der Gleichungen lassen sich aus den einzelnen Lésun- 
gen die Koeffizienten der Reihe des Temperaturfeldes herausheben, z. B. 


Tica = Pnm Y ER ge 
Die Losung wird dann 


dies pa (Pam Fs ima™ + Gam Dain | 


, (38) 
Uo S= oF (Gna Whans sive C Pam Werane ): 


mn 


Die GroBen F1,,* usw. lassen sich fiir eine bestimmte Dreiecksform ein fiir allemal 
ermitteln. 


Zur Frage der konformen Abbildung 
Zur Losung des Problems brauchen wir die Abbildungsfunktion 


z= (C). 


die das Innere des Dreiecks (z-Ebene) konform auf das Innere des Einheitskreises 
(¢-Ebene) abbildet. Ferner brauchen wir zur Lésung des Gleichungssystems Gln. (31) 
den Ausdruck 5 
WM \o 
Um die Abbildungsfunktion zu finden, konnen wir von der Schwarz-Christoffelschen 
Abbildungsformel ausgehen, die das Innere polygonal begrenzter Bereiche auf das 
Innere des Kinheitskreises abbildet®. 


Abb. 9 


eee oO eT (C—a)** i=1, 2,3,...1; (39) 


dabei bedeuten n die Anzahl der Ecken des Polygons, a; die Bildpunkte der Ecken auf 
dem KEinheitskreis. 


8 A. Betz: Konforme Abbildung, S. 223— 229. Berlin-Géttingen-Heidelberg: Springer-Verlag. 
1948. 
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Ferner gilt 


Cp = Bi; —, 
Wir wollen von vornherein folgende Spezialisierungen festlegen: 
Leaps d. h. wir beschranken uns auf gleichschenkelige 
Dreiecke 
2. @,=1, 
dg=— 3 +3 / 3%, 
a,=—}—F 31 


Bekanntlich kann man, um die Abbildung eindeutig zu machen, die Zuordnung dreier 
Punkte am Rand oder eines Punktes am Rand und eines Punktes im Innern vor- 


geben. 
Die Schwarz-Christoffelsche Formel liefert 


c 
z= (C)=d,+d, {(¢— 1)® (¢+4—4)/37) (¢+3+4)/32) de. 

0 
Die Integrale lassen sich nur fiir ganz spezielle Werte von «, «, geschlossen angeben 
(z. B. fiir das gleichschenkelig-rechtwinkelige Dreieck®) und fiihren dann auf elliptische 
Integrale. Die Abbildungsfunktion kann man aber auch so gewinnen, dafS man den 
Integranden in eine Potenzreihe von ¢ entwickelt und gliedweise integriert. Fir das 
gleichschenkelig-rechtwinkelige Dreieck erhalt man auf diese Weise 


z—dy _ se Ghee 461 .. 
en oS gS — 39S + 596 io240 6° t+: 


Auf dem Einheitskreis konvergiert die Reihe, wie man sieht, schlecht. Bricht man sie 
nach einem bestimmten Glied ab, so stellt das verbleibende Polynom nicht die Abbil- 
dung des Dreiecks, sondern einer ,,drei- 
ecksihnlichen“ Figur auf den Einheits- 
kreis dar. Welche Randfigur anstatt 
eines rechtwinkelig-gleichschenkeligen 
Dreiecks bei Beriicksichtigung der 
Glieder bis einschlieBlich ¢ entsteht, 
zeigt Abb. 10. Da das Polynom das 
Dreieck nur sehr schlecht annahert und 
es auBerdem sehr miihsam ist, aus der 
Schwarz-Christoffelschen Formel noch 
weitere Koeffizienten zu ermitteln, er- 
Abb. 10 scheint es ratsam, sich die konforme 
Abbildung auf experimentellem Wege 

zu beschaffen; etwa mit Hilfe des elektrolytischen Troges?® 11, 


Wir brauchen noch den Ausdruck ae . Da @ (c) nur numerisch (aus den Messun- 


A. Scherer: Allgemeine Herleitung singularer Lésungen der biharmonischen Gleichung 
sowie Aufstellung Greenscher Funktionen am Beispiel einer Dreiecksplatte mit freiem Rand. 
Ing.-Arch. XXV, 8. 255—272 (1957). 

Pp. A. Kennedy and K. Gordon: Electrolytic Tank, design and applications. Rev. of Scient. 
Instr. 27, No. 17, pp. 404—413 (1956). 

4“ C. Hachmeister: The application of the electrolytic tank to the solution of potential 
field problems. Annals of the New York Academy of Science 60, Art. 6, pp. 937—947 (1955). 
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gen) vorliegt, fihrt die Bestimmung von w’ (a) durch Differenzieren zu groBen Unge- 
nauigkeiten. Wir kénnen jedoch w’ (c) auch anders gewinnen, wenn wir uns folgendes 
_tiberlegen. Ist V die zum gemessenen Potential U konjugierte Potentialfunktion, so gilt 
d 2 d ; di d : ; 

q¢ WU +iV) =] (U+iV1 =] [UV +i Jo’ (0). 


dz 


- Nach den Cauchy-Riemannschen Gleichungen ist aber 


d ‘ oU .dU 
qe lU+t1V) = a ty 
Langs des Randes wird, Abb. 11, 
0U ella aU as aU an .[ aU as aU | 
We ag aeeee Ghee de oo oe by cen Oy | 


Da der Rand eine Potentiallinie darstellt, gilt 


also gemaB Abb. 11 


an 2 0U eu , . ' 
Feo ay gn (Sin «+ 2 COS a). 


Analog ist am Kreisrand 


Abb. 11 


ee [U+iV]= (cos #— 7 sin #) 


~Damit wird: 


(=) (cos # — isin #) ) i. (ee 

ar }R P'S ten tlamae 40 
o (OS 7} ee cree bus F F -) 
oe (sin « + 7 Cos a) an 


w (co) kann man der Zeichnung entnehmen, 
Z=ol(o)=h-e?: 


Hierbei sind h, O die Polarkoordinaten der Randpunkte. Der gesuchte Ausdruck lautet 
dann 


au 
8 (:) LC i(a+O-0-+ a 
0’ (0) (3) 
Or /R 
( o) ist eine Konstante, ae mu am Dreiecksrande gemessen werden. 
or JR 7 
Wir brauchen den Ausdruck Gl. (41) in Form einer Potenzreihe 
oo = ye “ (42) 
S die 


Um die Koeffizienten c, zu bestimmen, machen wir eine harmonische Analyse von 
Real- und Imaginarteil der Gl. (42) 
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|=« | ¢, cos?+ cy cos 20+. 
+¢, sind?+c¢, sin 20+... 


| =¢) +c,’ cos?+c¢,'’ cos20+... 
to,’snd+ec,”sn20+... 
Es muB gelten 
Cn O-" +6, 0" = C, COSNY+ C, SIN NI+74 (c,’ cos nF -+ Cc,’ sin n #). 


Daraus ergibt sich 
L( Cn Cn’) +4 (Cp — Cn) 
[( Gn —Cn'") +4 (Cn' > En). 


S 


1? 


| n #0 (43) 


V. Ein Beispiel 


AbschlieBend wollen wir noch ein einfaches Beispiel behandeln. Fiir einen Delta- 
fliigel mit gleichschenkelig-rechtwinkeligem GrundriB soll der erste Schritt der Ite- 
ration durchgefiihrt werden. Fiir das Temperaturfeld nehmen wir ein Glied der Reihen- 
entwicklung Gl. (18) mit m=1, n=1. Ferner setzen wir voraus, da8 das Temperatur- 


feld symmetrisch ist. Mit «= n (gleichschenkelig-rechtwinkeliges Dreieck) haben wir 


oe pa (44) 

=q r cos 4Q. J 

Das Temperaturfeld ist in Abb. 12 durch 

die Temperaturverteilung am Rand ange- 

deutet. 

. Wir brauchen zuerst Partikularinte- 
grale der Gln. (10). Diese lauten, Gl. (19) 


Ehp — Be 
Fe —* P cos 49 


= 
‘ 


(45) 


»_ edtu)g 
ae 7 


w r> cos 4. 


Wir betrachten zunachst die zweite Gleichung (Plattenproblem). Zur Lésung der 
homogenen Differentialgleichung brauchen wir langs des Randes den Ausdruck 


Diesen kann man aus dem Partikularintegral leicht ermitteln. Es ist 


Owe awP ia 1 dw = 

aE Fe Gis oe - ris COS @ | 
awh Owe = Le Ye TS 
“Sar & op 88 YT = ae sin —p. 


_ Fiir die weitere Rechnung brauchen wir die konforme Abbildung des gleichschenke- 
lig-rechtwinkeligen Dreiecks auf den Einheitskreis. In Abb. 13 ist diese durch die 


Zuordnung zweier orthogonaler Netze angegeben. Die MaBstabe sind so gewahlt, daB 
in der Umgebung des Nullpunktes gilt 


|2|=|¢]. 
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Die Ermittlung erfolgte mit Hilfe eines elektrolytischen Troges *. 


Um die Koeffizienten C, zu gewinnen, zeichnen wir unter Beniitzung der Zuordnung 


der Randpunkte die Funktionen = und = in Abhangigkeit von ? auf (Abb. 14 und 
Abb. 15) und machen harmonische Analysen 
= =O+C, cos 3+ CO, cos 20+... 
+0, sind+C, sin 20+... 
a =C,'+C;' cosd+C,' cos 20+... 


+C,"smd0+C," sn20+... 
| Daraus ergibt sich 
ow ow 


G(o)=—(F +i 


1 OPE Ore Ua... 
+0_,01+0C_,07+4+... 
) mit 
Oo=— [Co+i Oy’) i 
een Cs-F Or (Cy SO 
Ca —4 [(C,—C,") +2 (C, +6,)]. 
Wir wollen mit einer groben Naherungslosung zufrieden sein und die Reihe bereits 


nach dem sechsten Glied abbrechen. Die Koeffizienten lauten, abgesehen von einem 
dimensionslosen Faktor 0,0744 a (1+) q#?, der herausgehoben wurde: 


Co = +1 0,86 G,=—1 0,22 

C,=+10,51 C,=+70,12 a, 
C,= —i 0,08 C,= +4 0,07 

C,= +2 0,29 | 


* Die Abbildung wurde mir von den Herren Dr. R. Beer und Dipl.-Ing. H. Peterschinegg 
zur Verfiigung gestellt, wofir ich ihnen auch hier meinen herzlichen Dank ausspreche. 
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C_,=—1 8,20 C_,= —1 0,33 | 
C_,= +7 2,76 C_;= +1 0,97 (47) 
RUE Te 1 (he eat OO | 

Wir brauchen dann noch die Potenzreihen von oe Abb. 16 zeigt den Real- und 


Abb. 17 den Imaginarteil von a aus Messungen im elektrolytischen Trog, unter 
@ (Oo 


D042 (1+) gle 


Abb. 15 


Beniitzung von Gl. (41) ermittelt. Bei Abb. 14 bis 17 stellen die strichlierten Linien die 
Annaherung der entsprechenden Kurven durch die nach dem sechsten Glied abgebro- 
chenen Fourierreihen dar. Fiir die Koeffizienten c, und c_, erhalt man 


C= —0,117 C, = — 0,080 
Cc, = — 0,856 ¢; = + 0,021 , 
C= + 0,061 C, = — 0,009 (48) 


C3 = — 0,006 
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C_,= — 0,252 C_4= + 0,036 
C_»= +0,655 c_,=+0,101 (49) 
C_3= +0,108 c_,= + 0,019. 


Abb. 17 


Mit den Zahlenwerten Gl. (46), (47) und (48) erhalten wir fiir das Gleichungssystem (31) 


+ 1,856a, — 0,122a, + 0,018a, + 0,320a, — 0,105a, + 0,054a, = + 0,511 


— 0,061, + 0,012a,+ 0,240a,—0,084a,+0,045a, 0 = —0,08% 
+ 0,006a, + 0,160a,-+ 0,937a,+0,036a, 0 0 ..=-1.0,293 ab 
+ 0,080a, — 0,042a,-+0,027a,+1,000a, 0 6 = — 0225 
— 0,021a, + 0,018.4, 0 0° +£000a, 0 =+0,12: 


+ 0,009, 0 0 0 0 +1,000a,=+-0,07i 
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Dabei wurde angenommen, da® die Gré8en a, rein imaginar sind; dann gilt an = — On: 
Dies trifft deshalb zu, weil die SchnittgréBen beziiglich der y-Achse symmetrisch sein 
miissen. Das ist aber im Hinblick auf die Gln. (34) bis (37) und Gl. (33) nur dann der 
Fall, wenn gy (¢) und y (¢) langs der y-Achse rein imaginar sind. o (¢) und w’ (¢) sind 
dort, wie man aus Abb. 13 ersieht, ebenfalls rein imaginar. 


Die Lésung ergibt 1 


@,= +0,3114 @,=-— 0,261% 
d,= — 0,170% ds = +0,130% (51) 
dz = + 0,3472 Ag = + 0,067. 


Aus der Gl. (32) erhalt man 
yp (0) =b,= — 1,140%. (52) 


Fiir die Koeffizienten der Funktion y (¢), 6, ergibt sich mit Gl. (33) 


b, = +4,228% b, = —0,395% 

b= — 3,6591 b= + 0,0193 

b, = + 1,960i b, = —0,0601 is 
b, = — 0,169% bias = 00582 

b, = —-0,207i b= — 0,008% 

b, = +0,409% 


Abb. 18 


Nun wollen wir noch einige SchnittgréBen lings der Symmetrieachse des Fliigels 
ermitteln. Die Biegemomente erhalt man aus den Gln. (34). Zu diesen Momenten sind 
die vom Partikularintegral herriihrenden zu addieren. 


Die Momente langs der Symmetrieachse sind in Abb. 18 dargestellt. Die Funktionen 
@ (¢) und » (¢) kann man ebenfalls mit den Koeffizienten Gl. (51) bis (53) bilden, da ja 
die Randwertprobleme fiir w und F identisch sind. Statt des Faktors 0,0744 a (1+ 4) qe? 
miissen wir aber jetzt setzen: 0,0744 a Eh pt?. Die Schnittkrafte erhalt man aus GI. (36) 


© R. Zurmiht: Matrizen, 8. 248— 258. Berlin-Géttingen-Heidelberg: Springer-Verlag. 1950. 


(Hingegangen am 20. September 1960) 
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Tafeln fiir die Berechnung des elliptischen Paraboloids iiber 
rechteckigem Grundrif in der Membrantheorie 


Von M. Soare, aaror 


Zusammenfassung. Die Losung fiir den Membranspannungszustand in einem elliptischen Para- 
_ boloid tiber rechteckigem Grundri8 wird unter Annahme einer parabolischen Vertikalbelastung in 
analytisch-tabellarischer Form gegeben. 


Das elliptische Paraboloid als typisches Beispiel einer Schale mit positiver Gaub- 
scher Kriimmung ist mit Vorliebe durch verschiedene Forscher studiert worden, dank 
der analytisch zuginglichen Lésungsméglichkeiten. 

Der Zweck der vorliegenden Arbeit ist, eine Reihe von neuen Tabellen zur Berech- 
_ nung der Membranschnittkrafte infolge einer nach einem parabolischen Gesetz verteil- 
ten Belastung zu geben. 

1. Es sei 


gfe ede tee ailfiofes 0) (1) 


die Gleichung der Mittelflache, wo 2a - 2b die Abmessungen des rechteckigen Grund- 
) risses und f,, f, die Pfeile der Binderscheiben darstellen. 
Die Belastung pro Grundri®flicheneinheit lautet: 


pee 
4=g(1—4). (2) 
Die Grundri®projektionskrafte lassen sich bequem schreiben*, so wie sie bereits 
; von Pucher! abgeleitet wurden, mittels einer Spannungsfunktion /’; diese mu8 einer 


Differentialgleichung zweiter Ordnung geniigen, welche fiir die angenommene Mittel- 
flache und Belastung lautet 


2f, @F 27, &F a 
be” Ga? a oy (1 i) (3) 
Als Randbedingungen werden vorgeschrieben: 
a er : = ar 
inw= +4a:n,=0, alsop;=90; in y= +6; n, =0, also 7s = 0. (4) 


2. Eine Partikularlosung der Gleichung (3), die den Randbedingungen langs der 
Seiten « = + a geniigt, ist 


ga? b? ae 4a" i ( a ( = 7) 
F) > 124i, (5 6 a a oat, 1 aa 5 aE (5) 
' oder als einfach trigonometrische Reihe mit konstanten Koeffizienten geschrieben 
| ae 
64gaeb? ,(—-1)> sn we m 
FQ) = a m = COS 9, (= INS) (6) 
Die Lésung der homogenen Gleichung 
2 fe Fa) 2h, - 2 Fy) 0 
Fl oe a2 ayz 


kann als eine einfach trigonometrische Reihe mit veranderlichen Koeffizienten ange- 
setzt werden: 


Fa =D) (Anch Bry + Br sh Bry) Cos 


n 


Nx 
2a 


* Die Bezeichnungen fiir die Schalenkrafte sind die gleichen wie bei A. Pucher’, E. Tung]® 
und K. Girkmann’? (S. 390). 
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in der die Beziehung a 
Bn = m %/ to (7) 
eingefiihrt wurde. 


Aus Symmetriegriinden folgt B, = 0; aus der zweiten Bedingung (4) erhalt man 
leicht die Werte der Konstanten A, 


n—1 
64ga?b?  (—1) - 
7 fe n> ch Bnb * 


A, = 


Mit der Partikularlésung (5) kann man schlieBlich die Spannungsfunktion schreiben : 


n—1 


ga*b? | 1 / Ae oo 64 ., (—1) 4 ohBny | NUL 
F= (1 )(5—4) fa Ps, a chB, b 8 : (8) 


(j= Ty 3 Gace) 
Nach zweimaliger Ableitung folgen die Ausdriicke fiir die GrundriBprojektionskrafte : 


n—1 
er —_ 16ga? Sai eh Bag NUT ga 
v F 3 cs. UPd, 
ag ahh a ma hath iain. SEAia eaialy 
n—1 
er es et wy) tog) (=) = eh Bny M1 % gb? | (9) 
a pons pie he Bete = NS : . == : 
oa ae eS 2 fo i z) Nets - n3 ch Bnb se a ks tae { 
n—1 
er = lé6gab (1)? SIEDAY. IG L gab 
een rs = — ——- \ . . —_— —fp ° 
dx dy xe Vi hy ne chBnb ea Pye Vi 


Die Ausdriicke der GrundriBprojektionskrafte fiir einen gegebenen Punkt (z/a, 
y/b) hangen nur von dem Verhiltnis /,/f, ab, was ihre Tabulierung erméglicht. 
Die Schnittkrafte in der Schale werden durch die Beziehungen 


yy Ge 

Sie ele re 

b2 
8, = — hyd, (10) 

b 

PS oe 

mits 

ausgedriickt, wobei 

a= V(01+-8*) [4 244) (11) 


und die Koeffizienten k,, k., kj. den Tafeln 1, 2, 3 entnommen werden kénnen. 
Zwischen den Koeffizienten k, und k, besteht die algebraische Beziehung 
1 le 
hy + ky = 5(1— Gh. (12) 
Die Lésungen fiir die Schnittkrafte sind in den Ecken regular (haben endliche Werte), 
da die Last in diesen Punkten null ist. 
3. Die Loésung der Gleichung 


2 fe ‘ oF Pee Gale y? 
b2 ax wee aye Sie (1 — a 


(13) 
folgt aus der Lésung der Gleichung (3) durch Austausch von 2, a, f, bzw. y, b, fe. 
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Ks sei noch erwahnt, daB die Lésung der Gleichung 


of, @F of @F 
ia toy oF riper s (14) 


entsprechend einer konstanten Grundrifbelastung in der Fachliteratur ausfiihrlich 
behandelt worden ist. 


Auf Grund des Uberlagerungsprinzips kann man aus den drei vorhergehenden 
Fallen den von der Belastung 


Z= (Cy) + C7 


hervorgerufenen Spannungszustand bestimmen, wo C,, O,, C, beliebige Konstanten 
| sind. 

| Wenn man die Lésungen der Gleichungen (14), (13) und (3) durch die Indizes 0, 1, 2 
) bezeichnet, lauten Belastung und Spannungsfunktion 


Zi == lig — 2, — Ls, 


a 
oe 


+ 0,4) (15) 


(16) 
1 Ata oa ue 


2 \ 


02 
Z=9(C,+C,+4+ C,), Z,=90,(1—~3), Z,= 90, (i—$). (17) 


Eine andere Lésung fiir dasselbe Problem wurde von E. Tung] * * gegeben, wo die 
Belastung (15) nur in zwei Teile aufgespalten wird. 
4, Fiir den Sonderfall, daf8 die Konstanten Cy, C;, C, in (15) durch die Ausdriicke 


Co= 0, C*; es 6,2) 


| gegeben sind (/ eine beliebige Lange), lassen sich die Schnittkrafte in geschlossener Form 
| wie folgt ausdriicken: 


s,=—3f, (1-4), 


“7204 a 
rat 
pat 
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Tabelle 1. Werte des k,-Koeffizienten 


y/d 
x/a hilfe as 
0 0,25 0,50 0,75 1 
————— —————— 
1/0,4 0,334.20 0,344.10 0,374.19 0,425.59 
1/0,6 0,275.92 0,288.41 0,326.76 0,393.47 
1/0,8 0,238.31 | 0,252.95 0,298.02 0,378.14 
0 1,0 0,205.32 0,221.12 0,270.67 0,360.26 0,500.00 
0,8 0,172.87 0,189.63 | 0,242.86 0,341.47 
0,6 0,133.43 0,150.91 0,207.35 0,316.21 
0,4 0,085.49 0,102.47 0,159.92 0,279.26 
1/0,4 0,309.74 0,319.13 0,347.77 0,396.96 
1/0,6 0,255.43 0,267.18 | 0,303.37 0,366.77 
1/0,8 0,220.47 0,234.16 0,276.47 0,352.23 
0,25 1,0 0,189.87 0,204.61 | 0,250.96 0,335.45 0,468.75 
0,8 0,159.82 0,175.40 0,225.06 0,317.78 
0,6 0,123.32 0,139.54 0,192.03 0,294.04 
0,4 0,079.00 0,094.70 0,147.99 0,259.40 
1/0,4 0,238.97 0,246.70 0,270.60 0,312.23 
1/0,6 0,196.48 0,205.90 0,235.24 0,287.84 
1/0,8 0,169.30 0,180.12 0,213.89 0,275.82 | 
0,50 1,0 0,145.67 0,157.21 0,193.86 0,262.25 0,375.00 
0,8 0,122.52 0,134.64 0,173.61 0,248.05 
0,6 0,094.47 0,107.00 0,147.87 0,229.04 
0,4 0,060.49 0,072.56 0,113.73 0,201.40 
1/0,4 0,130.45 0,135.01 0,149.31 0,175.60 
1/0,6 0,106.90 0,112.27 0,129.29 0,161.22 
1/0,8 0,091.94 0,097.90 0,117.11 0,153.79 
0,75 1,0 0,078.73 0,084.89 0,106.02 0,145.91 0,218.75 
0,8 0,066.41 0,073.09 0,094.77 0,137.67 
0,6 0,051.18 0,058.03 0,080.56 0,126.65 
0,4 0,032.77 0,039.31 0,061.80 0,110.86 
1 ke 0 0 0 0 0 
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Tabelle 2. Werte des k,-Koeffizienten 


y/b 
x/a filfe 
0 0,25 0,50 0,75 | 1 
a omer ee ee ey ad | er A 
1/0,4 0,165.80 0,155.90 0,125.81 0,074.59 
1/0,6 0,224.08 0,211.59 0,173.24 0,106.53 
1/0,8 0,261.69 0,247.05 0,201.98 0,121.86 
0 1,0 0,294.68 0,278.88 0,229.33 0,139.74 0 
0,8 0,327.13 0,310.37 0,257.14 0,158.53 
0,6 0,366.57 0,349.09 0,292.65 0,183.79 
0,4 0,414.51 0,397.53 0,340.08 0,220.74 
1/0,4 0,159.01 0,149.62 0,120.98 0,071.79 
1/0,6 0,213.32 0,201.57 0,165.38 0,101.98 
1/0,8 0,248.28 0,234.59 0,192.28 0,116.52 
0,25 1,0 0,278.88 0,264.14 0,217.79 0,133.30 0) 
0,8 0,308.93 0,293.35 0,243.69 0,150.97 
0,6 0,345.43 0,329.21 0,276.72 0,174.71 
0,4 0,389.75 0,374.05 0,320.76 0,209.35 
1/0,4 0,136.03 0,128.30 0,104.40 0,062.77 
1/0,6 0,178.52 0,169.10 0,139.76 0,087.16 
1/0,8 0,205.70 0,194.88 0,161.11 0,099.18 
0,50 1,0 0,229.33 0,217.79 0,181.14 0,112.75 0 
0,8 0,252.48 , 0,240.36 0,201.39 0,126.95 
0,6 0,280.53 0,268.00 0,227.13 0,145.96 
0,4 0,314.51 0,302.44 0,261.27 0,173.60 
1/0,4 0,088.30 0,083.74 0,069.44 0,043.15 
1/0,6 “0,111.85 0,106.48 0,089.46 0,057.53 
1/0,8 0,126.81 0,120.85 0,101.64 0,064.96 
0,75 1,0 0,140.02 0,133.86 0,112.73 0,072.84 0 
0,8 0,152.34 0,145.66 0,123.98 0,081.08 
0,6 0,167.57 0,160.72 0,138.19 0,092.10 
0,4 0,185.98 0,179.44 0,156.95 0,107.89 
i _ 0 0 0 0 0 
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Tabelle 3. Werte des k,.-Koeffizienten 


y/b 
xa hilfe 
0,25 0,50 0,75 1 
0 es 0 0 0 0 

1/0,4 0,030.378 0,063.074 0,097.475 0,135.609 
1/0,6 0,032.298 0,066.940 0,106.044 0,151.310 
1/0,8 0,032.132 0,067.622 0,109.622 0,160.720 

0,25 1,0 0,031.250 0,066.779 0,110.767 0,166.783 
0,8 0,029.712 0,064.709 0,110.405 0,171.703 
0,6 0,026.854 0,060.334 0,107.763 0,176.414 
0,4 0,021.783 0,051.505 0,100.278 0,180.408 
1/0.4 0,058.457 0,119.871 0,186.602 0,261.817 
1/0,6 0,060.622 0,126.248 0,201.689 0,290.902 
1/0,8 0,060.036 0,126.959 0,207.691 0,308.305 

0,50 1,0 0,058.221 0,125.002 0,209.305 0,319.515 
0,8 0,055.230 0,120.814 0,208.123 0,328.607 
0,6 0,049.806 0,112.334 0,202.576 0.337.314 
0,4 0,040.319 0,095.610 0,187.843 0,344.693 
1/0,4 0,078.474 0,162.031 0,255.524 0,365.278 
1/0,6 0,080.578 0,168.931 0,273.576 0.403.373 
1/0,8 0,079.379 0,168.913 0,280.145 0,426.135 

0,75 1,0 0,076.734 0,165.699 0,281.251 0,440.789 
0,8 0,072.609 0,159.655 0,278.708 0,452.671 
0,6 0,065.323 0,147.972 0,270.217 0,464.047 
0,4 0,052.770 0,125.526 0,249.364 0,473.689 
1/0,4 0,086.01] 0,178.365 0,284.218 0,417.322 
1/0,6 0,087.888 0,184.933 0,302.334 0,458.599 
1/0,8 0,086.367 0,184.369 0,308.459 0,483.247 

1 1,0 0,083.369 0,180.532 0,308.928 0,499.111 
0,8 0,078.799 0,173.688 0,305.490 0,511.973 
0,6 0,070.768 0,160.736 0.295.493 0,524.287 
0,4 0,057.159 0,136.150 0,271.954 0,534.723 


(Hingegangen am 3. Oktober 1960) 
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Zum Problem der an den Rindern fest eingespannten, 
elastischen Platte* 


Von W. Gerisch, Berlin ** 
Mit 15 Textabbildungen 


Zusammenfassung. Die vorliegende Arbeit soll vorwiegend der Zusammenstellung numerischer 
Resultate iiber eingespannte, freischwingende, quadratische, isotrope, elastische, schubbeanspruchte 
Platten dienen. Beziiglich der Methode wird ein bereits in zwei friiheren Mitteilungen eingeschlage- 
ner Weg weiter verfolgt. Unter anderem geben wir mit Hilfe eines fiir die numerischen Rechnungen 
weniger geeigneten Orthonormalsystems einen Uberblick iiber den Verlauf des ersten Eigenwertes 
als Funktion eines Schubparameters sowie iiber die zugehérige Schar von Eigenfunktionen. Speziell 
sind in diesen Resultaten Aussagen zum Thema Schubbeulung enthalten. Unter Benutzung eines 
bereits friher untersuchten Orthonormalsystems von zugeordneten Kugelfunktionen, auf dessen 
vielseitige Anwendbarkeit bei derartigen Problemen schon hingewiesen wurde, werden die hinsicht- 
lich des Ausbeulens erhaltenen Resultate verfeinert. Dabei gelingt es insbesondere infolge der auBer- 
ordentlich einfachen und itbersichtlichen, zum Einsatz elektronischer Rechenmaschinen geeig- 
neten Verhaltnisse zu zeigen, daB gewisse Angaben anderer Autoren hinsichtlich einer unteren 
Schranke nicht stichhaltig sind. Umfangreiche Vergleiche der erhaltenen Resultate mit denen ande- 
rer Autoren, auch fir den Fall drehbar gelagerter Platten, sollen zur Abrundung des gewonnenen 
Uberblicks dienen. 


I. Einleitung 


In der vorliegenden Arbeit wollen wir hauptsachlich einige numerische Resultate 
zusammenstellen, die sich bei der Untersuchung eingespannter, freischwingender, qua- 
dratischer, isotroper, elastischer, schubbeanspruchter Platten ergeben haben. Gleich- 
zeitig sollen friihere Betrachtungen zum selben Gegenstand!? fortgesetzt und erweitert 
werden. Was die Behandlung des betreffenden Eigenwertproblems anbetrifft, so ist 
zu bemerken, daB die zu gebenden Resultate mit einer bereits in 1 erdrterten Methode 
gefunden wurden. Daselbst und in ? haben wir gezeigt, wie man zu dem hier mageb- 
lichen Differentialausdruck zunachst einen Differentialoperator A in einem passenden 
Hilbertschen Raum L? gewinnen kann. Weiterhin haben wir in ! unter anderem tiber 
ein Orthonormalsystem berichtet, mit dessen Hilfe dem Operator A in L? eine Kern- 
matrix (@mns:) in einem Hilbertschen Folgenraum /? zugeordnet wurde. Durch besonders 
gliickliche Wahl des erwahnten Orthonormalsystems lieB sich ein geschlossener Ausdruck 
fiir Gn>, angeben, so daB — wie an Hand eines Beispiels gezeigt — zu den jetzt in 
Form von unendlichen Gleichungssystemen vorliegenden Problemstellungen durch 
abschnittweises Vorgehen auf besonders einfachem Weg Naherungslosungen gewonnen 
werden kénnen. Bevor die Konstruktion dieser geschlossen darstellbaren Kernmatrix 
gelang, war mit einem anderen Orthonormalsystem versucht worden, die angedeuteten 
Gedankengange zu realisieren. Uber diese Bemiihungen, die zu den eingangs erwahn- 
ten Resultaten fiihrten, soll méglichst kurz berichtet werden. Seinerzeit ist es namlich 
nicht gelungen, eine geschlossen darstellbare Matrix zu finden, und die Naherungs- 
rechnungen fielen demzufolge verhaltnismafig kompliziert aus. Unsere diesbeziiglichen 
Ausfiihrungen werden wir abschliefend durch einige Bemerkungen tiber das Aus- 
beulen quadratischer Platten erginzen. Bei dieser Gelegenheit sollen auch weitere 
Einzelheiten und Ergebnisse der Rechnung mit dem oben erwahnten, besonders giin- 
stigen Orthonormalsystem mitgeteilt werden. 


* Durch den Senator fiir Wirtschaft und Kredit wurden Mittel aus dem ERP-Sondervermogen 
zur Verfigung gestellt. Der Verfasser dankt den Herren W. Boéhmert und H. Maus fir Mithilfe 
bei den numerischen Rechnungen und fir die Anfertigung der Figuren. 


** Bundesanstalt fir Materialprifung. 
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Il. Problemstellung und Bezeichnungen 


Freie Transversalschwingungen quadratischer, isotroper, elastischer, schubbean- 
spruchter, eingespannter Platten lassen sich durch die Differentialgleichung 


2 2 @ Tie GANG 
saw +3-(2)- 25 —G(sfiorw=o c 


mit den Randbedingungen 


é . ¢ 
w= 0, 5. = 0 fiir x = 0 baw. 2; w= 0, 5, = 0 fir y= 0 baw. x (2) 


beschreiben. Bezeichnungen: a, h Plattenabmessungen, w Durchbiegung als Funktion 
von (a, y), D Plattensteifigkeit, 9 Dichte, 7 Schubkraft am Rande pro Langeneinheit, 
o Kreisfrequenz. Mit Hilfe gewisser dimensionsloser GréBen 


sa (Sf an thlee ® 


14Bt sich die Gleichung (1) wie folgt iibersichtlicher schreiben: 


AAw+o-25 — lw=0. (4) 


Wegen der Konstruktion eines Differentialoperators A aus (2) und (4) in einem Hilbert- 


schen Raum JL? (Skalarprodukt von der Form i [ f (a, y) g (x, y) dw dy) verweisen 
0 0 

wir auf 2. Daselbst sowie in ! wurde allerdings mit A ein formal modifizierter Differen- 
41 41 

tialoperator in L? mit { | f (x, y) g (x, y) dw dy bezeichnet. Da dieser Operator in 
Sah Sa 

der vorliegenden Arbeit lediglich unter Punkt VI eine Rolle spielt, behalten wir uns 

die Bezeichnung A fiir den zu (2) und (4) gehorigen Operator vor und kennzeichnen den 

modifizierten Operator sowie dessen Parameter und Kernmatrix durch einen Strich. 


Den Ubergang von L? nach /? (Skalarprodukt von der Form 3’) Ann Bnn) voll- 


m n 


ziehen wir wiederum mit Hilfe der Schrodingerschen Formel unter Benutzung eines fiir 
A zulassigen, vollstandigen Orthonormalsystems {@,,, (x, y)}. Die Forderung nach 
Zulassigkeit eines Orthonormalsystems ist dabei bekanntlich nicht allein mit der Erfiil- 
lung der Randbedingungen (2) befriedigt. Wir werden auf die sich durch diesen Sach- 
verhalt ergebende Problematik an anderer Stelle noch zu sprechen kommen. Fiir die 
Elemente der Kernmatrix (@nn5:) in 2? gilt schlieBlich: 


ot 
— 


Gmnst = fifa Dyn Dg: dex dy. 
O.10) 


Auf eine auch fiir die numerische Behandlung besonders wichtige Eigenschaft des 
Operators A bzw. A’ wurde ebenfalls schon frither hingewiesen*. Es existieren namlich 
im Falle # + 0 Paare von Unterraéumen in L? und /?, die den betrachteten Operator 
und seine Kernmatrizen vollstandig reduzieren**. Die Untersuchung la&t sich damit 
in eine Untersuchung der Einschrankungen auf die einzelnen Unterraiume in den 
betreffenden Raumen iiberfiithren. Im Gegensatz zu } beschrinken wir uns in der vor- 


ES 2OSe 
** Fir den separat ausfithrlicher zu untersuchenden Fall # = 0 (Freischwingende, schubfreie 
Platte) erfolgt die vollstandige Reduktion sogar durch 6-tupel von Unterraumen. 
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lhegenden Arbeit unter Ausnutzung des beschriebenen Sachverhaltes insbesondere auch 
im Hinblick auf die Anwendungen auf eine Untersuchung desjenigen Teiles von A bzw. 
(mnsi), der den ersten Higenwert enthalt. Bei der Berechnung von Naherungen fiir © 
EHigenwerte und Kigenelemente des so definierten Teiles von A bzw. (Qnnst) gehen wir 
wie in ' vor. Wir betrachten zunachst einen gewissen Abschnitt dieses Teiles der Kern- 
matrix und erhalten Naherungen fiir die EKigenwerte in Gestalt der Abschnittseigen- 
werte. Naiherungen fiir die Eigenelemente sind die zugehérigen Abschnittseigenelemente 
in /? und die diesen Elementen entsprechenden Elemente in L?. Es wird sich unter ande- 
rem zeigen, da durch dieses Vorgehen mit ertraglichem Rechenaufwand ein zufrieden- 
stellender Uberblick iiber den ersten Eigenwert als Funktion von ? und die zugehorige 
Schar von Kigenfunktionen unserers Ausgangsproblems (2), (4) zu gewinnen ist. 


Ill. Orthonormalsystem und Kernmatrix 


Nach! * ist das Funktionensystem {¥,,,, (x, y)} mit den Elementen 
mn (XL, Y) = Wm(L)- nly) = (wx — x) sin ma: (ny — y*)-sinn y(m,n=1,2,3,...) (6) 


nau 


vollstandig im Raume L? mit { lf (x, y) g(x, y) dxdy. {Wrnn (x, y)} soll in ein voll- 
0 0 

‘stindiges Orthonormalsystem {@,,, (x, y)} iiberfiihrt werden. Wegen des symme- 

trischen Aufbaues von (6) erhalten die Elemente von {@,,, (x, y)} die Gestalt 


Pin (2%, Y) = Pm (%)*Gn(y) (m,n = I, 2,3, ...) (7) 

und die Aufgabe reduziert sich zunichst dahin, {y, (w)} zu {q (u)} in L? mit | } (uw) g (uw) du 
0 

zu orthonormalisieren. Der zuletzt genannte Raum ist in zwei Unterraume zer- 


legbar. Ein Unterraum wird durch das Teilsystem {y, (u), k = 1.3, 5,...} aufge- 
‘spannt, der andere durch {y, (w), & = 2, 4, 6,...}. Demzufolge kann man diese beiden 
Teilsysteme fiir sich orthonormalisieren. Unter Benutzung der Relation 


9 
a 


| (wu — w)2- sin kusinl udu = 


16 16 


_ kel 
| er Ve ie 
12 (& — 1) =0 (mod 2) 
| eG aK. ae 

0 (k —1)'=1 (mod 2) 


liefert das Orthonormalisierungsverfahren ** fiir die Elemente der beiden orthonormier- 
ten Teilsysteme folgendes Ergebnis: 


b. Ak (u) 
Pk (6) = feet (eu — U?) ee (9) 


mit 


* §. 232 und 228. 
** S. 193 bis 196. 
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4 
le a Abel ir ee 
(+3) (5-3) SE Qe a te 
py i ask a te 
(50 — 7) (sch gs) a> 7k OE 


A, = i, 2/1. a ee 
(5 5) (a-i) (e+e F 


(q = Hy (ek : my)” 


sin wu 


sp 7s ae | 
(F | eee 


sin wu sin 3 wu 


( 1 (5 ss x) 
Aways (55 os i) (5 a 3) - i) 
8) GH) & 


(aa Dee Gea) 


| sin. we sind sind... 


fiir die Elemente des ersten Teilsystems und 


Pee at Ck ek i 
an + 45) Bt dee See 


% (Gan WeE=Dy) 


sin k wu 


Age 


Ha ASR Bera yy ieee 


eee e ee Eo 
ge 98) N58 z) (ae +H sas 


1 1 ; 
Foes: (Fh ay) 


— 


1 1 : 
(q (K+ 2))* (4(k — 2) ¥) 


(10) 
bel 
bs 
k= 6,7, 9e0e 
ae 
2 
eee ME a 
(11) 
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{ sin 2 wu k= 2 
ietear) (si) k=4 
sn2u sin4u ei 
(a+ a) (ga) (em) (@erap- eee) 
Se ea | 
(8) (8) (+) | 
b= 6,8,10)... 
(q a (3 (k ss ay) 
sin2u sin4du sin6u ... sin ku 


fiir die Elemente des zweiten Teilsystems. Besonders aus den beiden letzten Glei- 
chungen ist ersichtlich, daB die hier vorliegenden GesetzmaBigkeiten, verglichen mit 
den Verhaltnissen bei den in 1 untersuchten Orthonormalsystemen, wesentlich kompli- 
zierter und unhandlicher sind. Bis zum AbschluB dieser Arbeit ist es deshalb auch nicht 
gelungen, eine geschlossene Darstellung fiir g, (w) zu finden. Ebenfalls ist nicht geklart, 
welche Differentialgleichung g, (w) erfiillt und welche Beziehungen zu anderen bekann- 
ten Funktionensystemen bestehen. (7) und (9) liefern unter Benutzung von (10) und 
(11) formal den gesuchten Ausdruck fiir @,,,, (z, y). Zur Erginzung unserer Ausfiihrun- 
gen geben wir nun zahlenmafige Ausdriicke fiir die ersten Glieder des Orthonormal- 
systems {q (w)} sowie fiir fiinf Glieder des Orthonormalsystems {@,,,, (a, y)} an 


Q, (u) = (nu — uw?) - 0, 366 212: sin u 
: = (1 u — u2) - 0, 436 537 - sin 2 u (12) 


oS 
i 
2 


(x y — y”)- 0, 1384 11 sin x sin y 

( y — y?)- (0, 051 28 sin x sin y + 0, 173 12 sin 2 sin 3 y) 
( 

( 


) 
y) = (xu — 2?) (x y — y?) - (0, 051 28 sin x sin y + 0, 173 12 sin z sin y) 
oo (%; y) = (mw — 2?) (wy — y”)- 0, 190 54 sin 2 xsin2 y (13) 
33 (X, Y) = (a & — a) (wy — y?): (0, 019 61 sin « sin y + 0, 066 19 (sin w sin 3 y + 


+ sin 3 asin y) + 0, 223 46 sin 3 2 sin 3 y) 


AuBerdem stellen wir in den Tabellen 1 und 2 Zahlenmaterial tiber die genannten 
Orthonormalsysteme zur Anwendung unter Punkt V bereit. Abb. 1 zeigt den Verlauf 
der ersten drei Funktionen des Orthonormalsystems {q, (u)}. Wir wollen abschlieBend 
unsere unter Punkt II gemachte Bemerkung iiber die Zerlegbarkeit von L? mit 


ij f f (x, y) g (x, y) dx dy in Unterraume ergiinzen. Die Zerlegung kann in vier und 
eo 0 
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somit auch auf verschiedene Weisen in drei und zwei Unterraume erfolgen. Die vier 
Unterriiume werden durch die Teilsysteme {Ynn (x, y) m, = 1, 3,5, .. as eee Chee, 
m =1,3,5,.. 2,0 = 2,4, 6,.5.}5 {Yon Os Wit enon 13,5, ne} ) ee ma ee 
m,n = 2,4,6,...} des Funktionensystems {Poin (a, y)} oder durch die entsprechenden 
Teilsysteme des Orthonormalsystems {®,n (x, y)} aufgespannt. 


Tabelle 1 
a a ee eee 
U Uh es P, (u) Po (U) Pz (U) 
. ie 2 1 2 1 | 2 | 1 2 
2 EEE SR 
0 I 0 0 0 
tee TM ae 0,071454 +] — 0,16456 0,27934 
OH eas NE Se 0,26999 a 0,51823 0,74397 
Ml ae We te +) + 0.47921 +} — 0,80783 +} + 0,8018] 
Lis Gas Wes ir 0,69558 eee (ae 0,82916 4. |e 0,26597 
Ws BG Ces ts a) 0,84853 +] =— 0,80215 = — 0,15004 
Wh, Gs he ae ae | aie 0,90359 0 =i 0,26280 
plu) | 


7 


als 


Abb. 1. Verlauf der ersten drei Funktionen des Orthonormalsystems {q, (u)} nach (12) 
Das beschriebene vollstandige Orthonormalsystem {@,,, («, y)} soll nun unter 


Voraussetzung der Zulassigkeit zur Darstellung des Operators A herangezogen werden. 
Mit (7) folgt aus (5): 


Oe nse Sey Qsdx + Ont + 2 - [ey Qs dx - |g gidy (0x: Kroneckersches Symbol) 
0 0 0 


9 


+ dns {OY ge dy +0 J pn g.da- |g, we dy. (14) 
49 0 0 


Fiir die in (14) auftretenden Integrale findet man unter Benutzung von (9) bis (11): 


{oigqidu = 0 (k —1)=0(mod2) (15) 
0 


[vk ordu = — |v du (16) 
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Tabelle 2 
Se 
7] 7] & v ®,, (2%, y) | P15 (”, y) | ®5, (2, y) | Dy. (x, Y) ®,, (x, Y) 

t a : ee EPPS eee Liesls Baer elcid las | 
eee 
0 Tt 0 0 0 0 

Oo | x 0 0 0 | 0 
Way 76 |!/,92¢+[+|+|+/0,0081/+|+|+|+|0,02001+|+|+1+|0,02001+|—|—|4 0,0780 
1/_7| 5/40 0,0179 0,0701|+|+|+|+10,05321+/—|—|+ 0,2079 
1/97 | */y970) toc] 2/q2c1+|+|+1+|0,034214 |+ 0,1339}+|+|+/+/0,0573|+|—|—|+ 0,2240 
Weer tela +/0,0497)/+ |+]4 0,1943}+ |+|+)/+/0,0190]+|—|—|+ 0,0743 
5/490 | 7/4970 |+|+|+|+]0,0606)+ |+)+/+/0,2370|—|—|—|—|0,0341}+]—|—|+ 0,1333 
1/y70| 3/570 | +|0,06441+|+|+/+/0,2524|-|—|—|-|0,0587 0,2293 
0 Ms 0 0 0 0 
Tf 9 |22/45 | +|+|+|+]0,0179]+|+|+/+|0,0532}+|+/+/+/0,0701}+|—|—|+ 0,2079 
1/.m| 5/_7¢|+|4 0,0630'+|+|+|/+/0,18671+)+|+|+/0,1867]+|—|—|+ 0,5535 
1/50 | 5/gx| 1/40} 8/470|+|+/+/+ 0,1203} 0,3565]+ + /0,20131+|—|-|+ 0,5965 
"/am| */am|+|+/+|+|0,1746 +/0,5175]+|+|+|+/0,0668]+|—|—|+ 0,1979 
5/159 | 7/4990 |+|-+|+|+|0,2130 +|0,6313|—-|—|—|—|0,1198}+|—|—|+ 0,3552 
1/.a| 1/500 |+|+|-+|+|0,2268]+|+|+/+|0,6723|—|—|—|--|0,2060 0,6107 
0 I 0) 0) 0 0 
1/457 |4/,,+/+|+|+/0,034214|+/+/+10,057314|+|+|4(0,1339}4|—|—|+ 0,2240 
1/,0| ®/go0|+|+|-+|+ 0,1203/ 0,20131+|4 0,3565]+|—|—|+ 0,5965 
1/42) */4| "/ar| ®/am|+|+| +) +|0,2296 +10,3842}+ 0,3842)+)—|—|+ 0,6429 
1/7 | 2/576(4 +10,3333/+|+|+/+|0,5577/+/+|+|+|0,1275}+|—|—|+ 0,2133 
5/9 | 7/1970 |+\+|+|+/0,4066 +|/+/+|+/0,6804I-|—|—|—/0,22881+|—|—|+ 0,3828 
1/,a0| 2/,%|+|+|+|+|0,4330}+ 0,7245)—|—|—|-|0,3934) | |_| 0,6582 
e 0 I 0 0 0 
I ee Be PE 0,0190}+|+/ +/+ nies ire 0,0743 
ve 0,06681+|+|+|4 +|—|—|+]0,4296}+)+]+/+|0,1979 
U/y0| 2/gac| *yoe| 2/qoe]+|+|-4 14+ 0,12751+|+|4+/+ -|-|-|+]0,6697}+/+|+|+|0,2133 
1/40 | 2/32 |+/4 0,1850)1+|+|+/+ +|—|—|+/0,6874+]+|/+|+|/0,0707 
5/9 | "ap [tlt +|+ 0,2257|—|—|—|— .|-|—|+]0,4341|-|—|-]—|0,1270 
1/30 | 157 +/+]+]; 0,2403/—|—|—|— S O f+ |-|-|-|0,2395 
0 0 
0 I 0 
1/4970 |1/,520]+|+|-+|+ 0,03411+|+/+/+|0,23701+|—|—|+/0,0862|—|—|—|—-]0,1333 
: _ 0,2713)—|—|—|—|0,38552 
eee Sliaaal alae ella 0,1198}+]+ > > 
5/57 7/497 ae ae ses fene |b 0,2288)+|-4 +|—]—|+/0,4229)—|—/—|—|0,3828 
Peat om | BE nb relict 0,332If¢/4/+/4 .|—|_|+]0,4341!-|-|-|-]0,1270 
5/9 | 7/4976 1+ |+|+|+ 0,4051|—|—|—|—|0,4051}+|—|—|+/0,2741}+ |+|+)+|0,2279 
Seat 11500 a ae eee eS ee 0 pata eas 
= 0 
0 1% 0 0 ! 
1/4570 [2/400] +|-+|+|/+ 0,0587}+-|+|+|+ 0 |-|-|-|-|0,2293 
Lege |) 4c t+ |+|+ 0,2060}+ |+/+|+ ; =|—|—\— ae oe 
ae a Dea Nie fermen) 
1/7 1/,a0| 1/420 8/47 1+|+]+|+ res banal a wae ae 
A mG | 2 [gat |+ |e 0,5710)4-|-+]+)+ 
5/97 | Wyatt! +|+|+ 0,6966}—|—|!—|- OQ jes Be 
1/4 m| 2/gu|+|+|+|+ OAT a= |= Opti +}t]4 [9 
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uw 


/ ‘ du 4 . U 
bene 3% VAg—1-A 


ke 


= vi fw — 2): ew — w)- Fil) Aiu) du + 
[AL aL ae 


uw 


+ [(ou— usp A; (au) - Ar (uw) d 


s 


(k —1) = 1 (mod 2) (17) 


{oe gidu=0 - (¢ = 2,4) (k-- 1) = 1 (mod 2) (18) 
0 

[of gy du = { oe \) du (i = 2, 4) (19) 
0 0 

fos? gidu == ; -|(— 4) ei eo (u) - Ar (w) du 
eae 3% VAp—1- Ag: Ai-1° Al 3 


+i-[(a— 2 wu) (2 wu — w2) » Ay’-D (w) + Ay (w) du 


a {a u— u?)- Ay) (u) - Ai (w) du 
10) 


(i = 2,4) (k —1) = 0 (mod 2) (20) 


Versuche, (17) und (20) zu vereinfachen, scheitern, solange fiir die auch hier auftreten- 
den Determinanten keine befriedigenden Ausrechnungen gefunden sind. Dies ist auch 
der Grund, weshaib fiir a,,,;, kein brauchbarer Ausdruck gegeben werden kann. Fiir 
die numerische Auswertung von (17) und (20) sind neben (8) die folgenden Relationen 
von Wichtigkeit: 


uw 
f@—- 2u): (xu — u?)-sinku snludu = 
10) 


L 1 1 1 (21) 


[ (au — u)+ cos ku sini udu = 
i) 


0 (k —1) =0 eM 

= Faia 1 ) fe Bi 1 22) 
—2a8- (Gop — waa) + 2-H — Baap) &—)=1 (mod 2)| 

[ru — v2) sin kw sind udu = 
| 0 (k —-l)=1 (mod 2) 
1 1 . 

Sa eE la alee = ee (23) 

| | (k — 1) =0 (mod 2) 
os 42 12 | =f 
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(x — 2u) (mu — u?)- cosku-sinludu= 


ot 
va 


0 (k — 1) =1 (mod 2) 
ik 1 
Os 6x: |p = ca steal, (24) 
(k — 1) =0 (mod 2) 


Mit (14) bis (24) und (8) sind wir nun in der Lage, fiir die Naherungsrechnungen einen 
Abschnitt des zu dem unter Punkt II definierten Teiles von A gehérigen Teiles der Kern- 
matrix (Qmns:) anzugeben. Die dort fiir #+ 0 erwahnten beiden Unterraume werden 
in L* durch die Teilorthonormalsysteme {@,,, (x, y), (m — n) = 0 (mod 2)} und 
{Brn (x, y), (m — n) = 1 (mod 2)} aufgespannt. Fiir den fiinften Abschnitt ,(a),,) 
((m — n) = 0, (s — t) = 0 (mod 2)) ergibt sich schlieBlich zahlenmaBig der folgende 
Ausdruck : 


13,55107  — 4,058355 — 4,058355 1,14977 - 9 2,07205 

—4,058355  185,00641 2,07211 —1,93394-9 —22,703369 

i(a®),,) =| —4,058355 2,07211 185,00641  —-1,93394- 9 —22,703369 
1,14977- 9 — 1,93394-9 — 1,93394-& 125,98852 3,25311-9 

2.07205  —22,703369 —22,703369  3,25311-9 524,14766 
(25) 


IV. Eigenwerte 


Mit Hilfe von (25) machen wir nun N&herungsaussagen tiber gewisse EKigenwerte 
unseres Problems (2), (4). Die zu (25) gehorige fiinfreihige Sakulardeterminante 1aBt 
sich aus Symmetriegriinden sofort in eine vierreihige Sakulardeterminante tiberfiihren. 
Diese liefert den folgenden Ausdruck fiir 3 (A): 


Aa — J4 + 43+ 850,765 77 — A2- 197 940,138 + A- 14731 216,74 — 163 470 014,00 (26) 
lies at — j?- 19,384 943 + 72-6 462,614 188 — 180 290,206 956 = 


Tabelle 3 gibt zunachst eine Zusammenstellung der Null- und Polstellen des Aus- 
drucks (26). Offenbar sind die Nullstellen von (26) Naherungen fiir entsprechende 
Eigenwerte von (2), (4) bei # = 0*. Da wir fiir die numerischen Rechnungen lediglich 
den in (25) gegebenen fiinfreihigen Abschnitt des einen Teiles der Kernmatrix zugrunde 
legten, ist nach dem Aufbau dieses Abschnitts insbesondere fiir den ersten und dritten 
Abschnittseigenwert eine gute Ubereinstimmung mit bereits berechneten Naherungs- 
werten zu erwarten. Wir haben deswegen in Tabelle 3 nur diese Werte mit den zugehori- 


Tabelle 3 
2° (0) = 13,361 | A (0) = 125,985 | 1 (0) = 185,268 | 4 (0) = 524,148 
cee Va (0) = 3,655 12 (0) = 13,611 
RP ety nee abi yrai it Gatti!) | si uig debi ole pe ake et 
S.Iguchi ....... 3,646 13,393 | 
i ay unaani. 2h, J (cc) = 30,73 | AW (co) = 302,65 


8© (0) = 30,1116 


* Freischwingende, schubfreie Platte. 
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gen von S. Iguchiin 4 ermittelten Werten verglichen. Die aufgefiihrte Naherung #°? (0) 
fiir den Beulwert diskutieren wir unter Punkt V. Zur Erginzung von Tabelle 3 werden 
in Tabelle 4 einige Zwischenwerte von # (A) zu- 
sammengestellt. Abb. 2 gibt eine Ubersicht tiber 
den Verlauf der wichtigeren zu (26) gehorigen 
Umkehrfunktion 4 (9). Auch hier messen wir den 
beiden im Intervall 0 < # < 9 (0) ausgezogenen 
Kurvenasten 4{° (#) und A? (8) auf Grund der 
benutzten Naherung besonderes Gewicht bei. 
Bemerkenswert ist, daB AY) (#) im Gegensatz zu 
A) (9) mit wachsendem # zunimmt. Naherungs- 
rechnungen zum ersten Kigenwert des gleichen 
Problems wurden bereits von M. Hasegawa in ° 


oF (ay 


0 5 10 15 20 25 30 35 


Abb. 2. Niaherungsweise bestimmte Abb. 3. Naherungsweise bestimmter erster Eigenwert 
Eigenwerte als Funktion von # nach als Funktion von #, Vergleich mit den Angaben von 
(25) und (26) M. Hasegawa und K. Munekata nach® 


mit Hilfe des Rayleighschen Prinzips angestellt. Daselbst wird zusatzlich iiber das 
Ergebnis weiterer Naherungsrechnungen von K. Munekata* berichtet. Die beiden 
mageblichen Kurven haben wir nach passender Umrechnung in Abb. 3 unserem 
Kurvenast /{°) (9) gegeniibergestellt. Der Vergleich fallt zufriedenstellend aus. 


Tabelle 4 
A | 3 | 6 | 9 | 12 | 100 | 190 | 210 | 550 
(A) 27,411 | 24,016 | 19,349 | 11,415 | 18,887 | 7,958 | 20,481 | 27,612 


V. Eigenelemente 


Wir wollen uns in diesem Abschnitt zunachst einen Uberblick iiber die zu dem unter 
Punkt IV naherungsweise bestimmten ersten Eigenwert unseres Problems (2), (4) 
gehorige Schar von Kigenfunktionen verschaffen. Zu diesem Zwecke stellen wir in Tab. 5 
sechs zu (25) gehérige auf 1 normierte Abschnittseigenelemente in /? zusammen. 0 
nimmt hierbei genau die sechs Werte aus dem Intervall 0 < 9 < 8® (0) an, zu denen 
in den Tabellen 3 und 4 4‘ (9) zu finden ist. Die sechs Punkte sind zusatzlich in den 


* Vortrag auf der Jahrestagung der Japanischen Physikalischen Gesellschaft 1948. 
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Abb. 2 und 3 markiert. Den in Tabelle 5 zusammengestellten Abschnittseigenelementen 
entsprechen in L? fiir die jeweiligen Werte von % die Elemente 
wi” (x, y5 9) == AT (P) > D, (x, y) + Af'is (By (Dy, (, y) + Ds, (x, y)) a 
+ Ales (9) - Bop (x, y) + AM3s (8) - Day (zx, y) (27) 
Die Funktionen (27) sind, 


wie unter Punkt II be- Tabelle 5 
merkt, die gesuchten aA 
Naherungen. Eine Tabel- od AY (9) 3: oP a Ase (9) As (8) 
lierung dieser Funktionen eae 
a Benutzung von 0 + 0,999 + 0,023 0 — 0,002 
abelle 2 findet sich fiir 
5 3 11,415 0,994 0,009 —0,112 +0,005 
die maSgeblichen Werte . a: 
von % in Tabelle 6. Aus 19,349 + 0,977 — 0,019 — 0,209 + 0,020 
Symmetriegriinden konn- 24,016 + 0,958 — 0,046 — 0,278 + 0,034 
ten wir uns hierbei auf 27,411 +0,936 | —0,070 | —0,334 | +0,047 
ee 30,112 0,912 0,092 0,382 0,060 
0 < y < 2/2) beschran- y as om ies at 


ken. In den Abb. 4 bis 9 

wird durch graphische Auswertung der rechnerischen Ergebnisse aus Tabelle 6 der ange- 
strebte Uberblick iiber die betreffende Schar von Eigenfunktionen gegeben. Bevor wir 
das erhaltene Resultat physikalisch interpretieren, wollen wir noch eine Bemerkung 
tiber die in Abb. 4 in Grund- und Aufri& wiedergegebene aus zwei Kurvenasten beste- 
hende Raumkurve einfiigen. Diese Raumkurve entsteht, wenn man die in den Abb. 4 
bis 9 dargestellten naherungsweise bestimmten Eigenfunktionen der betrachteten 
Schar zum Schnitt bringt. Auf Grund unserer Naherungsrechnungen steht demnach 
zu vermuten, daB allen Kigenfunktionen der Schar eine solche Raumkurve gemeinsam 
ist. Die diesbeziiglichen Verhaltnisse in den zu héheren Eigenwerten gehorigen Scharen 
von Eigenfunktionen sollen unter Umstanden Gegenstand einer besonderen Unter- 
suchung werden. Physikalisch 1iBt sich auf Grund unserer Abb. 3 bis 9 hinsichtlich 
des Grundtones der in der vorliegenden Arbeit untersuchten freischwingenden, schub- 
beanspruchten Platte folgendes aussagen. Der Grundton wird mit Zunahme der Schub- 
krafte an den Randern tiefer. Die Platte selbst deformiert sich in anschaulicher Weise. 
Wenn die Schubkrafte so groB werden, daB # = 24 wird, treten Knotenlinien auf 
(Chladnische Klangfiguren). Mit weiterhin zunehmenden Schubkraften verschieben 
sich die Knotenlinien in Richtung Plattenmitte. Die sich etwa bei stroboskopischer 
Beleuchtung abzeichnende Form des Platteninneren zwischen Knotenlinien und 
Berandung 1aBt sich treffend als Beule charakterisieren. Bei derjenigen GrédBe der 
Schubkrafte, bei der sich 3 ~ # (0) = 30 ergibt, wird der Grundton der Platte zu Null. 
Die Platte beult unter diesen Verhaltnissen mindestens bei der geringsten Storung aus. 
Die betreffenden Schubkrafte nennen wir kritisch, den zugehorigen Wert %,;, von # 
Beulwert und die Darstellung der Eigenfunktion nach Abb. 9 Beulfigur. Die hier be- 
gonnene Diskussion des Beulvorganges fithren wir unter Punkt VI weiter. Daselbst 
werden wir unter anderem auch naherungsweise bestimmte Beulfiguren anderer Autoren 
zum Vergleich mit unserer Abb. 9 zusammenstellen. 

Zur Erginzung unserer Ausfiihrungen iiber die zu dem ersten Kigenwert von (2), 
(4) gehorige Schar von Kigenfunktionen geben wir abschliefbend noch eine Ubersicht 
iiber die zu dem durch 4° (#) naherungsweise bestimmten Eigenwert gehérige Schar 
von Eigenfunktionen. Tabelle 7 enthalt wiederum Abschnittseigenelemente in EB. 
Wegen der Werte von # wird auf die Tabellen 3 und 4 sowie auf Abb. 2 verwiesen. Die 
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Tabelle 6 


wy (x, Y) 
# = 24,016 


0 


0 0 


0 0 


0 0 


4 


fer) 0 


0 0 


0 0 


0 0 


0 


5/4 a 
3/4% 
7/50 
"yom 
*/y 7 


"12% 


11/, 71] + 0,006|+0,003 


+ 0,020)+0,010 
+ 0,038]+ 0,022 


+0,009} ~0 


+0,011]—0,002 


+0,013}]—0,003 


+0,030}]+ 0,002 
+ 0,052|+ 0,007 


+ 0,037}—0,005 
+ 0,062}—0,005 


0,054]+ 0,036 
+ 0,066}]+ 0,052 
+0,069]+0,065 


+0,067|+ 0,018 
+0,071]+ 0,035 
+0,065]+ 0,055 


+0,074]+ 0,003 
+0,071)+0,020 


0,042]—0,010 
+0,069]—0,015 
+0,078]—0,010 
+ 0,068]+ 0,007 


+0,055}+ 0,045 


+0,045]+ 0,036 


+0,015}—0,005 
+0,047]—0,015 
+ 0,074|—0,024 
+0,081}—0,022, 
+0,065]—0,004 
+ 0,036}+ 0,027 


+0,016 
+0,050 
+0,078 
+0,082 
+0,062 
+0,027 


% 


) 0 


0 0 


0 0 


0 0 


0 0 


to) 


9/4 % 
3/40 
2/47 
“/1a% 
*/y% 


Say 
/6% 


11/,.77 + 0,020] + 0,010 


+ 0,070]+ 0,039 
+0,132]+0,081 
+0,188]+ 0,132 
+ 0,224]+ 0,184 
+0,239|+ 0,227 


+0,030|+0,002 


+0,037}—0,005 


+0,042])—0,010 
+0,137}/—0,031 


+0,099}]+ 0,009 
+0,175]+ 0,031 
+0,228}+ 0,074 
+ 0,245}+0,135 
+0,227}+0,201 


+0,122/—0,012 
+0,206]—0,006 
+0,253}+ 0,028 
+0,248]+ 0,093 
+0,201}+0,175 


+ 0,226]—0,038 


+0,047]—0,015 
+0,149]—0,046 
+0,242]—0,066 
+ 0,276}]—0,047 


+ 0,267|]—0,012 
+0,244/+ 0,056 
+0,175)+0,151 


+0,238]+ 0,022 
+0,151}+ 0,128 


+0,050 
+0,159 
+0,254 
+0,282 
+0,230 
+0,128 


LA 


0 0 


0 0 


0 0 


0 0 


0 0 


0 


5/4% 
3/4% 
*/3.% 
"19% 
Mom 


"4% 


11/,,7|+0,038}+ 0,022 


+ 0,052]+ 0,007 


+0,132|+0,081 
+0,246]+ 0,165 
+0,349]+ 0,264 
+0,418}]+ 0,359 
+0,442/+ 0,430 


+0,175}+0,031 
+0,3121+0,086 
+0,4141+0,177 
0,454}+ 0,293 


+0,062]—0,005 
+0,206]—0,006 
+0,359]+ 0,026 
+0,457/+ 0,109 
+ 0,469}+0,238 


+ 0,430]+ 0,404 


+0,404)+ 0,377 


+0,069]—0,015 
+0,226]—0,038 
+0,388]—0,026 
+0,481]+0,050 
+0,473)+0,189 
+0,377/+0,351 


+0,074|—0,024 
+0,242|-—0,066 
+0,410|—0,072 
+0,498]—0,002 
+0,472|+0,145 
+0,3511+ 0,325 


+0,078 
+0,254 
+0,427 
+0,510 
+0,468 
+0,325 


1 


0 0 


0 0 


0 0 


0 0 


0 0 


0 


5/4 a 
3/40 
3/57 
"em 
"am 


117, 71+ 0,054]+ 0,036 


+0,188]+ 0,132 
+0,349]+ 0,264 
+0,492]+ 0,407 
+0,588}+ 0,536 
+0,621}+ 0,620 


+0,067|+0,018 
+0,228]+ 0,074 
+0,414)/+0,177 
+ 0,562]+ 0,324 
+0,634]+0,486 
+0,620]+ 0,616 


+ 0,074)+0,003 
+0,253]+0,028 
+0,457/+0,109 
+0,611)+0,258 
+0,667]+0,445 
+0,616}+0,609 


+ 0,078}—0,010 
+0,267]—0,012 
+0,481}+0,050 
+0,640]+ 0,201 
+0,686]+ 0,409 
+0,609}+ 0,600 


+0,081}—0,022 
+0,276]—0,047 
+0,498]—0,002 
+0,660]+0,149 
+ 0,699}+ 0,375 
+0,600]+ 0,590 


% 


0 0 


0 0 


0 0 


0 0 


0 0 


+0,082 
+0,282 
+0,510 
+0,675 
+0,707 
+0,590 


0 


5/6% 
3/40 
2/37 
"lem 
"ym 


Wee 


Sa cg +0,066]+ 0,052 


+0,224]+ 0,184 
+0,418]+ 0,359 
+0,588]+ 0,536 
+0,700]+ 0,678 
+0,739]+ 0,753 


+0,071}+ 0,035 
+0,245]+ 0,135 
+0,454]+ 0,293 
+0,634]+ 0,486 
+0,740]+ 0,667 
+0,753]+ 0,779 


+0,0711+0,020 
+0,248]+ 0,093 
+0,469}+ 0,238 
+0,667]/+0,445 
+0,7811]+0,658. 
+0,779}+0,799 


+0,068]+ 0,007 
+0,244/+ 0,056 
+0,473]+0,189 
+0,686]+ 0,409 
+0,810]+ 0,650 
+0,799]+ 0,815 


+0,065]—0,004 
+0,238]+ 0,022 
+0,472]+0,145 
+0,699]+ 0,375) 
+0,833]+ 0,641 
+0,815]+ 0,827 


+0,062 
+0,230 
+0,468 
+0,707 
+0,850 
+0,827 
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0 0 


0 0 


0 0 


0 0 


0 0 


0 


5/g% 
3/40 
?/3.7 


"1% 


1 
/s a 


Se est +0,069 + 0,065 


+0,239]+ 0,227 
+0,4421+ 0,430 
+0,621]+ 0,620 
+0,739)+ 0,753 


1/, 7 |+0,780}+0,800 


+0,065}]+0,055 
+0,227}+0,201 
+0,430}+ 0,404 
+0,620}+ 0,616 
+0,753]+0,779 


+ 0,055}+ 0,045 
+0,201)+0,175 
+0,404]+ 0,377 
+0,616]+ 0,609 
+0.779]+0,799 


+0,045]+ 0,036 
+0,175}+0,151 
+0,377/+ 0,351 
+0,609]+ 0,600 
+0,799}+0,815 


+0,800}+ 0,840; + 0,840 +0,873) + 0,873 +0,900 


+0,036]+ 0,027 
+0,151]+0,128 
+0,351]+ 0,325 
+0,600}+ 0,590 
+0,815}+ 0,827 
+0,900}+ 0,922 


+0,027 
+0,128 
+0,325 
+0,590 
+0,827 
+0,922 


\ 


Abb. 4. Naherungsweise bestimmte erste Eigen- Abb. 5. Naherungsweise bestimmte erste Eigen- 
funktion nach (25) und (26), ®—=0. Schnittkurve funktion, = 11,415, Richtung der Schubkrafte 
der in den Abb. 4 bis 9 wiedergegebenen Figuren ation Random ok oe 


Abb. 6. Naherungsweise bestimmte erste Abb. 7. Naherungsweise bestimmte erste 
Eigenfunktion, # = 19,349 Eigenfunktion, # = 24,016 
5* 
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Abb. 8. Naherungsweise bestimmte erste Abb. 9. Naherungsweise bestimmte erste 
Eigenfunktion, # = 27,411 Eigenfunktion, # = 30,112 (Beulfigur) 


in Analogie zu (27) gebildeten Funktionen w{ (x, y;@) sind unter Benutzung von 
Tabelle 2 in Tabelle 8 tabelliert. Die Abb. 10 bis 12 enthalten das gewiinschte Ergebnis 
in graphischer Form. Wie aus Abb. 10 zu sehen ist, besitzt bereits die mit dem betrach- 
teten Oberton schwingende schubfreie Platte eine sogar geschlossene Knotenlinie. Die 


Tabelle 7 
Asi, (0) = 
o A®,, (8) oO, (6) A®,. (#) A®, (3) 
0 + 0,032 — 0,704 0 — 0,094 
7,958 + 0,044 — 0,673 —0,114 +. (),284 
20.481 + 0,078 — 0,614 —0,185 + 0,454 


entsprechende von 8. Iguchi in * angegebene Figur geben wir in Abb. 13 wieder. Beide 
Figuren weichen insbesondere im Verlauf und in der Umgebung der Knotenlinien von- 
einander ab. Wir haben diesen Sachverhalt durch eine von der in der vorliegenden Arbeit 
beschriebenen Rechnung vollig unabhangige Rechnung iiberpriift und sind zu der Uber- 
zeugung gekommen, da unsere in Abb. 10 gegebene Figur den Dingen am nachsten 
kommt. Daf der hier untersuchte Oberton mit zunehmenden Schubkraften an den 
Randern im Gegensatz zum Grundton hoher wird, wurde bereits unter Punkt IV ange- 
deutet. Diese Tatsache steht offensichtlich in engem Zusammenhang mit der aus den 
Abb. 11 und 12 ersichtlichen weniger anschaulichen Deformation der Knotenlinie bei 
zunehmenden Schubkraften. 
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Tabelle 8 
pee at CS a aS SEE et ea a a 
(5) . 
i] be! wz Bees wi (a, y); d= 7,958 w® (x,y); 0 = 20,481 
1 2 1, 2 1 2 1 2 
a a et is a I See (a 
0 0 0 0 0 0 
0 I 0 0 0 0 0 
Mio Wl 9% — 0,035 — 0,028 — 0,043 — 0,027 — 0,050 
Hees ects — 0,106 — 0,082 — 0,130 — 0,074 — 0,152 
Uso 9% yx 3/47 —0,154 —0,115 —0,190 — 0,096 20,217 
1/5 a 2/5 7 —0,155 O.0IE —0,188 — 0,079 = 0,208 
Was ate Ot: — 0,128 — 0,094 — 0,143 — 0,056 — 0,134 
My 7 Pal EL —0,112 —0,101 —0,101 —0,071 —0,071 
0 I 0 0 0 0 0 
Sc 1Yf5 JG — 0.106 —0,082 — 0,130 —0,074 — 0,152 
1/7 5/4 7 2.314 — 0,234 —0,389 — 0,205 — 0,449 
bieset 4% 31a — 0,445 — 0,319 — 0,556 — 0,254 — 0,634 
pss 2/3 I — 0,426 — 0,286 — 0,530 — 0,187 —0,577 
5/49 7 "hy I — 0,320 0,217, Seyi —0,108 —0,355 
1/. x Ll, x — 0,265 — 0,234 — 0,234 —0,155 —0,155 
0 I 0 0 0 0 0 
Vie ae py ke — 0,154 —0,115 — 0,190 — 0,096 —0,217 
1], 5/4 7 — 0,445 —0,319 — 0,556 — 0,254 — 0,634 
1,7 Lym 8/4 7 — 0,594 — 0,395 — 0,765 0a — 0,869 
Ls ee Ls ep — 0,491 — 0,280 — 0,660 —0,130 — 0,738 
at Coste — 0,269 — 0,122 — 0,362 —0,017 — 0,366 
ar ee Sts — 0,157 — 0,129 —0,129 — 0,048 — 0,048 
0 mt 0 0 0 0 0 
1% | 14/12 —0,155 =—0,111 —0,188 —0,079 — 0,203 
1/5 2 5/5 — 0,426 — 0,286 — 0,530 — 0,187 =), 57% 
ea Uy Rae 3/4 % — 0,491 — 0,280 — 0,660 —0,130 — 0,738 
5 a Cs ae — 0,256 — 0,040 — 0,431 + 0,109 — 0,514 
ON esis Cie + 0,106 + 0,235 —0,011 + 0,332 — 0,062 
es ate 1), 0 + 0,269 +0,277 + 0,277 + 0,296 + 0,296 
0 It 0 0 0 0 0 
Loe PLAS — 0,128 — 0,094 — 0,143 — 0,056 — 0,134 
Ts 5 5/_ 7 — 0,320 —0,217 — 0,371 — 0,108 — 0,355 
Hi OLE Lhe wre — 0,269 — 0,122 — 0,362 —0,017 — 0,366 
Seas 2/9 7 + 0,106 + 0,235 —0,011 + 0,332 — 0,062 
ile sats "le I + 0,572 + 0,628 + 0,473 + 0,635 + 0,387 
1]. om 1]. 3 +0,781 +0,748 +.0,748 + 0,679 + 0,679 
* 0 It 0 0 0 0 0 
Tie tt ee —0,112 —0,101 —0,101 — 0,071 —0,071 
gaa Ie Be — 0,265 — 0,234 — 0,234 — 0,155 — 0,155 
a Ll it Sage — 0,157 — 0,129 — 0,129 — 0,048 — 0,048 
hs 7 Ws 4 
1/5 2 2/5 7 + 0,269 +0,277 + 0,277 + 0,296 + 0,296 
OW iy 6 Hepege +0,781 + 0,748 + 0,748 + 0,679 + 0,679 
12 12 ? 
weece OS + 0,999 + 0,957 + 0,957 + 0,849 + 0,849 


Abb. 10. Naherungsweise bestimmtedritte Eigen- Abb. 11. Naherungsweise bestimmte dritte Eigen- 
funktion des untersuchten Teiles nach (25) und  funktion, ®= 7,958. Richtung der Schubkrafte 


(26), d=0 an den Randern: i 


Abb. 138. Vergleichsfigur 
zu Abb. 10 von 8. Iguchi aus 4 


Abb. 12. Naherungsweise bestimmte dritte Eigen- 
funktion, ® = 20,481 
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VI. Zum Ausbeulen quadratischer Platten* 


Unser Ziel ist zunaichst eine Verbesserung der unter Punkt IV erhaltenen Naherung 
9 (0) fiir den Beulwert %;,;,. Wir wollen uns hierbei, wie in der Kinleitung erwahnt, 
der zu A’ gehorigen mittels des vollstiindigen Orthonormalsystems von auf 1 normierten 
zugeordneten Kugelfunktionen {J// (w) - IZ‘ (y)} geschlossen darstellbaren in !** 
angegebenen Kernmatrix (a’,,,;,) bedienen. Zu beriicksichtigen ist, daB unsere in (3) 


definierten Parameter 9 und 4 mit den zu A’ gehorigen Parametern # und 1’ wie folgt 
in Beziehung stehen: 


i (e 16 <, 


In Analogie zu (25) geben wir in (29) den zahlenmaBigen Ausdruck fiir den neunten 
Abschnitt ,(@’7)..) ((m —n) = 0, (s — t)=0 (mod 2); m, n, s, t = 4) des mabgeb- 
lichen Teiles der Kernmatrix (a@’»,,;;) an. Die zahlenmaBige Ermittlung der Koeffizien- 
ten konnte wiederum ohne Benutzung von Rechenmaschinen erfolgen. Alle Koeffizien- 
ten sind im Gegensatz zu (25) in (29) ohne irgendwelche Rundungen eingetragen. Fiir 
die Wahl des benutzten Orthonormalsystems von zugeordneten Kugelfunktionen 
spricht auBerdem, da in (29) vier Koeffizienten verschwinden. Die gesuchte Naherung 
&#°) (0) fiir 0;,;, erhalt man, indem man die zu (29) zu bildende Sakulardeterminante 
und den dort auftretenden Parameter /’ gleich Null setzt. Die genannte Determinante 
1aBt sich aus Symmetriegriinden diesmal in eine sechsreihige Determinante tiberfiihren. 
Durch passende Multiplikationen von Spalten und Zeilen sind simtliche auftretenden 
Wurzeln zu beseitigen. Auf dieselbe Weise laBt sich erreichen, daB lediglich die Glieder 
eines zweireihigen Abschnittes der Determinante mit 1/9’? behaftet sind. Die Anwendung 
des Eliminationsverfahrens fiihrt auf eine quadratische Gleichung in #’?. In Tabelle 9 
vergleichen wir den erhaltenen Naherungswert #’) (0) fiir den Beulwert #’;,;, nach 
passender Umrechnung mit den von anderen Autoren angegebenen Werten. Die dort 
aufgefiihrten Werte von S. Tomotika, K. Hidaka, 8. Tomotika und I. Imai und 
K. Munekata wurden der bereits erwahnten Arbeit von M. Hasegawa entnommen. 
Zu den von den ersten drei Autoren benutzten Rechenverfahren kann nicht Stellung 
genommen werden, da die Originalarbeiten nicht zuganglich waren. Der Wert von 
S. Iguchi findet sich in § und. 8S. Iguchi verallgemeinert seine bereits auf den Fall der 
freischwingenden, schubfreien Platte in * angewandte Entwicklungsmethode auf den 
hier vorliegenden Fall. Der angegebene extrapolierte Wert wurde aus Abschnittsergeb- 
nissen der betrachteten Determinantengleichung mit Hilfe der Newtonschen Inter- 
polationsformel ermittelt. Den Untersuchungen von W. Moheit'+ ” lag das Ditffe- 
renzenverfahren zugrunde. Die Extrapolation auf den Fall einer unendlichen Deter- 
minante erfolgte mittels eines Potenzgesetzes. R. C. T. Smith untersuchte orthotrope 
Platten mit Hilfe des Rayleigh-Ritzschen Verfahrens. Als Spezialfall findet sich in ¥ 
der in Tabelle 9 aufgefiihrte Wert fiir isotrope Platten. Bei diesen Untersuchungen 
wurden dieselben Funktionen verwendet, die S. Iguchi in * und ?° fiir seine Entwick- 
lungsmethode benutzte. Ein Vergleich der beiden Verfahren fiel zugunsten des von 
R.C.T. Smith benutzten aus. B. Budiansky und R. W. Connor verwendeten in“ zur 
Ermittlung einer oberen und einer unteren Schranke fiir den Beulwert die Lagrangesche 

* Wir wollen an dieser Stelle nochmals darauf hinweisen, da durch ein Versehen bei der Druck- 
legung von ! die beiden mafgeblichen Beulfiguren 2 und 3 vertauscht worden sind. Unsere daselbst 
gemachten Ausfihrungen tiber Eigenwerte und Kigenelemente mogen auch zur Erganzung der 
vorziiglichen Darstellung von C. F. Kollbrunner und M. Meister 6 (insbes. 8S. 180) dienen. Eine wei- 
tere Bemerkung zur Deutung der verhaltnismabig groBen Abweichung des Ergebnisses von S. Iguchi, 
verglichen mit den iibrigen in * S. 239 zusammengestellten Werten, wurde von R. C. T. Smith in’ 
gemacht. 
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Multiplikatormethode. Es zeigt sich bei Betrachtung unserer Tabelle 9, da% der von 
uns gefundene, als obere Schranke fiir den Beulwert anzusehende Wert 14,6395 unter- 
halb der von den beiden Autoren zu 14,64 angegebenen unteren Schranke liegt. Der 
auftretende Unterschied diirfte bei Benutzung eines groferen Abschnittes unserer 
Kernmatrix unter Hinzuziehung elektronischer Rechenmaschinen mit Leichtigkeit zu 
verscharfen sein. Wir haben die von uns angestellten Rechnungen unter Beriicksichti- 
gung einer groBen Zahl von Dezimalstellen mehrmals auf verschiedenen Wegen iiber- 
priift und sind zu der Uberzeugung gekommen, da der von uns angegebene Wert 
zuverlassig ist. Damit ist durch Gegenbeispiel gezeigt, daB der Wert 14,64 nicht untere 
Schranke fiir den Beulwert sein kann. 

Zur Abrundung des sich hinsichtlich des Ausbeulens quadratischer Platten ergeben- 
den Bildes stellen wir zur Verdeutlichung der Einspannungseinfliisse in Tabelle 9 den 
Werten fiir die an den Randern fest eingespannte quadratische Platte den von 
St. Bergmann* und H. Reissner in ® gegebenen analogen Wert fiir eine drehbar 
gelagerte quadratische Platte (Naviersche Randbedingungen) gegeniiber. Es zeigt sich, 
da} die drehbar gelagerte Platte bereits bei einem um 36% geringerem Wert im Vergleich 
zur eingespannten Platte ausbeult. 

Wie unter Punkt V angekiindigt, stellen wir nun in Abb. 14 simtliche bisher bestimm- 
ten Beulfiguren fiir die an den Randern fest eingespannte, quadratische Platte zum 
Vergleich mit unserer in Abb. 9 gegebenen Figur zusammen. Es handelt sich um Figuren 
der Autoren S. Iguchi ® 1, R.C.T. Smith ® und W. Moheit 1 1. Auf eine Neu- 


Abb. 14. Vergleichsfiguren (Beulfiguren) zu Abb. 9 von S. Iguchi nach ® und , R. C. T. Smith 
nach 13 und W. Moheit nach 1! und 


* Der Verfasser ist Herrn Prof. Dr. St. Bergmann, Stanford, fiir freundliche Hinweise und fi 
die Uberlassung von Sonderdrucken sehr zu Dank verpflichtet. 
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berechnung unserer Figur mit Hilfe des Kernmatrizenabschnitts (29) haben wir ver- 
zichtet, da insbesondere mit den Figuren von S. Iguchi und R. T. C. Smith gute 


Ubereinstimmung besteht. 
In Abb. 15 stellen wir zum Abschlu8 unserer 


Ausfiihrungen die von St. Bergmann und H. Reissner 
in 5 angegebene Beulfigur fiir die drehbar gelagerte, 
quadratische Platte den Beulfiguren fiir die fest ein- 
gespannte, quadratische Platte aus den Abb. 9 und 
14 gegeniiber. Die zuerst genannte Figur ist gegen- 
iiber den anderen Figuren insbesondere durch eine 
fast das ganze Platteninnere ausschépfende Beule und 
einen grundsatzlich anderen Verlauf der Knotenlinien 
ausgezeichnet. 
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Ein einfaches Wartezeitproblem bei einem Poissonschen Verkehrsflu8 


Von H. Kremser, Wien 
Mit 3 Textabbildungen 


Zusammenfassung. Es wird die Verteilung der Wartezeit berechnet, die verstreicht, bis bei 
einem Poissonschen ProzeB erstmalig eine Zeitliicke >¢ auftritt. Das Problem ist von Interesse 
fir das Einbiegen von Kraftfahrzeugen in einen Verkehrsstrom und fiir das Kreuzen einer Vorrang- 
straBe. Die Ergebnisse koénnen Verwendung finden zur Entscheidung der Frage, ob eine Kreuzung 
geregelt werden soll oder nicht. 


I. Einleitung 


Betrachtet man an einer FernverkehrsstraBe mit nicht zu dichtem Verkehr die 
zeitliche Aufeinanderfolge der Fahrzeuge, so kann man diesen Vorgang in guter Nahe- 
rung als Poissonschen Proze8 beschreiben. Die Wahrscheinlichkeit dafiir, daB in einem 
Zeitintervall der Lange ¢ die Anzahl N (t) der passierenden Fahrzeuge gleich m ist, ist 
also gegeben durch die Formel: 


(¢ t)m “ect 
m! 


PNG) = y= 


Die Anwendung dieses Modells kann man rechtfertigen, indem man nachweist, da 
die folgenden kennzeichnenden Bedingungen fiir einen PoissonprozeB erfiillt oder 
wenigstens naherungsweise erfiillt sind (vgl. [3], [4]): 
a) Die Anzahlen der Fahrzeuge, die in sich nicht tiberschneidenden Zeitintervallen 
eine bestimmte StraBenstelle passieren, sind unabhangig verteilt. 
b) Die Wahrscheinlichkeiten fiir das Auftreten einer bestimmten Anzahl von Fahr- 
zeugen in Zeitintervallen vorgegebener konstanter Lange sind gleich. 
ec) In kleinen Zeitintervallen ist die Wahrscheinlichkeit fiir das Passieren eines Fahr- 
zeuges proportional zur Lange des Zeitintervalls: 
P(N (4t) = 1)=c-At, tim A=) _ 
At— 0 
d) In kleinen Zeitintervallen treten zwei oder mehr Fahrzeuge kaum aut (an einer be- 
stimmten StraBenstelle wird man also selten Uberholungen feststellen kénnen): 
Piven a3) 70 (Ai Minn ES 2 26) 


At—+>0 


Diese Bedingungen sind auf Fernverkehrsstrafen mit nicht zu dichtem Verkehr an- 
nahernd erfiillt. Untersuchungen dariiber sind in [1] enthalten (vgl. ferner die in [1] 
angegebene Literatur). 

Miindet an der Stelle der FernverkehrsstraBe (HauptstraBe), an der der Verkehr 
beobachtet und als PoissonprozeB beschrieben wird, eine NebenstraBe, so wird ein 
Einbiegen von derselben in die Hauptstrafe nur méglich sein, wenn im HauptstraBen- 
verkehr eine geniigend groBe Liicke eintritt. Dasselbe gilt fiir das Kreuzen des Haupt- 
straBenverkehrs. Es ist daher die Verteilung der Wartezeit von Interesse, die verstreicht, 
bis bei einem Poissonschen ProzeB erstmalig eine Zeitliicke gegebener Mindestlange auf- 
tritt. Man kann sich von den Verhaltnissen etwa folgende genauere Vorstellungen 
machen: 

Bei der Einmiindung der NebenstraBe befindet sich eine Stoptafel. Die Kinreihung 
kann nur in die rechte Fahrtrichtung erfolgen. In der Nahe der Einmiindung gilt auf 
der HauptverkehrsstraBe ein Uberholverbot. Ein NebenstraBenfahrzeug wartet so 
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lange vor der Stoptafel, bis im Hauptverkehrsstrom eine Liicke eintritt, die dem 
Fahrer als ausreichend zur Durchfiihrung der Einordnung erscheint. Diese Vorstellungen 
entsprechen der Zufahrt auf eine Autobahn. Fiir die Autobahn kommt natiirlich nur 
Einordnung, aber nicht Kreuzung in Frage. 

Ist auch eine Einordnung in die linke Fahrtrichtung von Interesse, so kann man 
annehmen, daB die Linksabbieger sich vor der Stoptafel bereits auf der mittleren Fahr- 
spur der NebenstraBe befinden. Man kann die Links- und Rechtsabbieger dann getrennt 
betrachten. Der linksabbiegende Fahrer wartet auf eine hinreichend lange Zeitliicke in 
einem ProzeB, der eine Summe zweier unabhangiger Poissonprozesse darstellt, deren 
mittlere Zeitliicken — , + = 1, 2 im allgemeinen voneinander verschieden sein werden. 
Dieser SummenprozeB ist aber wieder ein Poissonprozels mit der mittleren Zeitliicke 
pee Man kommt also auf das vorher geschilderte Problem zuriick. 

1 2 

Untersucht man das Kreuzen des Hauptverkehrsstroms, so betrachtet man ebenfalls 

die Liicken des Summenprozesses. 


Man kann diese Aufgabe als ,,einfaches‘‘ Wartezeitproblem bezeichnen zur Unter- 
scheidung von einem ,,zusammengesetzten‘‘ Wartezeitproblem, bei dem auch die 
Warteschlange der Fahrzeuge vor der Kreuzung beriicksichtigt wird. Ob man bei der 
Berechnung der mittleren Wartezeit das einfache oder das zusammengesetzte Problem 
zugrunde zu legen hat, hangt vom Verhaltnis der Verkehrsdichten auf der Haupt- 
und der NebenstraBe ab. 

Die Arbeit [1] enthalt naherungsweise Lésungen dieser Probleme. Mégen diese 
Naherungen fiir die praktischen Erfordernisse durchaus geniigen und dariiber hinaus 
wegen ihrer EKinfachheit der exakten Lésung vorzuziehen sein, so ist es doch von Inter- 
esse, dafi man auch die exakte Verteilung der einfachen Wartezeit mit einem tragbaren 
Rechenaufwand ermitteln kann. Das ist der Gegenstand dieser Arbeit. Interessanter- 
weise sind iibrigens die exakten Formeln fiir die ersten beiden Momente — insbesondere 
also fiir die mittlere Wartezeit — einfacher als die entsprechenden Formeln der genaher- 
ten Rechnung, obwohl die Verteilung selbst nur rekursiv berechnet werden kann. Ohne 
Zweifel kann man die Berechnung dieser Verteilung bei Anwendung weittragenderer 
mathematischer Hilfsmittel kirzer gestalten. Doch besitzen bei mathematischen 
Themen, die fiir die Anwendungen bedeutsam sind, auch elementare Methoden ihre 
Vorziige. Uber die Verteilung der zusammengesetzten Wartezeit hofft der Verfasser 
in nachster Zeit ebenfalls eine Arbeit verdffentlichen zu kénnen. 


In der Arbeit [1] findet man auch die Ergebnisse empirischer Untersuchungen zu 
diesen Fragen und weitere Literatur. Die grundlegenden Tatsachen iiber den Poisson- 
prozeB, kann man beispielsweise den Biichern [2], [3], [4], [5] entnehmen. 

Die Rechnungen dieser Arbeit kann man dadurch iiberpriifen, da man auf einem 
Rechenautomaten einen Poissonprozef} ablaufen laBt und die zugehérige ,,empirische‘‘ 
Wartezeitverteilung ermittelt. Der Vorteil dieses Vorgehens liegt darin, da&8 man ohne 
groBen Aufwand sehr groBe Stichproben erhalt. Die Ermittlung wirklich empirischer 
Wartezeitverteilungen zu groBen Stichprobenumfangen st6B8t auf Schwierigkeiten, die 
nicht nur in den Kosten der Erhebung bestehen. Man miif&te dann namlich die Beob- 
achtung tiber langere Zeitintervalle erstrecken, und die Annahme, dafi der VerkehrsfluB 


ein PoissonprozeB mit konstanter mittlerer Liicke — ist, ware dann sicher ungerecht- 


fertigt. Die obige Methode wurde an dem Elektronenrechner IBM 650 des mathemati- 
schen Labors der Technischen Hochschule Wien angewendet*. Die Ergebnisse sind in 


* Die Durchfihrung wurde erméglicht durch verstiandnisvolle Unterstiitzung durch Herrn 
Prof. Dr. R. Inzinger, wofir ich auch hier meinen Dank aussprechen méchte. 
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Abb. 1 wiedergegeben. Die Hauptergebnisse der Arbeit sind in den folgenden Formeln 
enthalten : 


(59), (60), (73), (74): Wartezeitverteilung. 

(70): mittlere Wartezeit. 

(72): Dispersion der Wartezeit. 

(79), (81), (83), (85), (87): mittlere Wartezeit zu Fahrzeugkollektiven. 


Die Abbildungen 2, 3 geben einen Vergleich der Ergebnisse mit den entsprechenden 
Naherungen der Arbeit [1]. 


Abs. Haufigkeiten 


Simulation 
| Theorie east 
500 | 
| ss 
| Say 
Se 
N 
x 
ri 
SS 
7. 
+— ++ + EEE EEE —— | $+ + + + -¢ 
0 OS 1.0 16 


Abb. 1. Vergleich absoluter Haufigkeiten (c= 1, t= 0,5): 
a) ,,.Empirische* Haufigkeiten, erhalten durch Simulation auf einem Rechenautomaten. 
b) Theoretische Haufigkeiten, angenahert durch die Funktion N - 9 (a, t) - Ao (Stichprobenumfang 
N = 34719, Klassenbreite do = 0,05, Wahrscheinlichkeitsdichte » (o, t) fiir t= 0,5) 


II. Berechnung der Wartezeitverteilung 


Es wird ein Poissonscher ProzeB betrachtet und die Verteilung der Wartezeit w 
ermittelt, die verstreicht, bis erstmalig ein Ereignis auftritt, auf das eine ereignislose 
Zeitliicke = t folgt. Dabei bestehen offensichtlich die einander ausschlieBenden Méglich- 
keiten, daB diese erste Zeitliicke =¢ vor dem ersten Ereignis oder nach dem ersten, 
zweiten, ..... n-ten Ereignis auftritt. Dementsprechend wird definiert (P (M) Wahr- 
scheinlichkeit des Ereignisses 7): 

Ha (t2) = P(i<t A is REET BE aie mine), 
ial (1) 
HA, (t,t) = P (t, 2 t), 


worin t; die Zeitliicke zwischen dem (7 — 1)-ten und dem 7-ten Ereignis bedeutet (¢, Zeit 
bis zum ersten Ereignis). Es gilt dann: 


Pat) = as Gert (tet). (2) 


In den obigen Ausdriicken kann man fiir n> 0 statt > auch =, statt < auch < 
schreiben, weil bei diesen Abainderungen nur Ereignisse mit Wahrscheinlichkeit Null 
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unberiicksichtigt bleiben. Zur Berechnung der Wahrscheinlichkeiten H, (t, t) bendtigt 
man die Wahrscheinlichkeitsdichte der Zeitliicken (¢,, ... tm), die durch 


Chie (fiir x; = 0) (3) 


gegeben ist. Durch diese Dichte kann man sich namlich den Poissonschen ProzeB defi- 
niert denken (vgl. [2], S. 403). Daraus ist insbesondere auch die Unabhangigkeit der ¢; 
zu erkennen, aus welcher sich fiir H, (t, r) sofort die folgende Beziehung ergibt: 


Bil) = P(h<t Wo Bist) PGuet= 


i=1 (4) 
= Ob re 
mit 
a ea Cee Suxd). (5) 
7=1 
(5), (2) fiihrt auf 
Pigs = err yO, @, ze) (6) 


Um die Wartezeitverteilung zu ermitteln, mu8 man also die GréBen Q, (t, t) berechnen. 
Dies erfolgt in rekursiver Weise. 

Die Berechnung der Q, (t,t) gelingt ohne Schwierigkeit fiir + => n¢. Dann folgt nam- 
lich aus ¢, <7, ...% <t die Beziehung 


weshalb gilt: 


Ce es Sih Pe a aya By PS ea ae ae) 


also: 
OQ, Ain T) c= Cb — 6. ine eee ia. (7) 


Fiir beliebige t ergibt sich in gleicher Weise: 
On (6.2), SAL = e)", (G2) 


woraus die gleichmaBige Konvergenz der Reihe (6) in t und fiir alle t < ¢, folgt. 
Fir0<t< nt erfordert die Berechnung der Q, (¢, 7) mehr Mihe: 


n 


Q, (t, 7) =“) a { oP + exp (- oy «| [ax (8) 


t=1 


n 
OSn<1,0S Sy ast 
—o 


Fiir n = 2 wird nun in folgender Weise vorgegangen: 


Die Integrationsvariablen (x, ... %,) =z werden einer orthogonalen Transfor- 
mation ) = 2x unterworfen, wobei speziell 
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gesetzt und YY sonst nicht genauer ae wird. Der Integrationsbereich 


OAS a; = t, Rea YS ans 
n 


MSS ae 
a 


geht dabei iiber in OS 4.5 = 
n 


(Cen SEI. 


worin die Menge M nicht genauer beschrieben werden mu8, da die Integration iiber 
Yo, --. Yn Offenbar den Inhalt von 


(ai: Sai = ys |n) (m0 5% < ti = | een 3) 


ergibt. Bezeichnet man diesen (n — 1)-dimensionalen Inhalt mit J, (y, \/n n,t), so gilt 
wegen Det 2 = 1 und nach Venue dnng, der Substitution \/n 2 Yiu 


Qu (t,t) =f... [orem dy, ff dye = 


1=2 


ee, 


ie 

=fos n (Ys|/n, t)e- clan. dy, = ris J, (a, t)e-°? do fiir n = 2. 
e 2 
0 


S 
= 
ford 
Se 


Darin ist J, (o, t) der Inhalt von 


n 


tet) = Cae == 0) Was Oe = t,t I ae eee) 
Die rekursive Fee eet ie 8 (t, ¥ ist dadurch zuriickgefiihrt auf die rekursive Be- 
rechnung der Inhalte J, (c, t). Jn+1 (0,0) ist definitionsgema8 der Inhalt der Punkt- 
menge 
M,+1(6,t): (11) 
Oe w thea os Beth 


Indem man x, +1 durch x*, +; = 0 ersetzt, kommt man zu einer Projektion JN, (a, ¢) der 
obigen Menge M,,+: (0, #): 


Sa, =.s; 02-3 = G; 
N, (oe, t): 4 = (12) 


ee ech Moe la the 


Betrachtet man jenen Teil der Menge N, (a, ¢), der sich fiir einen speziellen Wert s 
ergibt, so erhalt man offensichtlich J, (s, t): 


| ge (13) 


i=1 


ee 
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Der Inhalt von M, (s, t) ist J, (s, t). Daher ist der Inhalt von N, (o, t) gegeben durch: 


o 


te 
TNs ah =) 4 Max; = {7 [M, (y; Vn, t)] dy, = 
ed = 0 (14) 
lr : 


— 1] 
= | Jn (ys Yn, 8) dy = if Toate ae 
0 


0 


Zu (14) ist prazisierend noch folgendes zu sagen: Der Integrand J, (s, ¢) kann fiir Teile 
des Integrationsbereiches [0, ¢] verschwinden. Das tritt naimlich dadurch ein, daB die 
Menge M, (s, t) fiir entsprechende s-Werte leer ist. AnschlieBend werden die s-Bereiche 
ermittelt, fiir die VM, (s, t) nicht leer ist. Beziiglich o sind dabei Fallunterscheidungen vor- 
zunehmen. 


03,0 St: Man kann n — 1 der Koordinaten 2; gleich Null setzen und die n-te 
Koordinate gleich s. Dadurch ist also 
yy) “LS (15) 
ee 


fiir beliebige s mit 0 < s <o erfillbar und es gilt dann 
t+1=—d0— $§, 0 SUH ye SS Ue 
Die Grenzen lauten somit hier: 


0<s <a. (16) 


t=6 3 2k: (15) ist fir s <o—#, also s=o—t—e (e>0)=nicht erfillbar. 
Aus 
n+1 n 
Du=o, Sy So Sta 
i=1 ee 


folgt namlich 
Un+1 = t ++ E, Unt+1 => ihe 
Fir 
6—ts8 SG (17) 
sind alle Bedingungen erfiillbar, wenn n = 2 ist: Man kann n — 2 der 
Koordinaten 2; gleich Null setzen und die restlichen zwei x; gleich s t 


Diese beiden x; liegen dann zwischen 0 und ¢ und dasselbe gilt fiir x, +1: 


OESep | 90 == Se ate (18) 


ttso <(¢+1)-tfirn 21+ 1: Wie vorher ergibt sich, daB (15) fir s<o—i 
nicht erfillbar ist. Fiir 
o—t<s<o * PGES) 
sind alle Bedingungen erfiillbar: Man setzt » — 1 — 1 der Koordinaten 
x; gleich Null und die restlichen 7 + 1 Koordinaten gleich ee 
Wieder gilt dann auch (18). 


= ts 
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nt<o<(m-+ 1)t: Wie oben folgt, daB (15) fiir s<o—t nicht erfiillbar ist. 
s mit o —t<-s< nt sind realisierbar: Man setzt 


8 


Eso os Wie oh: tee enti =9—8, (20) 


weshalb jedenfalls 

Oe 2 (20) 
gilt. Zu zeigen bleibt die Richtigkeit dieser Beziehung fiir i = n +1: 
AUS Gta ss ni-folet s=a—ttk nteso2 So 4 
+k [(n+ 1) —co] fir ein k mit 0< kl. 


Aus nt<o<(n-+ 1)t folgt c=nt+4 xt fiir ein x mit 0<x< 1. Daraus 
ergibt sich: s=(n—1)t+x“t+kt—kxet, o—s=t[1—k(1—~»)], 


also 
ORGS = eae 
s>nt sind nicht realisierbar, da aus 7; <t,i=1,... mn folgt: 
ee (21) 
Gant: Wie in (21) folgt hier J, (s, t) 7 fiir alle s mit O< s<o. 


Formel (14) kann nun folgend prazisiert werden: 
b 


J (AZo, t)) = [Jn (s, t) ds, (22) 
in 
is ge fir 02 6274 
a OE betioefitit2 oe Ge lin si +1 (23) 
re nt firnt<os (n+ 1)t. 


Formel (22) wird nun durch die folgende Uberlegung zu einer Rekursionsformel 
fiir J, (co, t) ausgebaut: 
: ae! 
M,+1(c, t) ist Teil der Hyperebene » x; =o und schlieBt mit seiner Projektion 
i=1 


N, (co, t) einen Winkel «,+1 ein, dessen Kosinus gegeben ist durch 


1 1 eeaee fet mern ; 
(ie el, bye eh, 2 10 ==) PCOS o, 1, 
aa See ) eee ea I pa (24) 
cos a eae 
Daher folgt: 
J [Mn +41 (o, t) | Nd se (o, t) = cos Paes LN n (o, t) | ee |n+ Ls J [NV (o, t)], (25) 
b 
Jen) 22 yee [J (s, t) ds, (26) 


wobei a,b durch (23) gegeben sind. 

(26) ist die gesuchte Rekursionsformel fiir J, (o, ¢). Zur Vereinfachung des Weiteren 
werden an Stelle von J, (a, t) die GroBen f, (a, t) eingefiihrt. In der Formel (10) wird 
namlich 


fa (a, #) = 2 (27) 


bendtigt. Fir f, (c, ¢) gilt nach (26): 


Ingenieur-Archiv XVI/1 6 
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fn+1 (6, t) = ( fn (8, t) ds (28) 
mit a, b gemaB (23). (10) erscheint nun a der Form: 
QO, 7) = oft (6, i) 6 aa. (29) 
6 
Der Fallunterscheidung fiir o Rechnung tragend, kann man J, (a, ¢), fn (o, t) mit einem 


zweiten Index versehen: 


J (o, t) = er (o, t) 


(30) 
fn (9, t) = fa,j (o, t) 
fir ( —1) ts o<jt,j =1,... », woraus folgt: 
tng (9 — 1) t,t) = fa—-a (GG — N48) (30,1) 
fir 7 <n. ; 
(30), (28), (23) ergibt: 
iru at )= faa (8, t) (31) 
fur) <o' = ¢: 
(EaA)4 
faxay Xe = fies 1 Ce avidaicial hategh (32) 
oat (i=1)t 
fir @— ltsaeatltitan; 
nt 
fda (Ost eats t) ds (33) 


Ci 
fir nt<o< (n+ l)t. 

Daraus konnen die f,,; (o, ) rekursiv berechnet werden. Da die Formeln (9) bis 
(33) fiir n= 2 gelten, sind noch die Falle n = 0, 1 zu betrachten. Bei der Berechnung 
von P (uw, <7) gemaB (2) oder (6) bendtigt man zunachst 

Het, 0): 2 Gye 8) = Cee Oe oe (34) 
Man kann also 
Qo (é, tT) = 1 (35) 
setzen. Ferner gilt: 
fi, (,%) = P(,<t4, ote" =P (= min (7t))*e* =O, ee 


min (t, tT) 


As (t, T) — cfe- do le e-min (ts < 
0 
wesheib man 


fia(o, t) = 1, fie (o, t) = 0 (36) 

setzen kann und (29) dann auch fiir » = 1 gilt. Fiir n = 2 erhalt man: 
1 

(o, ere J [ (ay, 0g: ty 4+ 4, = ) ) (xy, ey 0S 2B <1, OS a= Fe 


t) es 
o fir O..0 2h, 
2t—o fir tso<2t, 


0 fir2t < o, 
woraus folgt, dafS auch (31), (32), (33) fiir n = 1 gilt. 
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Nun ergibt sich durch Induktion: 


1 
fn,1 (G, t) = pT i 
ni ¥ n 
wegen fn+1,1 (6, )=|o ds = ae 
0 
( — n—t1 
hv (t) = ba 
nt 0 
(nt—s)"—-1 peas aE) te Ol” 
wegen fe esa dt nth] nH! ds = he ee - 
a (n+1)t—6 
Inn + (0, t) = 0 (39) 


fiir k > 0 wegen J, (s, t) = 0 fiir o= nt. (38), (39) ergibt die Giiltigkeit von (30, 1) auch 
fiir 7 >.n, wobei laut (36) nur der Fall n = j -- 1 = 1 auszunehmen ist. 

Die bisherigen Ergebnisse werden nun zusammengefaft. 

Wegen (6), (29), (30), (35) gilt: 


(40) 


k—1 09 It 23 t 
ae es CS aS Del hi (G,i)e-"de Doe st face tonne =< i) 
=! g~41)¢ n=k (k-1); 

fiir (b —l)t <1 < bt. 

Der weitere Vorgang ist nun folgender: 

a) Es wird nachgewiesen, da man in (40) Summation und Integration vertauschen 
kann. 

b) Fir die Summen unter den Integralen werden Differentialgleichungen aufgestellt, 
die ihre Berechnung [ohne Ermittlung der f, (c, t)]| ermdglichen. 


a) Die Formel (40) wird fiir & =: 1, 2, ... nacheinander betrachtet. Aus (31), (32), 

(33) folgt induktiv, daB die Integrale 

o° | fa & (a, t).e* do 

(kK —1)¢t 
regulare Funktionen in 1 sind (fiir festes t). Die Reihe, die durch Summierung tiber n 
entsteht, ist nach (7, 1) gleichmaBig konvergent. Sie stellt daher ebenfalls eine regulare 
Funktion von + dar. Ihre Ableitung ergibt sich durch gliedweise Differentiation. Die 
abgeleitete Reihe ist wieder gleichmabig konvergent. Man kann daher in (40) Summa- 
tion und Integration vertauschen: 
| CL a ee eo(i ++ e-*% >’ t.,. (a, 8) to = 

6 n=1 


U1 (41) 


k =a 1 a cs ¥ os 
= ee( ey | eos» ott,4 (o,.t) do + | e-2* > i470, 2) da} 
ve Meee ia Oe Gots gt=* 
fir (k—l)tstsakt. 
(Man kann den Nachweis der gleichmaBigen Konvergenz auch ohne Benititzung des 
; err 1 
WeierstraBschen Doppelreihensatzes erbringen. Ohne Miihe gelingt dies fiir ¢< ~~ auf 


Grund der folgenden Abschatzung: 


1 1 jn-1 
fr, 1 (a, t) Sati, is 


eC =f) > @far(o,t) SoD (ett 


n n 


6* 
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Fir t= = ist — nach Weglassen eines geniigend langen Anfangsteiles der Reihe — 


1 relativ klein gegen n und es diirfte dann eine ahnliche Abschatzung moglich sein, die 
aber wohl etlichen Aufwand an ,,Geometrie“ erfordern wiirde.) 
b) Im Grunde genommen sind also nicht die Funktionen f, (, ¢), sondern ihre 


Summen: 


S (0, t) = d'c* fr (6, #) (42) 
n= 
bzw.: 
Si (a, t) = S7e* fr,1 (6, 8) fir (I — l)tsoslt (43) 


von Interesse: 


LOS e) ae [ + ioe S (oc, thd |= 
: (44) 


1 lt t \ 
=e (145 | e-** §,(a, 1) do + | e-** 8, (a, t) io] 
\ tt it @= ie 
fir (n — l)tst<nt. 
Aus (31), (32), (33), (37), (43) folgt: 


(iN: ee (l—1)t 


Si (a, t) = ott? frase (at) =| f-11-1 (3,8) de+ Sort {fn 1a (8, t)ds + (45) 
pli ot var 4 
6 (—1)t Go 


+ Sort | fi (3, t)ds=c| 
oO 


C— 


DY otfyi-1(s,t)ds tefl Set fa1(s, t) ds, 


n=l (=1)4 3 n=I1—1 d—1)t | 


(i= BD) t o 


Si (9, natalie Sa (s, t) ds ah c| 8 (s, t) ds fiir ja: (46) 
ont Ui 
oder 
Si(o,t 
eS t,t) (47) 


ine (Ae Ih sway 
(I—J)t 


S; ((L— 1)t, t) =e | S;—1 (s,t)ds (Anfangsbedingung). (48) 


(i—2)t 


Aus (30,1) folgen einfachere Anfangsbedingungen fiir (47): 


Si (@—=1)t,8) = Det fr (l= DED) aD fala) e = 


/ 
n=l n=l 


= 8-1 ((L— 1) t,t) — =) freagrey  =4)e, a) = (49) 
= 8-1 (0 — 1) t,t) —e-?} + boa; 
(0;, Kroneckersymbol) 
ASS NE ieee ~) Anfangsbedingungen 
S: (— 1) 4, t) = Si, (@— 18,8) fir 132° (51) fiir (47). 
Aus (45), (47) fclgt durch Induktion: 


Si (o, t) = ef T1 (0, t), (52) 
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worin T 1 (a, t) ein Polynom in o vom Grad / — 1 ist. Im weiteren werden nur mehr die 

Funktionen 7; (oc, t) verwendet. Differentialgleichungen und Anfangsbedingungen fiir 

T; (o, t) ergeben sich unmittelbar aus den entsprechenden Beziehungen fiir S; (c, t): 
eT i (a, t) 


i ae ee Ce a GO edt) furl >. (53) 
1, @,.t) = ¢ (1 —-e=*),’ T, (o, t) =e. (54) 
Pe yh) = Tis (= 1) 4 8) tiie FS (55) 
Batt) —-¢ e-t fc ds + ¢ (l—e-#) =e [1 —e-«! (1+ ce(6—1))] (56) 
o ot 
Pet) = = ¢ et [Ty (st, t) ds F724 (0 — 1) t,t) = oe! | Py 5 (s, t) ds + 
(l—1)t (1— 2) 
Tee ihe, ty, fare 1S. 2) 
pial Fe B 
DP (wiS T) = E°% [ oF oe | T (o,t) do a | 1; (0, t) i (58) 
ee ry ea (n—1)t 


fir (n—1l)t<t< nt. 
Die Wahrscheinlichkeitsdichte ist somit fiir o > 0 gegeben durch 


GG.) = eT a(art). (59) 
Die 7’, (t,t) konnen rekursiv aus (56), (57) berechnet werden: 
Pilot = € 
T, (o, t) =c {1 —e-*' [1 + ¢ (o — #)}} (60) 
Pato, t) == {1 =e HF (1 - et) — e="* |\( — ee) (6 = 21) —e-**- . . 


Die Funktionen f,,; (o,¢) sind nicht negativ [fir (|) — 1)i<o <It], weil sie 
proportional zu Inhalten sind. Dasselbe gilt daher fiir die Funktionen S; (a, ¢), 71 (a, t). 
Die Funktion 7, (c,¢) ist nach (54) konstant. Beim Ubergang zu 7’, (o, t) an der 
Stelle o = ¢ tritt ein Sprung nach unten von der GroBe c e-*' auf. Nach (55) schlieBen 
fiir ] > 2 die Funktionen 7’, (c, t) stetig aneinander. Laut (53) sind sie monoton ab- 
nehmend. 

Die Wahrscheinlichkeitsdichte (59) ist fiir einige Werte des Parameters ¢ in Abb. 2 
dargestellt. Die Formel (58) fiir die Verteilungsfunktion der Wartezeit wird weiter unten 
mittels partieller Integration noch vereinfacht. 

Die Konsistenzbedingung 

P (w= oo) = 1 oder 


nt 


et_j]=)’ T,,(c, t) do 81) 
nae (n—1)t 
ist erfillt: 
nt 
T,(0,t)do = <|7, Git) P41 (ae 1s, ‘)| fare dts (tN 

(n—1)¢ (62) 

aes ed ct £2 = | 
Ps: { T2(6,. dom ct + <| Sato t) = Tri ((n+ Ni t)| = 
(et! n= n=: 


n—1)t 


=ct+ “7, (24,1) =e! —1. 
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So wie eben die Konsistenzbedingung verifiziert oder, anders ausgedriickt, das 0-te 
Moment berechnet wurde, kann man auch die hdheren Momente berechnen. Man braucht 
dazu die Polynome 7,, (a, t) nur fiir die niedrigen n-Werte explizit zu berechnen. An- 
schlieBend wird die mittlere Wartezeit H (w,) ermittelt: 


nt 
E(w) = alo. 1+ Dd) |sT,(s,4 aay (63) 
n= 1(q “1s }} 
Zur Berechnung von 
nt 
i} s 7, (8, t) ds (64) 
(@=—1)% 
wird die mit (57) aquivalente und fiir n= 2, nt<o< (n+ 1)¢ giiltige Beziehung 
ot 
[Pn (st) ds = [Py (nt, ) = Tne (0,0) (65) 
(n—1)t 


mittels o = z+ ¢ umgeformt in 
4 
1, (8,1) ds = — (7, (nt, t) — Trai (z+ #,8)] (66) 
(n—1)¢ 
fir n= 2, (n—1)t<z< nt. Fir (64) ergibt sich dann durch partielle Integration: 


n nt nt z 


t 
67, (8,0) ds=nt | Ty (8,0) ds — [ dz 1\1'.{8,t) ds = 
(n—1)t (n—1)t (n—1)t (n—1)t 


[T, (nt, t) — Tra ((n + 1)t,8)] —{=IT, (nt, t)— Trail +t,8)]dz = 


(n—t)t 


ntest 
c 


(n+1)¢ 
—]) tect tect ct 
Se nh Tir ((m+1),t) +S] Taos (Z, t) dd = (67) 
nt 
(n — 1) tect tect 
at 1) t, t) 
2 ct 
+ [PT a+a ((n + 1), t) — Tra ((m + 2) 8, 4)] 
far ie )2: 
Ct 
fs Ts (8, ids Josds=S, (il), (68) 
0 0 


nt 


E(w) =e-«t « » [s Ts tlds = 


ele Pe hy 
ct e—«t 1 sf : 
a ge om | — WED e (mt, t) — mt Pars (Cm + 1)4, 4) + 
[Pats ((m + 1) ¢, )— Trea (w+ 2)49]| = (69) 
ee ata, ) += 7,(86,9)]. 
E(w) = =~ fet (1-2er)i. (70) 


Durch eine weitere partielle Integration ergibt sich auf gleiche Weise das zweite 
Moment* : 


* Ahnlich wie die Momente kann man auch EF (eiwes) — (¢ — 48) + (ceist — ig ect)-1 h 
d.i. die charakteristische Funktion von wz. ae ) Beige 
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Bat 1 2 e¢ ¢ 
2 2 42 
= Slee - ] — 2'C.t) tot+ SI. 
Aus (70), (71) kann man die Dispersion berechnen: 
D (w,) = (et — 2ete!— 1). (72) 


Mit Hilfe von (66) kann man fiir die Verteilungsfunktion (58) eine einfachere Darstellung 
gewinnen: 


n—1 ai 
Pamsy=ei tet SS —-[ Ti (lt, t) — Trai (b+ 1)4,4)) + 
1=2 
+ *<[ (nt, cee Ty 41 +o) = 


=e (Ltof)t+— [Ps (2t,t) — Teor (+40), 


P< t) = Fy (t,t) =1—— Trai (e +40) (73) 
fir (n— l)t<t<nt,n 2 2. 


Plot) 


as 0.50 


0s Naherung nach [1]: -__—-~~~ 
0.75 


044 


0-3} 


0 1 2 3 4 5 6 


Abb. 2. Vergleich der Wahrscheinlichkeitsdichte mit den entsprechenden Naherungen der Arbeit [1] 


Fir 0<1<t, (n = 1) gilt natiirlich: 
P(wsae y= Fy, t) = e-*' (1+ ct). (74) 
Die Konsistenzbedingung (61) ergibt mittels (73): 
lim Piya (ee 2G 1) = 0 fit (ev — Dita ant, (75) 


n> oO 
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was besagt, daB unendliche Wartezeiten nur mit Wahrscheinlichkeit Null auftreten. 
Die Abbildungen 2, 3 geben einen Vergleich der Wartezeitverteilung und der Mo- 
mente mit den entsprechenden Naherungen der Arbeit [1]. Zur Normierung der graphi- 
schen Darstellung wurde c = 1 gesetzt. Das kommt auf eine geeignete Festsetzung der 
Zeiteinheit hinaus. Die niherungsweise 
Wartezeitverteilung der Arbeit [1] 
unterscheidet sich von der exakten Ver- 
teilung vor allem fiir kleine ¢, und zwar 
hauptsaichlich im Bereich t<o< 21. 
om Die Momente weichen um so mehr 
voneinander ab, je groBer ¢ ist. 


Momente Diy 


Ill. Beriicksichtigung der Verteilung 
der Operationszeiten 


Die oben berechnete Verteilung kann 
auch als bedingte Verteilung 


PikZ oer 
=P wat ee 


aufgefaBt werden, namlich als Vertei- 
lung der Wartezeit unter der Hypo- 
1 2 these, daB die Operationszeit t (= ge- 
Abb. 3. Vergleich der mittleren Wartezeit H (wz) schatzter Zeitbedarf fiir die Durchfih- 
und der Dispersion D (w:) mit den entsprechenden rung einer Hinreihung bzw. Kreuzung) 
Naherungen der Arbeit [1] gleich 7 ist. Fiir den einzelnen Ver- 
kehrsteiinehmer wird oft nur diese be- 
dingte Verteilung von Interesse sein. Wenn dagegen beispielsweise zu entscheiden ist, 
ob eine Kreuzung geregelt werden soll oder nicht, so wird man mit der Gesamtheit der 
Fahrzeuge rechnen miissen und die Operationszeit ¢ wird ebenfalls als zufallige Variable 
aufzufassen sein. Aber auch schon beim einzelnen Verkehrsteilnehmer ist die Operations- 
zeit im Grunde genommen keine determinierte GroBe. Sie ergibt sich ja durch Liicken- 
abschatzung und auBerdem wird sich der Fahrer je nach seiner Gemiitsverfassung und 
der jeweiligen Kondition seines Fahrzeuges eine schnellere oder langsamere Durchfiih- 
rung der Operation zutrauen. AnschlieBend wird nun die Wartezeitverteilung unter 
verschiedenen Annahmen iiber die Verteilung P (¢ < 7) = G(T) der Operationszeit ¢ 
behandelt. Allgemein seien dabei die folgenden beiden Voraussetzungen zugrunde 
gelegt: 


0 


a) Die Operationszeit ¢ sei unabhangig von den Zeitliicken im Verkehr der Vorrang- 
straBe verteilt. 


b) Auf Grund technischer Bedingungen seien Operationszeiten < 7, (7, > 0) 
unmoglich: PG sod) = 0; 


Es gilt dann offensichtlich: 


H(w) =| fr-dG@)- de F(t.) =[F () dG) = [Bw dE) = 
0 0 0 To (76) 
=—E(e—(1+et)) =+H(e) —+_ By. 
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co W ae 

P(w< W)=| |aG()-d,F (t,t)=8,(F (Wd) = » J ROW.nae (77) 
en) n=1 Ww 

wegen F (W, t) = F, (W, t) fiir a he, 


#, bedeutet den Erwartungswert beziiglich t. Fiir 0 < W < 7, gilt wegen b): 
P(w< W)=(1+c W)- E (e-*?). (78) 


Die Reihe (77) ist wegen b) fiir jedes feste W eine endliche Summe. Die Anzahl der 
Summanden nimmt aber mit wachsendem W zu. 

(76), (78) werden nun fiir einige spezielle Verteilungen G (t) ausgewertet: 

1. Kine erste Verallgemeinerung der determinierten Operationszeit stellt eine 
diskrete Verteilung 


P(T) =p, i=1,...N, SY) p-=1, T;= T, (N eventuell co) 


dar: 
1 % 
E(w) = 2) pilee™—(U+e T)] (79) 
i=l 
x 
tee W) ee (l oe W) 2 pe * fur 0= W =a 7,. (80) 
i= 
2. Die Exponentialverteilung 
0 bod ied Bi cea la 
Pits 7) = 
Pee Oe tur 


entspricht der Vorstellung, da durch Verkehrsvorschriften ein Minimalwert 7, fiir 
die Operationszeit festgelegt ist: 


A = ‘igo 
E(w) =sGr=5e" (2; L = =) (81) 
fiir c < A [andernfalls existiert H (w) nicht]. 
P(w< W)=(1+eW)poce'™ (82) 


PO eee Lg Sal 4. 


3. Setzt sich in Verallgemeinerung von 2. das Fahrzeugkollektiv aus mehreren 


i 


3 : : 2 1 
Kategorien mit den Minimaloperationszeiten 7;, den mittleren Operationszeiten — 


und den Anteilen 7; > p= | zusammen, so ergibt sich unter den Voraussetzungen 


(Head f Mtoe Bal kee 
hi 1 1 
Be) =P sas (ele Mig ha (83) 
7 Api ect ... 
PwsW)=(1+e¢eW)-»d Sia fir 0<W< T,. (84) 


4. Die Annahmen 2., 3. sind insofern unrealistisch, als vorgeschriebene Minimal- 
operationszeiten nicht streng beachtet werden. In Wirklichkeit werden auch zu t-Werten 


90 Buchbesprechungen 


mit 7, <t< Tmin positive, wenn auch kleine Wahrscheinlichkeitsdichten gehoren. 


Erlangverteilungen entsprechen — zumindest qualitativ — diesen Umstianden: 
—kh (T —T) 
E | tea Cer ye ed T > Ty. 
ap b= Py 
| lo tind 2 ay 
k ee 
Bw) =+((Ge) e—A+e7)|—7Hhir o< kA (85) 
7 ki i Sy as 86 
PwsW)=(1+eW)-(Gz2) \ fir 0O< WT). (86) 


5. In den Fallen 2. bis 4. wird die Wahrscheinlichkeitsdichte erst in oo gleich 0. 
Der Wirklichkeit niher kommt in dieser Hinsicht wohl die Betaverteilung: 


Oftr-? 245. i ; 
ee Lee 
1 T—T,\? ( T—T,\?2 ar T T<T 
d ae , ‘(q=t So 3] fiir = = 1 
qpitsT)= (T,— Ts) BF, 5) 1 ) ; 1 ft) 
fam ere, 


Fiir m = n = 2 ergibt sich speziell die Gleichverteilung im Intervall [7, 7,], wofir 
(76), (78) zu den folgenden Formeln fihren: 
ecT,1—ecTo 1 To +- 7; 


Ey = ap Ty ng eee (87) 


@—6To — e—edn. 
e(T, — T) 


P(ws W) = (l1+cW) rr OF aT: (88) 
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Stresses in Shells. Von W. Fliigge. Mit 244 Figuren, XI, 499 S. Berlin-Gottingen-Heidelberg: 
Springer-Verlag. 1960. Geb. DM 58,80. 


Das vorliegende Buch ist nicht, wie man vielleicht zundachst vermuten kénnte, einfach eine 
englische Ausgabe des bekannten Standardwerkes ,,Statik und Dynamik der Schalen“ des gleichen 
Verfassers, das 1957 erschienen und bereits wieder vergriffen ist. Es handelt sich vielmehr um ein 
vollkommen neues Buch, das mit dem urspriinglichen nur den grundsatzlichen Aufbau gemeinsam 
hat. Der Umfang ist von 286 auf 500 Seiten gestiegen und die Zahl der Abbildungen hat sich ver- 
doppelt. Die theoretischen Abschnitte sind ausfihrlicher gehalten, die behandelten Beispiele sind 
weitaus zahlreicher geworden. Es wurde auch manches neu hinzugenommen: eine sehr einpriig- 
same zahlenmaBige Gegentiberstellung von Ergebnissen der Membran- und der Biegetheorie, die 
Aufstellung der Biegegleichungen der allgemeinen (nichtkreisférmigen) Zylinderschale und die 
Verallgemeinerung der Biegetheorie der Rotationsschalen auf nicht axialsymmetrische Belastung, 
um nur einiges zu nennen. 
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Weggelassen wurden dagegen die Ausfihrungen tiber Stahlbetonschalen und das Kapitel iber 
Schwingungen. 

Was an dieser Stelle* tiber die deutsche Ausgabe gesagt wurde, gilt in vollem Umfang auch fir 
die englische. Nirgends tiberwuchert der mathematische Apparat, tiberall stehen die anschauliche 
Deutung und die Anwendungsmiglichkeit im Vordergrund. Das Buch ist fiir den Ingenieur geschrie- 
ben und man kann nur wiinschen, da8 auch die hoffentlich bald zu erwartende N euauflage der 
deutschen Ausgabe im gleichen vergréSerten Umfang erscheinen moge. H. Parkus, Wien 


Impulstechnik. Erzeugung und Anwendung von Kondensatorentladungen. Von F. Friingel. 
Mit 256 Abb., 27 Tab., X, 575 S. Technisch-physikalische Monographien: Band 10. Leipzig: 
Akademische Verlagsgeselischaft Geest u. Portig K. G. 1960. Geb. DM 54,—. 


Es wird einmal als kennzeichnend fir unser Zeitalter angesehen werden, da sich zufolge des 
intensiven Fortschritts der Naturwissenschaften in rascher Folge einzelne Anwendungszweige 
zu umfangreichen, eigengesetzlichen Gebieten entwickeln, die man als Spezialtechniken bezeichnet 
und die nicht selten dank ihrer Perfektion entscheidenden Einflu8 auf andere Gebiete der Technik 
oder Forschung gewinnen. 


Die Erzeugung und Anwendung von Kondensatorentladungen — besonders dann, wenn es 
sich um extrem groBe Strom- oder Spannungsimpulse oder um extrem rasche Anderungen handelt — 
ist ein solches Gebiet. Praktische Anwendungen sind die , Umwandlung von Kondensatorenergie 
in Réntgenblitze‘‘, die ,, Umwandlung kapazitiv gespeicherter Energie in Warme zum Schweifen 
oder zur Mikroinduktionshartung“’, die ,,Umsetzung kapazitiv gespeicherter Energie in kurz- 
dauernde sehr hohe Magnetfelder‘‘, wie sie in der Kernphysik oder zur Herstellung von Dauer- 
magneten erforderlich sind, die ,,Umwandlung kapazitiv gespeicherter Energie in akustische 
Impulse“ zur Verkehrssicherung (LuftschallstéBe) oder im Wasser, die ,,Materialbearbeitung 
durch hochfrequente Kondensatorentladungen‘‘ (Funkenerosion) und die ,,Lichtblitzerzeugung“, 
um nur die diesbeziiglich wichtigsten Schlagworte des Inhaltsverzeichnisses zu nennen. 


Allen diesen Anwendungsgebieten ist der elektrische Mechanismus der gesteuerten Konden- 
satorentladungen gemeinsam; eine einheitliche Behandlung aber ist zumindest dem Bericht- 
erstatter in der Fachliteratur bisher nicht bekannt. Einzelheiten sind selbstverstaindlich zerstreut 
im Zusammenhang mit den speziellen zum Teil sehr wichtigen Anwendungen zahlreich veroffent- 
licht. Dr.-Ing. Friingel ist mit dem ganzen Anwendungsvolumen der gesteuerten Kondensator- 
entladungen, wie es oben geschildert, praktisch und theoretisch eng verknipft. So konnte eine 
Monographie dieser Technik im besten Sinn des Wortes von ihm geschrieben werden, die vielerlei 
Bedirfnisse befriedigen wird: Forscher, die schnell verlaufende Vorgange analysieren wollen, 
Meteorologen und Verkehrsstrategen, die Lichtblitzsignale einsetzen wollen, werden sich alle 
ein gutes Bild tiber vorhandene Moglichkeiten und notwendigen Aufwand schaffen. Aber auBer 
diesen Spezialisten werden sich ebenso Studenten und alle technisch Interessierten gerne des 
Buches bedienen, um sich ein grundsatzliches Bild zu schaffen, was es auf dem Gebiete der Erzeu- 
gung und Anwendung von Kondensatorentladungen heute gibt. Zwei wichtige Eigenschaften 
des Buches wollen wir hier betonen, die wir konsequent durch das ganze Buch durchverfolgen 
konnten: 1. Bei jedem Unterabschnitt der klaren Gliederung findet sich mindest ein Abschnitt, 
der den praktischen Zweck und die zugrunde liegenden Voraussetzungen umreift. 2. Fast fir 
jeden behandelten Spezialfall sind gut verstandliche Konstruktions- und Schaltbilder gebracht. 
3. Durchgehend sind entscheidende energetische Uberlegungen klargelegt. 

AbschlieBend méchte ich dieses Buch als das Musterbeispiel einer Monographie fir eine moderne 
Spezialtechnik bezeichnen. Hiezu gehért die sachliche Beschrankung bei der Wiedergabe des 
mathematischen Riistzeugs. Der Verfasser setzt allerdings die Beherrschung der elektrotechnischen 
Grundlagen voraus, wie sie dem Stoff einer heutigen technischen Hochschule entsprechen; dadurch 
ist es ihm moéglich, das Wesentliche in jedem Fall rasch herauszuarbeiten. Die Plastik des Eindrucks 
ist durch die kraftvolle und eigenwillige Sprache meist stark unterstiitzt. Daher verzeiht man 
gerne gelegentlich die Verwendung einzelner weniger, im allgemeinen Sprachgebrauch sonst 
anders gebrauchter Ausdriicke. Vor allem betrifft dies den Ubertitel Impulstechnik. Erzeugung 
und Anwendung von Kondensatorentladungen sind nur ein Teilgebiet davon; ja, man wird in 
erster Linie bei Impulstechnik zuerst an Telegraphie hoher Geschwindigkeit, Richtfunkstrecken 
mit Impulsmodulation, Funkpeilung und Funkortung, Fernsehen, Oszillographen- oder Digital- 
Zahltechnik, kaum aber an Kondensatorentladungstechnik denken. Abgrenzung gegen diese 


Techniken sollte in der Einleitung der nachsten Auflage vorgenommen werden. 
H.von Bertele, Wien 


* Osterreichisches Ingenieur-Archiv 11, 328, 1957. 
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Begriindung der Funktionentheorie aut alten und neuen Wegen. Von L. Heffter. Zweite, wesent- 
lich verbesserte Auflage. Mit 13 Abbildungen im Text, VIII, 64 S. Berlin-Géttingen - Heidel- 
berg: Springer-Verlag. 1960. DM 19,80. 


Der Nestor der deutschen Mathematiker hat in den letzten Jahren eine Reihe von Arbeiten 
geschrieben, die einer méglichst einfachen Charakterisierung der analytischen Funktionen gelten. 
Hier wird nun eine historische Zusammenstellung dieser und anderer solcher Satze und Beweis- 
anordnungen gegeben. H. Hornich, Wien 


Konstruktive Geometrie in der Technik. Von F. Hohenberg. Zweite, neubearbeitete und erwei- 
terte Auflage. Mit 459 Textabb., IX, 319 S. Wien: Springer-Verlag. 1961. Steif geh. 5 150,—, 
DM 25,—, sfr. 25,60, $ 5,95, £ 2/2/6. Ganzln. 8 174,—, DM 29,—, sfr. 29,70, $ 6,90, £ 2/9/6. 


Dieses ausgezeichnete Lehrbuch der darstellenden und kinetischen Geometrie in ingenieur- 
maiger Sicht, dessen Qualitaéten an dieser Stelle bereits bei seinem ersten Erscheinen gewiirdigt 
worden sind (vgl. Osterr. Ing.-Archiv 10/1956, S. 416), hat sich als durchschlagender Erfolg erwie- 
sen und schon jetzt eine Neuauflage notwendig gemacht. Diese Gelegenheit wurde auBer zu einer 
kleinen Titelinderung zu einigen Erweiterungen des Inhalts beniitzt, die den Umfang nicht unbe- 
trichtlich vergréBerten. Solche Erginzungen finden sich u. a. bei der raumlichen Zentralkollineation 
und Reliefperspektive, bei der Schraubtorse und beim Frasen von Schraubnuten im darstellend- 
geometrischen Teil, ferner bei den Koppelkurven, Radlinien und Kriimmungskonstruktionen im 
kinematischen Teil. Hier wurde auch ein neuer Abschnitt tiber ,, Trochoidenmaschinen* aufgenom- 
men, als deren bekanntester Vertreter der Wankel-Motor von NSU sich gegenwartig besonderen 
Interesses erfreut. 

Das vielseitige Werk hat durch den vorgenommenen Ausbau zweifellos noch an Wert gewonnen 
und kann vor allem den jungen Ingenieurstudenten wairmstens empfohlen werden. Dem Verneh- 
men nach sind derzeit Ubersetzungen ins Schwedische und Spanische in Vorbereitung, was als 
sichtbare Anerkennung der sinnvollen und seit jeher hochstehenden Pflege der Geometrie an den 
Technischen Hochschulen Osterreichs aufgefaBt werden darf. W. Wunderlich, Wien 


Quasikonforme Abbildungen. Von Hans P. Kiinzi. Mit 35 Abb., VIII, 182 8. Ergebnisse der 
Mathematik und ihrer Grenzgebiete: Heft 26. Berlin-Gottingen-Heidelberg: Springer-Verlag. 
1960. DM 39,—. 


Bei der groBen Bedeutung der quasikonformen Abbildungen war eine zusammenfassende 
Darstellung dieses Gebietes sehr an der Zeit. Zunachst werden die konformen und die quasi- 
konformen Abbildungen (nach Grotzsch) und deren Anwendungen in der Funktionentheorie 
behandelt; fiir zwei allgemeinere K-quasikonforme Homéomorphismen wird der Aquivalenzsatz 
bewiesen. Ein weiterer Abschnitt gilt den von Teichmiller aufgegriffenen extremalen quasi- 
konformen Abbildungen. Endlich werden die pseudoanalytischen Funktionen behandelt und der 
Lavrentieffsche Fundamentalsatz ftir quasikonforme Abbildungen gegeben. 

H. Hornich, Wien 


Grundbau-Dynamik. Von H. Lorenz. Mit 302 Abb., VIII, 308 8S. Berlin-Gottingen-Heidelberg: 
Springer-Verlag. 1960. Geb. DM 46,50. 


Das Buch gliedert sich in zwei Hauptteile. Der erste, der etwas mehr als die Halfte des Umfanges 
einnimmt, bringt eine — nach den Worten des Verfassers auf die Denkweise des Bauingenieurs 
zugeschnittene — Schwingungslehre fir Systeme mit konzentrierten und mit verteilten Massen. 
Neben exakten Losungen werden auch einige der wichtigsten in der Literatur bekannt gewordenen 
Naherungsverfahren angegeben. Bei der Behandlung der Stabwerke stiitzt sich der Verfasser im 
wesentlichen auf das bekannte Buch von Hohenemser und Prager. Hier wire vielleicht eine 
Erganzung durch neuere Verfahren, wie z. B. die Methode der Ubertragungsmatrizen, vorteilhaft 
gewesen. Ein Abschnitt ist den nichtlinearen Schwingungen gewidmet. 

Die Aufstellung der Bewegungsgleichungen, der erste und wichtigste Schritt bei der Unter- 
suchung eines schwingenden Systems, wird leider nirgends, nicht einmal am einfachsten Beispiel, 
vorgefiihrt. Auch erscheint dem Besprecher das Bemithen des Verfassers, die Bezeichnungsweise 
der Statik bei der Behandlung dynamischer Probleme beizubehalten, nicht glicklich. So wird z. B. 
haufig von Lasten gesprochen, wo Massen gemeint sind. Sogar Ausdriicke wie , Eigenfrequenz 
einer Belastung“ (Seite 134) kommen vor. Der Besprecher ist nicht der Ansicht, da8 damit dem 
Leser, auch wenn er in erster Linie Statiker ist, ein Dienst erwiesen wird. Ein solches Verfahren kann 
doch nur das wirkliche Erfassen und Verstehen der Zusammenhiinge erschweren. Auch sonst 


treten gelegentlich recht unscharfe Formulierungen auf, wie etwa bei der Besprechung der Ritz- 
schen Methode. 
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Der Wert des Buches liegt nach Ansicht des Besprechers wohl vor allem im zweiten Teil, der 
sich mit Baugrunddynamik beschiftigt. Hier wird, ausgehend von E. Reissners Berliner Disser- 
tation, erstmalig ein zusammenfassender Uberblick tiber die Ansitze gegeben, mit denen man bis 
jetat dem so auSerordentlich schwierigen Problem des dynamischen Verhaltens des Baugrundes 
beizukommen versuchte. Freilich bringt der Verfasser im wesentlichen nur Resultate; er sichtet 
und vergleicht aber kritisch. Besprochen werden die Arbeiten von E. Reissner (die Kapitel- 
uberschrift ist irreftthrend; es handelt sich dabei nicht um eine Punktmasse auf dem Halbraum), 
ferner die Ehlerssche Naiherung und dann erfreulicherweise die neuere, im deutschen Sprachraum 
kaum bekanntgewordene Arbeit von Novak. Ein weiterer Abschnitt ist einer Ubersicht tiber die 
eimschlagigen Arbeiten des Shell-Laboratoriums in Amsterdam gewidmet. Zusammenfassend 
zeigt sich freilich, da8 wir von einer wirklich befriedigenden Theorie noch weit entfernt sind. 

Zwei Kapitel iber Baugrunduntersuchung und iiber Bodenverbesserung durch Schwingungs- 
verdichtung sowie ein kurzer Abschnitt ttber Schwingrammung beschlieBen die baugrunddynami- 
schen Betrachtungen. 

Die Literatur tiber Fundamentschwingungen und Bodendynamik ist nicht allzu reichhaltig. 
Das vorliegende Buch bildet daher, vor allem in seiner zweiten Halfte, eine willkommene Erganzung 
und wird von einem kritischen Leser mit Gewinn studiert werden. H. Parkus, Wien 


Mechanik der festen Kérper. Von H. Parkus. Mit 191 Textabb., VIII, 264 8. Wien: Springer-Verlag. 
1960. Steif geh. S 156,—, DM 26,—, sfr. 26,60, $ 6,20, £ 2/4/6. Ganzleinen S 174,—, DM 29,—, 
sir. 29,70, $ 6,90, £ 2/9/6. 


Die in diesem Buch niedergelegten Vorlesungen iiber Technische Mechanik fiir Studierende des 
Maschinenbaues, der Elektrotechnik und der Technischen Physik sind schon insofern beachtlich, 
als sie vom traditionellen Aufbau der Mechanik Abstand nehmen, indem bereits die Grundlagen 
auf die Kontinuumsmechanik hin strukturiert werden. Damit erreicht der Verfasser eine gedank- 
liche Okonomie, die es gestattet, auf knappen 260 Seiten die Mechanik des starren und elastischen 
Korpers in sehr umfassender Art zur Darstellung zu bringen. 

Das Werk beginnt in unorthodoxer Weise mit der Kinematik, woran die Behandlung der Krifte 
und Kraftegruppen und die Massengeometrie anschlieBen. Die Grundlagen der Dynamik werden 
gleich fiir das Kontinuum formuliert; es folgen Schwerpunkt- und Drallsatz, der Arbeits- und Energie- 
satz, das D’Alembertsche Prinzip und die Lagrangeschen Gleichungen. Die folgenden Abschnitte 
bringen die Elastizitatstheorie und ihre Anwendung auf stabformige Korper, Kreisplatten und 
Rotationsschalen. Es folgen die Satze ttber Formanderungsarbeit und ein Einblick in die Theorie 
der Warmespannungen. Das Werk schlieBt mit der Behandlung der Stabilitatskriterien, dem 
Verfahren von Ritz und einer Theorie der StoBvorgange. Im Anhang werden die wichtigsten Formeln 
der Vektor-Algebra und -Analysis angefthrt. 

Damit ist ein den modernen Anforderungen der Technik voll gerecht werdendes Kompendium 
der Technischen Mechanik der festen Korper geschaffen, das jedem Kenner gewii Freude bereiten 
wird. Man muB sich allerdings fragen, ob der padagogische Optimismus des Verfassers gerechtfer- 
tigt ist, der bereits dem Anfanger dieses Niveau der Darbietung zumutet. Wenn zum Beispiel die 
dynamische Grundgleichung gleich im Zusammenhang mit den Spannungen am Substanzelement 
formuliert wird, dirfte dies die Erfassung und Erarbeitung der Grundlagen nicht erleichtern. Kine 
vorangehende, elementare Vermittlung des physikalischen Unterbaues wirde von den Horern 
wahrscheinlich dankbar aufgenommen werden. Der Fortgeschrittene kann natitrlich eher auf solche 
Konzessionen an sein Verstandnisvermégen verzichten. G. Heimrich, Wien 


Uber die totale Stabilitit erzwungener Reibungsschwingungen. Von RF. Reissig. Mit 10 Abb., 28 S. 
Berlin: Akademie Verlag. 1959. Brosch. DM 5,—. 
Abhandlungen der Deutschen Akademie der Wissenschaften Berlin, Klasse fir Mathematik, 
Physik und Technik, 1959, Nr. 1. 


Der Verfasser untersucht die erzwungenen Schwingungen eines Schwingers mit nichtlinearen 
Riickstell- und Dampfungskraften und Coulombscher Reibung. Vor allem wird die Stabilitat 
untersucht. Es ergibt sich bei jeder Bewegung unter gewissen einschrankenden Bedingungen 
totale Stabilitit. Fir gewisse Stérungen ist sie sogar asymptotisch. Der Stabilitatsbeweis weicht 
von den bekannten Ljapunowschen Untersuchungen ab, weil die hier gebildete Ljapunowsche 
Funktion nicht den Bedingungen fir totale Stabilitat gentigt. (Die Funktion andert sich sprunghaft, 
wenn die Geschwindigkeit ihr Vorzeichen wechselt.) Weiters wird ein Schwinger mit linearer 
Federcharakteristik, der eine Storung in Form einer Anderung des Reibungskoeffizienten erfahrt, 
untersucht. Mit diesem Modell kann auch der Unterschied zwischen Haft- und Gleitreibung erfaBt 
werden. Ein Zahlenbeispiel hiezu wird gebracht. Bemerkenswert ist, dal dies die erste geschlossene 
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Behandlung eines derartigen Problems darstellt. Von praktischer Bedeutung ware noch die Beant- 
wortung der Frage, inwieweit mit Hilfe des Verfahrens von Krylow-Bogoljubow brauchbare 
Naherungslosungen gefunden werden konnen. W. Peschka, Wien 


Ingenieur-Mathematik. Von R. Sauer. Erster Band: Differential- und Integralrechnung. Mit 
178 Abb., VIII, 304 S. Berlin-Géttingen-Heidelberg: Springer-Verlag. 1959. Geb. DM 24,—. 


Der vorliegende erste Band des geplanten zweibandigen Werkes umfa8t im wesentlichen den 
Stoff fiir das erste Jahr der an unseren Technischen Hochschulen gehaltenen mathematischen 
Grundvorlesung. Da sich die ,,Ingenieur-Mathematik‘’ vornehmlich an angehende Techniker 
wendet, wurde — bei aller mathematischen Strenge — auch auf Anschaulichkeit gréBter Wert 
gelegt. Deshalb sind die grundlegenden Begriffe, auf deren Herausarbeitung mit besonderer Sorgfalt 
eingegangen wurde, auch stets an konkreten Beispielen erlautert und die im Text angefiihrten 
Lehrsitze diskutiert, die meisten ihrer Beweise aber in einem Anhang am Schlu8 des Buches 
zusammengestellt. 

Was den Inhalt betrifft, so gliedert er sich in zwei Kapitel; Differential- und Integralrechnung 
fiir Funktionen von einer und mehreren Veranderlichen. Hierin sind aber auch enthalten: Vektor- 
algebra, lineare Gleichungen, lineare analytische Geometrie, Differentialgeometrie der Kurven 
und Flichen im Raum, Nomographie sowie graphische und numerische Methoden der Analysis. 
Insbesondere ist die Darstellung der numerischen Methoden — auch an Hand einfachster Bei- 
spiele — auf die Verwendung moderner Rechenmaschinen zugeschnitten. 

Der didaktisch sehr zweckmaBige Aufbau und die vielen durchgerechneten Beispiele sowie 
die groBe Anzahl instruktiver Abbildungen werden sicher dazu beitragen, diesem Buch einen 
weiten Leserkreis zu sichern. H. Fieber, Wien 


Ingenieur-Mathematik. Von R. Sawer. Zweiter Band: Differentialgleichungen und Funktionentheorie. 
Mit 93 Abb., XI, 180 S. Berlin-Gottingen-Heidelberg: Springer-Verlag. 1961. Ganzleinen 
DM 18,—. 


Der Stoff der Mathematik-Vorlesungen des zweiten Studienjahres an unseren technischen 
Hochschulen wird hier — ahnlich wie beim ersten Band des Werkes — kurz, aber klar und einprag- 
sam gebracht, wobei die Theorie nur auf das zum Verstandnis Notwendige beschrankt und einige 
Beweise wieder auf den Anhang verwiesen werden; vielfach werden numerische Methoden vorge- 
fiihrt. Nach einem kurzen Abschnitt tber Vektoranalysis werden gewohnliche und am Beispiel der 
Wellengleichung auch partielle Differentialgleichungen behandelt. Dabei werden auch die Fourier- 
schen Reihen eingeftthrt. Die zweite Halfte des Buches gehort einer Einftthrung in die Funktionen- 
theorie, vornehmlich in die konforme Abbildung. — Die beiden Bande der ,, Ingenieur-Mathematik“ 
kommen einem wirklichen Bediirfnis an den technischen Hochschulen entgegen und werden sicher 
starke Verbreitung finden. H. Hornich, Wien 


Grenzschicht-Theorie. Von H. Schlichting. Dritte erweiterte und neubearbeitete Auflage. Mit 
374 Textabb., XV, 603 S. Karlsruhe: Verlag G. Braun. 1958. 


Das bekannte Standardwerk liegt nun in dritter Auflage vor. Im Gegensatz zur zweiten Auf- 
lage handelt es sich nicht um einen unveradnderten Nachdruck, sondern dem Verfasser ist es gelun- 
gen, alles Wesentliche aus der tiberreichen, seit der Herausgabe der ersten Auflage erschienenen 
Literatur tiber dieses Fachgebiet (etwa hundert Arbeiten pro Jahr) in den Text einzuarbeiten. 
Dabei wurde die Kapiteleinteilung des Buches nicht geandert, die einzelnen Abschnitte jedoch 
vielfach erganzt und erweitert. 

Als Erganzung hinzugekommen sind insbesondere neue exakte Lésungen fiir laminare Grenz- 
schichten und weitere Ausfiihrungen tiber das praktisch wichtige Gebiet der Grenzschichtbeein- 
flussung. Das sehr moderne Thema der kompressiblen laminaren Grenzschichten wurde ebenfalls 
erheblich erweitert. In dem Abschnitt tiber Turbulenz-Entstehung wurden viele Ergiinzungen in 
den Anwendungen der Stabilitatstheorie hinzugefiigt. Auch dem Einflu8 des Warmeiiberganges 
und der Kompressibilitat auf den Umschlag der laminaren in die turbulente Strémung wurde nach- 
gegangen. Die turbulenten Grenzschichten mit Druckanderung lings der Stromung finden nunmehr 
eine dem gegenwartigen Stand der Forschung entsprechende Darstellung. Uber den Einflu8 der 
Kompressibilitat auf die turbulente Grenzschicht finden sich nur wenige Angaben, da auf diesem 
Gebiet klare und eindeutige Ergebnisse noch sparlich anzutreffen sind. 

Die Meisterschaft des Verfassers, ein kompliziertes und komplexes Stoffgebiet klar und wber- 
sichtlich zu gliedern und dem Leser leicht faflich und anschaulich darzubieten, zeigt sich auch in 
dieser Auflage. Die auferst reichhaltige Literaturibersicht vermittelt einen bequemen Zugang 
zu dem weitverstreuten Schrifttum, was sowohl der Forscher als auch der praktische Ingenieur 
mit Dank quittieren werden. G. Heinrich, Wien 
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Der durehlaufende Bogentriger auf elastischen Stiitzen mit und ohne Versteifungstriger. Von 


W. Stampf. Mit 206 Textabb., VIII, 196 S. Berlin-Gottingen-Heidelberg: Springer-Verlag. 
1960. Geb. DM 37,50. 


Far das hochgradig statisch unbestimmte System des durchlaufenden Bogentragers auf biege- 
elastischen Stiitzen, dessen Berechnung nach dem KraftgroSenverfahren groRen Arbeitsaufwand 
erfordert, wird hier die Deformationsmethode verwendet. Durch die Einfithrung der Horizontal- 
verschiebungen und Verdrehungen der Bogenknotenpunkte als Unbekannte wird nicht nur die 
Anzahl der zu lésenden Gleichungen reduziert, sondern auch die statische Behandlung des Systems 
anschaulich und iibersichtlich gehalten. Die lotrechten Knotenverschiebungen durch Zusammen- 
driicken der Pfeiler werden als sehr klein nicht beriicksichtigt, da vor allem auch nur die gegenseiti- 
gen Verschiebungen von Nachbarknoten einen Einflu8 auf das Kraftespiel ausiiben. So ergeben 
sich fir ein System mit » Knoten 2 Gleichgewichtsbedingungen als Bestimmungsgleichungen 
fir die unbekannten Verschiebungen und Verdrehungen. Zur Auflésung dieser Gleichungen, 
deren Anzahl bei langen Bogenreihen auch hier stark ansteigt, wird eine abgekiirzte Methode ange- 
geben, die das Gesamtlgeichungssystem in zwei dreigliedrige Systeme aufspaltet und durch auf- 
einanderfolgende Korrekturen lést. 

Im Teil A des Buches werden zunachst fiir die Reihe von ,,freien‘‘ Bogen auf biegeelastischen 
Pfeilern die fir die Deformationsmethode nétigen Grundwerte abgeleitet. Neben dem allgemeinen 
Belastungsfall des symmetrischen Bogens sowie dessen Eigengewichtsbelastung sind auch Tempe- 
raturanderung und Schwinden des Betons in die Betrachtung mit einbezogen. Zur genauen Erfas- 
sung der ungiinstigsten Verkehrsbelastung wird die Ermittlung der EinfluBlinien angegeben. 

Der Teil B erweitert diese Methode auf durchlaufende Bogen mit aufgestinderten Versteifungs- 
tragern, wobei das System auf dieselbe Form der Bedingungsgleichungen zuriickgefiihrt wird wie 
das Grundtragssystem des freien Bogens. Dadurch bleibt die Behandlung auch dieses an sich kom- 
plizierten Systems einfach und anschaulich. 

Im I. Abschnitt des Teiles C geht der Verfasser vom Briickenbau zum Gebiet des Hochbaues 
uber, wo solche und ahnliche Systeme als Bogen- oder Rahmenreihen, z. B. bei Hallenbauten, 
vielfach Anwendung finden. 

Im II. Abschnitt stellt er sich die Aufgabe, zum Zweck der Vorberechnung und Uberschlags- 
dimensionierung solcher Tragwerke eine rationelle Naherungsmethode zu entwickeln. Unter ver- 
einfachten Voraussetzungen werden explizite Formeln und Tabellen angegeben, die fiir den ent- 
werfenden Ingenieur besonders wertvoll sein werden. Im Anhang wird noch auf ein besonderes 
Auflosungsverfahren der Knotenpunktsgleichungen mit Hilfe von ,,Ubergangszahlen“ einge- 
gangen. 

Die zu jedem Abschnitt durchgefiihrten ausfihrlichen Zahlenbeispiele erlautern deutlich die 
jeweils dargelegten theoretischen Entwicklungen und erleichtern so die Aufnahme und Anwendung 
der hier entwickelten Methode. Die sauberen und klaren Zeichnungen sowie der deutliche Satz 
aller Zeichen, Formeln und Zahlen — ein Verdienst des Verlages — sind besonders hervorzuheben 
und tragen zum leichteren Verstandnis bei. 

Das vorliegende Buch wird sowohl vom Studierenden, der sich mit Sonderfragen der Baustatik 
beschaftigt, als auch von dem auf diesem Gebiet tatigen Ingenieur mit Freude und Nutzen studiert 
und angewendet werden. FE. Resinger, Graz 


Repertorium und Ubungsbuch der Technischen Mechanik. Von J. Szabé. Mit 254 Abb., VII, 273 S. 
Berlin-Géttingen-Heidelberg: Springer-Verlag. 1960. Geb. DM 24,—. 


Der Verfasser setzt sich das Ziel, nebst einer Sammlung von Ubungsbeispielen auch die wesent- 
lichen Satze und Formeln der Technischen Mechanik knapp zur Darstellung zu bringen. Fir die 
Erstellung der Ubungsbeispiele wird jeder, der in diesem Fach Priifungen abzuhalten hat, dankbar 
sein. Die Beispiele umfassen alle wesentlichen Teile der elementaren, angewandten Mechanik und 
sie sind mit Geschick und padagogischem Spirsinn ausgesucht. 

Bei der Rekapitulation der Lehrsitze sind hingegen die Formulierungen nicht immer ganz 
gegliickt; so wird z. B. der Inhalt der Galileischen Relativitatsprinzipes durch die Aussage zu 
erklaren versucht, da& es gleichgiiltig sei, ob sich das System in Ruhe oder in gleichformiger Bewe- 
gung befindet. ; ! 

Das Beispiel itber das Abfangen einer fallenden Masse durch einen elastischen Faden ware 
besser bei den Schwingungs- und nicht bei den StoBaufgaben einzureihen, da hier keine echte 
(idealisierbare) StoSkraft vorliegt. 

Diese und ahnliche kleine Unzulanglichkeiten, die leicht in den kommenden Auflagen beseitigt 
werden konnen, beeintrachtigen nicht den Wert des Repertoriums, das die vorhandenen Ubungs- 
biicher in wertvoller Weise erganzt. G. Heinrich, Wien 
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Elements of Materials Science. Von Lawrence H. Van Vlack. Mit 498 Abb., 528 S. Reading, Mass. — 
London: Addison-Wesley Publishing Company, Inc. 1959. $ 8,50, £ -/64/-. 


Der Inhalt dieser ,,Grundlagen der Werkstoffkunde‘‘, der eine Einfiihrung fiir Studierende 
an Technischen Hochschulen geben soll, besitzt einen kennzeichnenden Aufbau, weil in neun von 
insgesamt fiinfzehn Abschnitten der Zusammenhang zwischen Werkstoffeigenschaften und Werk- 
stoffstruktur systematisch behandelt wird. Ausgehend von den interatomaren Bindungskraften 
und den Kristall- bzw. Molekilstrukturen werden einphasige Werkstoffe (Metalle, keramische 
Werkstoffe und organische Werkstoffe), mehrphasige Werkstoffe (Zustandslehre) und die Zusam- 
menhiinge zwischen .Werkstoffeigenschaften und Mikro- bzw. Makrostruktur dargestellt. Die 
folgenden Abschnitte sind dem Werkstoffverhalten bei deren Verwendung gewidmet, ihren 
mechanisch-technologischen Eigenschaften, ihrem Verhalten gegentiber thermischer Beanspru- 
chung, bei Korrosion, in elektromagnetischen Feldern und gegentiber Korpuskular- sowie Wellen- 
strahlung. Am Schluf eines jeden Abschnittes findet sich ein Schrifttumsverzeichnis fir ein 
weitergehendes Studium sowie eine Aufgabensammlung. Das Buch besitzt einen Anhang mit 
Zustandsdiagrammen, Umwandlungsschaubildern und Tabellen. Der Autor hat auf verhaltnis- 
mifig kleinem Raum erstaunlich viel aus dem Stoffgebiet der Werkstoffkunde in klarer und 
iiberaus anschaulicher Weise dargestellt. Der Text ist durch zahlreiche instruktive Abbildungen 
erginzt und beriicksichtigt auch die neuere Entwicklung, wobei er sich auf das Grundsatzliche 
beschrankt. K. Lotsch, Wien 


Algebra I. Von B. L. van der Waerden. Finfte Auflage der modernen Algebra. VIII, 292 S. 
(Grundlehren der mathematischen Wissenschaften: Band 33.) Berlin-Gottingen-Heidelberg: 
Springer-Verlag. 1960. Ganzleinen DM 22,—. 


Eine griindliche Kenntnis algebraischer Methoden ist heute nicht nur fiir die reine Mathematik, 
sondern auch fiir viele Anwendungen in Physik und Technik unerlaBlich. Das vorliegende sehr 
bekannte Werk gibt eine ganz hervorragende Einfiihrung in die abstrakte Algebra und ihre Denk- 
weise und bringt nach Definition der grundlegenden Begriffe Menge, Gruppe, Ring, Korper usw. 
die Steinitzsche Kérpertheorie, Galois’sche Gleichungstheorie fiir kommutative Grundkorper, 
algebraische Theorie der reellen Kérper und Bewegungstheorie. Daneben werden aber auch Ergeb- 
nisse der klassischen Algebra, wie Resultantentheorie, Partialbruchzerlegung und Interpolations- 
formeln, beriicksichtigt. Die vorbildlich klare Darstellung und Gliederung macht das Buch als 
Nachschlagewerk und als Lehrbuch bestens geeignet und auch fiir den Anfanger leicht lesbar. 

Die vorliegende Auflage ist gegeniiber der vierten unverandert. A. Florian, Wien 


Mechanik. Von H. Ziegler. Band I: Statik der starren und flissigen Korper sowie Festigkeitslehre. 
Dritte, neubearbeitete Auflage. Mit zahlreichen Textabb., 244 8. Basel und Stuttgart: Birk- 
hauser-Verlag. Geb. sfr. 28,50. 


Die vorliegende dritte Auflage des ersten Bandes des bekannten Lehrbuches behandelt die 
elementare Statik und Festigkeitslehre. Begonnen wird in traditioneller Weise mit dem starren 
Korper und der Reduktion von Kraftegruppen. Vektoren sind dabei — eine Konzession an den angel- 
sachsischen Gebrauch — mit fetten Buchstaben bezeichnet. Dann werden die Begriffe Kraftfeld, 
Arbeit und Potential eingefiihrt, nach Ansicht des Besprechers vielleicht nicht immer mit der 
wunschenswerten Klarheit; man vergleiche etwa die beiden Ausdriicke (13,13) und (13,17). Es 
folgen dann der Spannungsbegriff und der Ubergang zur Hydrostatik. Der Rest des Buches ist 
der elementaren Festigkeitslehre gewidmet. Eine Ausnahme machen die Abschnitte 18 bis 20, in 
denen der réumliche Spannungs- und Verzerrungszustand (unter Beschrankung auf kleine Verfor- 
mungen und Weglassung umfangreicher Beweise) behandelt wird und an die dann als sehr begrii- 
Benswerte Neueinschaltung ein Abschnitt iber Flie8- und Bruchbedingungen anschlieBt. Damit 
kann auch die Traglastberechnung an einfachen Beispielen vorgefihrt werden. 

Das bewahrte Buch wird den Studierenden sicher auch weiterhin gute Dienste leisten 

H. Parkus, Wien 


Higentiimer, Herausgeber und Verleger: Springer-Verlag in Wien I, Mélkerbastei 5. — Fiir den Inhalt verant- 
wortlich: Prof. Dr. H. Parkus, Wien IV, Technische Hochschule, Karlsplatz 13. — Druck: Steyrermtihl, Wien VI, 
GumpendorferstraBe 40 —44. 
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In English Language / In englischer Sprache 


From the Preface 


This book was prepared by distinguished experts in the field of computers for the non-expert 
who does not intend to become an expert but who wants to acquire a general understanding of 
the problems and solutions which he needs to effectively perform the work in his proper field. 

So we address this book to the engineer and scientist, who want to know the performance of com- 
puters, as well as to the manager, who is mainly concerned with their economic aspect. 

The technical development of computing machines is going on at a rapid rate. Any detailed des- 
cription of computer components would therefore become obsolete within a few years. The general 
principles underlying their operation will probably remain unchanged for a longer time. Accord- 
ingly, this book puts main emphasis on principles and methods rather than on engineering details. 


Kernreaktortheorie 


Eine Einfiihrung 


Von ; 
S. Glasstone und M. C. Edlund 
United States Atomic Hnergy Commission Oak Ridge National Laboratory 
Ins Deutsche tibersetzt und bearbeitet von 
Dr. W. Glaser und Dr. H. Grimm 
weiland o. Professor Reaktor-Interessengemeinschaft 
an der Technischen Hochschule Wien Reaktorabteilung der SGP, Wien 


Mit 82 Textabbildungen. X, 341 Seiten. Gr.-8°. 1961 
Ganzleinen § 239.—, DM 38.—, sfr. 40.90, $ 9.50 
Der ,,Glasstone-Edlund“‘, wie man dieses Buch in Fachkreisen nennt, ist bereits ein klassisches Werk 


geworden. Von der 1952 erschienenen englischen Originalausgabe liegt heute der 7. Neudruck vor. 


Durch das Erscheinen dieses Buches sind die friedlichen Aspekte der Kernenergie — bis dahin 
uberlagert von militaérischen Restriktionen — zum erstenmal in wissenschaftlich klarer und erstaun- 
lich umfassender Weise freigelegt worden. Seit langem ist dieses Werk die wichtigste, allgemein 
zugangliche Quelle fir die Fragen der Reaktortheorie. An ihm wurde eine ganze Generation von 
Reaktorphysikern und Kernenergetikern erzogen. . 


In der Zwischenzeit ist die Kerntechnik von der Theorie und vom Experiment fortgeschritten zum 


Gro8-Kernkraftwerk. Trotz dieser stiirmischen Entwicklung hat der ., Glasstone-Edlund*‘ seinen 


Rang als eines der besten Lehrbiicher gehalten, das fiir den fachlich interessierten Hochschiler, _ 
den Wissenschaftler, den Ingenieur in der Entwurfs- und Konstruktionsarbeit unentbehrlich ist. 
Die vorliegende deutsche Ausgabe soll das Buch dem rasch wachsenden Interessentenkreis in den 


deutschsprachigen Lindern naherbringen. Im Einverstandnis mit den Autoren der Originalausgabe 


wurden in der deutschen Bearbeitung neuere Entwicklungen, vor allem durch weiterfiihrende 


Literaturhinweise, beriicksichtigt. 
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